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.

CCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCG  >  REL606/4157596‑4157666
                                   |                                   
ccTGGTGCGGGGCGTACCGCGTGTGATGGGCTTCAt                                     <  1:1502251/36‑1 (MQ=37)
ccTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                                     <  1:826601/36‑1 (MQ=255)
ccTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                                     >  2:2487183/1‑36 (MQ=255)
 ttGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATc                                    >  1:1431629/2‑36 (MQ=255)
 cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATc                                    <  1:3691591/36‑1 (MQ=255)
  tGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcg                                   <  2:1850802/36‑1 (MQ=255)
  tGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcg                                   <  2:3691640/36‑1 (MQ=255)
   ggTGCGCAGCGTACCGCGTGTGATGGGCTTCATcgg                                  <  1:1439329/36‑1 (MQ=255)
   ggTCCGCAGCGTACCGCGTGTGATGGGCTTCATcgg                                  <  2:268689/36‑1 (MQ=255)
    gTGCGCAGCGTCCCGCGGGTGATGGGCTTCATcggc                                 <  2:3285704/36‑1 (MQ=37)
     tGCTCAGCATACCGCGTGTGAGGGGCTTCATcggcg                                <  1:263605/36‑1 (MQ=25)
       cgcAGCGTACCGCGTGTGATGGGCTTCATcggcggg                              >  1:241148/1‑35 (MQ=255)
        gcAGCGTACCGCGTGTGATGGGCTTCATcggcgggg                             >  2:3145682/1‑34 (MQ=255)
        gcAGCGTACCGCGTGTGATGGGCTTCATcggcgggg                             >  2:2501921/1‑34 (MQ=255)
        gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTa                             <  2:809211/36‑1 (MQ=255)
         cAGCGTACCGCGTGTGATGGGCTTCATcggcgggag                            >  1:2425187/1‑33 (MQ=255)
         cAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTAc                            <  1:3261049/36‑1 (MQ=255)
          aGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACt                           <  2:494431/36‑1 (MQ=255)
           gCGTACCGCGTGTGATGGGCTTCATCGGCGGTACtt                          <  2:3415879/36‑1 (MQ=255)
             ttACCGCGTGTGATGGGCTTCATCGGCGGTACTTcc                        <  2:2278758/35‑1 (MQ=255)
              tACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCg                       <  2:2846775/36‑1 (MQ=255)
               tccGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa                      <  2:502784/35‑1 (MQ=255)
               aCCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa                      <  1:402234/36‑1 (MQ=255)
                 cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATc                    <  1:2448522/36‑1 (MQ=255)
                 cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATc                    <  1:1487504/36‑1 (MQ=255)
                 cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATc                    <  2:2325288/36‑1 (MQ=255)
                  gggtgtGATGGGCTTCATCGGCGGTACTTCCGATCg                   <  1:492168/34‑1 (MQ=255)
                   cgTGTGATGGGCTTCATCGGCGGTACTGCCGATCGt                  <  2:853277/36‑1 (MQ=255)
                     tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTcc                <  2:2124747/36‑1 (MQ=255)
                     tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTcc                <  2:1456311/36‑1 (MQ=255)
                     tgtgATGGGCTTCATCGGCGGGGGTTCCGTTCGttc                >  1:2416163/1‑34 (MQ=16)
                      gtgATGGGCTTCATCGGCGGGGATTGCGATCGTCCt               >  1:2853273/1‑36 (MQ=19)
                       tgAGGGGCTTCATCGGCGGTACTTCCGATCGTCCTg              <  1:3732744/36‑1 (MQ=255)
                        gATGGGCTTCATCGGCGGTACTTCCGATCGTCCTgc             <  2:557402/36‑1 (MQ=255)
                        gATGGGCTGCATCGGCGGTACTTCCGATCGTCCTgc             <  2:3064921/36‑1 (MQ=255)
                         aTGGTCTTCATCGGCGGTACTTCCGATCGTCCTgcg            <  2:1191329/36‑1 (MQ=255)
                           gggTTTCATCGGCGGTACTTCCGATCGTCCTGCGcc          <  2:2318373/36‑1 (MQ=255)
                           gggCTTCATCGGCGGTACTTCCGATCGTCCTGCGcc          <  1:2399836/36‑1 (MQ=255)
                            ggCTTCATCGGCGGCACTTTCCGTCGTCCTGCCCCa         >  1:2516364/1‑36 (MQ=12)
                             gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCaa        <  2:336922/36‑1 (MQ=255)
                             gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCaa        <  1:2725992/36‑1 (MQ=255)
                              cTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAAt       <  1:423584/36‑1 (MQ=255)
                              cTTCATCGGCGGTACTCCCCATCGTCCTTCGCCact       >  1:2028682/1‑34 (MQ=21)
                               ttCATCGGCGGTACTTCCGATCGTCCTGCGCCAATc      <  2:3236058/36‑1 (MQ=255)
                                tCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCa     <  2:2204208/36‑1 (MQ=255)
                                 cATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAg    <  1:1277124/36‑1 (MQ=255)
                                  aTCGGCGGTACTTCCGATCGTCCTGCGCCAGTCAGc   <  2:2707049/36‑1 (MQ=255)
                                  aTCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGc   <  2:1115154/36‑1 (MQ=255)
                                   tCGGCGGTACTTCCGAACCTCCCGCGGCAATCAACg  >  2:442280/1‑36 (MQ=12)
                                   |                                   
CCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCG  >  REL606/4157596‑4157666

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG