New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 381306239 (0.850)8 (0.180) 4/58 5.2 16.7% intergenic (‑143/‑791) ivbL/emrD ilvB operon leader peptide/multidrug efflux system protein
?REL606 = 4157631 41 (0.740)coding (278/546 nt) nusG transcription antitermination protein NusG
Rejected: Position hash score below cutoff.

TTGTACTGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAC  >  REL606/3813028‑3813097
                                  |                                   
ttGTACTGAATTGCAGATAACAAAAAACCCCGCCgg                                    >  2:2784754/1‑36 (MQ=255)
  gTACTGAATTGCAGATAACAAAAAACCCCGCTgggg                                  >  1:1775477/1‑34 (MQ=21)
  gTACTGAATTGCAGATAACAAAAAACCCCGCCGGAg                                  <  2:1546924/36‑1 (MQ=255)
    aCTGAATTGCAGATAACAAAAAACCCCGCCGGAggg                                >  1:2125446/1‑34 (MQ=255)
    aCTGAATTGCAGATAACAAAAAACCCCGCCGGAGCg                                <  2:3261625/36‑1 (MQ=255)
     cTGAATTGCAGATAACAAAAAACCCCGCCGGGGCGa                               >  2:1214421/1‑36 (MQ=255)
       gAATTACAGATAACAAAAAACCCCGCCGGAGCGAgg                             <  1:379808/36‑1 (MQ=255)
        aaTTGCAGATAACAAAAAACCCCGCCGGAGCGAGGt                            >  1:200920/1‑36 (MQ=40)
         aTTGCAGATAACAAAAAACCCCGCCGGAGCAAAGtt                           >  1:2891109/1‑36 (MQ=14)
          ttGCAGATAACAAAAAACCCCGCCGGAGCGAGGttt                          <  2:1364394/36‑1 (MQ=35)
           tGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTc                         <  1:607001/36‑1 (MQ=35)
            gCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCg                        >  1:365025/1‑36 (MQ=35)
            gCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCg                        >  2:2893441/1‑36 (MQ=35)
            gCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCg                        >  1:197647/1‑36 (MQ=35)
            gCAGATAACAAAAAACCCCGCCGGAGCGAAGTTTCg                        >  2:3749347/1‑36 (MQ=255)
             cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGt                       >  2:3478420/1‑36 (MQ=38)
              aGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTc                      <  1:1694341/36‑1 (MQ=40)
               gATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCa                     >  2:3252775/1‑36 (MQ=255)
                 tAACAAAAACCCCCGCCGGAGCGAGGTTTCGTCAGt                   <  1:2055748/36‑1 (MQ=255)
                 tAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGt                   >  1:2317463/1‑36 (MQ=255)
                   aCAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCg                 >  2:1220231/1‑36 (MQ=255)
                   aCAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCg                 >  2:1095725/1‑36 (MQ=255)
                       aaaaCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTg             <  1:2967145/36‑1 (MQ=255)
                        aaaCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGc            >  2:1129947/1‑36 (MQ=255)
                         aaCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCg           >  1:2328782/1‑36 (MQ=255)
                           ccccGCGGGAGCGAGGTTTCGTCAGTCGCCTGCGGc         <  2:2273406/36‑1 (MQ=255)
                           ccccGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGc         >  2:2059815/1‑36 (MQ=255)
                            cccGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCt        >  2:1093156/1‑36 (MQ=255)
                            cccGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCt        >  1:750801/1‑36 (MQ=255)
                            cccGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCt        >  1:1970487/1‑36 (MQ=255)
                            cccACCTGAGCAAGGGTTCGTCAGTCGCCTGCGGCt        <  2:320426/36‑1 (MQ=255)
                             ccgccgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTg       <  2:3353172/36‑1 (MQ=255)
                             ccgccgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTg       >  2:3099231/1‑36 (MQ=255)
                               gccgGAGCGAGGTTTCGTCAGTCGCCTGCGGATgag     >  1:2598037/1‑34 (MQ=37)
                                ccgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTggga    >  1:2480777/1‑34 (MQ=255)
                                ccgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGtt    >  1:1951811/1‑35 (MQ=255)
                                 cgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTgggga   >  2:2284311/1‑33 (MQ=255)
                                 cgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTaa   <  2:1900944/36‑1 (MQ=255)
                                 cgGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTaa   <  2:1938244/36‑1 (MQ=255)
                                  ggAGCGGGGTTTCGTCAGTCGCCTGCGGCTGGTAAc  <  2:2587617/36‑1 (MQ=255)
                                  ggAGCGAGGTTTCGTCAGTCGCCTGCGGTTGGGaaa  >  1:3174987/1‑35 (MQ=37)
                                  ggAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc  <  2:1622501/36‑1 (MQ=255)
                                  ggAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc  >  1:3279982/1‑36 (MQ=255)
                                  ggAGCGAGGTTTCGTCAGCCGCCTGCGGCTGGTAAc  >  2:905864/1‑36 (MQ=255)
                                  |                                   
TTGTACTGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAC  >  REL606/3813028‑3813097

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG