New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 72407 =27 (0.610)4 (0.090) 4/56 5.1 9.1% coding (446/1701 nt) araB ribulokinase
?REL606 3650765 = 53 (0.960)intergenic (+46/+141) cspA/hokA major cold shock protein/toxic polypeptide, small
Rejected: Position hash score below cutoff.

AAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCG  >  REL606/3650730‑3650800
                                   |                                   
aaaGCACAGAATCTAAGATCCCTGCCATTTGGCggg                                     >  1:160644/1‑36 (MQ=255)
aaaGCACAGAATCTAAGATCCCTGCCATTTGGCggg                                     >  2:941518/1‑36 (MQ=255)
 aaacacaGAATCTAAGATCCCTGCCATTTGGCgggg                                    <  2:3500310/33‑1 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                    >  2:3569126/1‑36 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                    >  1:3037861/1‑36 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                    <  2:3030203/36‑1 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                    <  1:2029564/36‑1 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                    <  2:2758643/36‑1 (MQ=255)
 aaGCACAGAATCTAAGAGCCCTGCCATTTGGCgggg                                    <  2:2926292/36‑1 (MQ=255)
  aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                   >  2:992280/1‑36 (MQ=255)
  aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                   >  1:2231190/1‑36 (MQ=255)
   gCACAGTATCTATGATCCCTGCCATTTGGCGGGGAt                                  <  1:432610/36‑1 (MQ=38)
   gCACAGAGTCTAAGATCGCTGCCATTTGGCGGGGAt                                  <  2:1863790/36‑1 (MQ=38)
    cacaGAATCTAAGATCCCTGCCATTTGGCGGGGAtt                                 >  2:399208/1‑36 (MQ=255)
    cacaGAATCTAAGACCCCTGCCATTTGGCGGGGAtt                                 <  2:2848463/36‑1 (MQ=255)
    cacaGAACCTAAGATCCCTGCCATTTGGCGGGGAtt                                 <  2:3069094/36‑1 (MQ=255)
     acaTAATCTAAGATCCCTGCCATTTGGCGGGGAttt                                <  2:3174518/36‑1 (MQ=255)
     acaGAAGCTAAGATCCCTGCCATTTGGCGGGGAttt                                <  2:2560154/36‑1 (MQ=255)
      caGAATCTAAGATCCCTGCCATTTGGCGGGGAtttt                               <  2:2635444/36‑1 (MQ=255)
      caGAATCTAAGATCCCTGCCATTTGGCGGGGAtttt                               <  1:734242/36‑1 (MQ=255)
      caGAATCTAAGACCCCTGCCATTTGGCGGGGAtttt                               <  1:453350/36‑1 (MQ=255)
       aGAATCTAAGATCCCTGCCATTTGGCGGGGAttttt                              <  1:3223623/36‑1 (MQ=255)
       aGAATCTAAGATCCCTGCCATTTGGCGGGGAtattt                              >  1:3536707/1‑36 (MQ=255)
        gAATCTAAGATCCCTGCCATTTGGCGGGGAtttttt                             >  2:3699973/1‑36 (MQ=255)
         aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                            <  1:256819/36‑1 (MQ=255)
         aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                            >  1:2103776/1‑36 (MQ=255)
           tCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAt                          <  1:1600616/36‑1 (MQ=255)
            cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAtt                         <  2:583916/36‑1 (MQ=255)
             tAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                        <  2:1446709/36‑1 (MQ=255)
              aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                       <  1:1757348/36‑1 (MQ=255)
               aGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGt                      >  1:654258/1‑36 (MQ=255)
               aGAACCCGGCCATTTGGCGGGGATTTGTTTATTTGt                      <  2:3393164/36‑1 (MQ=25)
                gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                     <  2:623166/36‑1 (MQ=255)
                 aTCCCGGCCATTTGGCGGGGATTTTTTTATTTGttt                    <  1:3723570/36‑1 (MQ=255)
                  tCCCTGCCATTTGGCGGGGATTTTTTTATTTGtttt                   >  2:3435298/1‑36 (MQ=255)
                   cccTGCCATTTGGCGGGTATTTTTTTATTTTTTTTc                  >  1:1244419/1‑36 (MQ=37)
                    ccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                 <  1:3304454/36‑1 (MQ=255)
                    ccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                 <  1:3013070/36‑1 (MQ=255)
                      tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAgg               <  1:1264075/36‑1 (MQ=255)
                      tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAgg               <  1:2305609/36‑1 (MQ=255)
                      tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAgg               <  2:2833752/36‑1 (MQ=255)
                        ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaa             <  1:3534400/36‑1 (MQ=255)
                         cttttGGCGGGGATTTTTTTATTTGTTTTCAGGaaa            <  1:518049/34‑1 (MQ=255)
                         cATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa            <  1:2027269/36‑1 (MQ=255)
                         cATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa            >  2:971304/1‑36 (MQ=255)
                         cATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa            >  1:3318880/1‑36 (MQ=255)
                          aTTTTGCGGGGATTTTTTTATTTGTGTTCAGGaaat           <  2:703784/36‑1 (MQ=37)
                          aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaat           >  2:3675244/1‑36 (MQ=255)
                           tgtGGCGGGGATTTTTTTATTTGTTTTCAGGaaata          <  1:1229913/34‑1 (MQ=255)
                            ttGGCGGGGATTTTTTTATTTGTTTTCAGGaaataa         >  2:2577816/1‑36 (MQ=255)
                             tGGCGGAGATTTTTTTATTTGTTTTCAGGaaataaa        <  1:288296/36‑1 (MQ=255)
                               gCGGGGATTTTTTTATTTGTTTTCAGGAaataaata      >  2:13804/1‑36 (MQ=255)
                                cGGGGATTTTTTTATTTGTTTTCTGGaaataaatta     >  1:232259/1‑34 (MQ=25)
                                cGGGGATTTTTTTATTTGTTTTCAGGAAataaataa     <  2:616541/36‑1 (MQ=255)
                                  gggATTTTTTTATTTGTTTTCAGGAAATAAATAatc   <  2:254336/36‑1 (MQ=255)
                                   ggATTTTTTTATTTGTTTTCAGGAAATAAATAatcg  >  1:2624313/1‑36 (MQ=255)
                                   |                                   
AAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCG  >  REL606/3650730‑3650800

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG