New junction evidence
  seq id position reads (cov) reads (cov) score skew annotation gene product
* ? REL606 72407 =24 (0.510)6 (0.130) 4/60 5.3 coding (446/1701 nt) araB ribulokinase
?REL606 3944634 = 38 (0.690)intergenic (+23/‑124) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity
Rejected: Position hash score below cutoff.

CTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCC  >  REL606/3944599‑3944668
                                   |                                  
ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                                    <  1:2413784/36‑1 (MQ=255)
ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                                    <  1:4240207/36‑1 (MQ=255)
   tcGGCGACCAATTTTGAGCCCCGCTTCGGCGGGGtt                                 <  1:5085411/36‑1 (MQ=21)
    cGCCGACCAATTTTGAACCCCGCTTCGGCGGGGttt                                <  1:582895/36‑1 (MQ=255)
     gCCGACCAATTTTGAACCCCGCTTCGGCGGGGtttt                               >  1:4240092/1‑36 (MQ=255)
     gCCGACCAATTTTGAACCCCGCTTCGGCGGGGgttt                               >  1:6779967/1‑36 (MQ=255)
      ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                              <  1:1422353/36‑1 (MQ=255)
      ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                              <  1:1890278/36‑1 (MQ=255)
        gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                            >  1:7342553/1‑36 (MQ=255)
        gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                            >  1:6044326/1‑36 (MQ=255)
        gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                            <  1:3786579/36‑1 (MQ=255)
        gACCAATTTTGAACCCCGCTTCGGCGGGGTTTCTTg                            >  1:5720689/1‑36 (MQ=255)
          ccAATTTTGAACCCCGCTTCGGCGTGGTTTTTTGtt                          <  1:3800389/36‑1 (MQ=255)
           cAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGttt                         >  1:7438438/1‑36 (MQ=255)
            aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                        <  1:906144/36‑1 (MQ=255)
            aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                        >  1:4847695/1‑36 (MQ=255)
               tttGAACCCCGCTTCGGCGGGGTTTTTTGttttttt                     >  1:3247583/1‑33 (MQ=255)
               tttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtg                     >  1:4276127/1‑36 (MQ=255)
                ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt                    <  1:4824476/36‑1 (MQ=255)
                  gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc                  <  1:238937/36‑1 (MQ=255)
                  gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc                  <  1:2959295/36‑1 (MQ=255)
                     ccccGCTTAGGCGGGGTTTTTTGTTTTCTGCGCAtt               >  1:3190502/1‑36 (MQ=18)
                      cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt              <  1:4891557/36‑1 (MQ=255)
                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCTTTTCg            >  1:3321398/1‑36 (MQ=255)
                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAtttag            >  1:2719336/1‑34 (MQ=255)
                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCg            <  1:5444578/36‑1 (MQ=255)
                          cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc          >  1:6234844/1‑36 (MQ=255)
                           ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa         <  1:1837966/36‑1 (MQ=255)
                           ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa         >  1:5518855/1‑36 (MQ=255)
                           ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa         <  1:2655250/36‑1 (MQ=255)
                           ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa         >  1:5817332/1‑36 (MQ=255)
                             cggcggGGTTTTTTGTTTTCTGTGCATTTCTTCACt       >  1:4178239/1‑36 (MQ=255)
                              ggcggGGTTTTTTGTTTTCTGTGCATTTCGTCACtt      <  1:557686/36‑1 (MQ=255)
                                cggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTc    >  1:6504966/1‑36 (MQ=255)
                                cggGGTTTTTTGTTTTCTGTGCATTTAGTCCCTTTc    >  1:5774838/1‑36 (MQ=21)
                                 ggggTTTTTTGTTTTCTGTGCATTTCGTCACttttc   >  1:1580103/1‑34 (MQ=255)
                                 ggggTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc   <  1:3485908/36‑1 (MQ=255)
                                  gggTTTTTTGTTTTATGTGCATTTCGTCACTTTccc  >  1:6273621/1‑36 (MQ=255)
                                   |                                  
CTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCC  >  REL606/3944599‑3944668

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG