New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 374957240 (0.880)8 (0.180) 4/58 5.2 16.5% intergenic (‑65/+33) mutM/rpmG formamidopyrimidine‑DNA glycosylase/50S ribosomal protein L33
?REL606 = 4157631 41 (0.740)coding (278/546 nt) nusG transcription antitermination protein NusG
Rejected: Position hash score below cutoff.

ATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTA  >  REL606/3749540‑3749607
                                |                                   
atGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAg                                  >  1:1243954/1‑36 (MQ=255)
atGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAg                                  <  2:1024283/36‑1 (MQ=255)
  gggggCAAGCAGATAACAAAAAACCCCGCCGGAGCg                                >  1:1008651/1‑36 (MQ=255)
  gggggCAAGCAGATAACAAAAAACCCCGCCGGAGCg                                <  1:547847/36‑1 (MQ=255)
  gggggCAAGCAGATAACAAAAAACCCCGCCGGAGCg                                <  1:1950985/36‑1 (MQ=255)
    gggAAAGCAGATAACAAAAAACCCCGCCGGAGCGAg                              <  1:216026/36‑1 (MQ=255)
     ggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAgg                             <  2:1750599/36‑1 (MQ=255)
       cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtt                           >  1:3640510/1‑36 (MQ=38)
       cAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtt                           >  2:137049/1‑36 (MQ=38)
       aaaGCGGATAACAAAAAGCCCCGCCTGAGCGAGGtt                           <  2:732741/35‑1 (MQ=11)
        aaGCAGATAACAAAAAACCCCGCCGGAGCGAGGttt                          >  1:2034896/1‑36 (MQ=35)
        aaGCAGATAACAAAAAACCCCGCCGGAGCGAGGttt                          >  1:1500986/1‑36 (MQ=35)
        aaGCAGATAACAAAAAACCCCGCCGGAGCGAGGttt                          >  1:1264984/1‑36 (MQ=35)
           cAGATAACAAAAAACCCCGCCGGAGCGAGGtttttt                       <  2:151412/36‑1 (MQ=38)
            aGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                      <  1:934185/36‑1 (MQ=40)
            aGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTg                      >  2:2702632/1‑36 (MQ=40)
             gATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGt                     >  1:43552/1‑36 (MQ=255)
                aaCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTAc                  >  1:2392621/1‑36 (MQ=255)
                 aCAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACa                 >  2:626704/1‑36 (MQ=255)
                  cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACAt                >  2:506770/1‑36 (MQ=255)
                  cAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACAt                <  2:1243928/36‑1 (MQ=255)
                   aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATc               >  1:399963/1‑36 (MQ=255)
                   aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATc               >  2:3291881/1‑36 (MQ=255)
                   aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATc               <  1:1518436/36‑1 (MQ=255)
                    aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCa              <  1:82791/36‑1 (MQ=255)
                    aaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCa              <  1:1343447/36‑1 (MQ=255)
                    aaaaaCCCCGCCGGAGCGACGTTTTTTCTTACATCa              >  1:2676553/1‑36 (MQ=18)
                     aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCaa             >  1:1335986/1‑36 (MQ=255)
                     aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCaa             <  2:2010168/36‑1 (MQ=255)
                      aaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCaaa            >  1:3652427/1‑36 (MQ=255)
                      aaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCaaa            >  2:219051/1‑36 (MQ=255)
                      aaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCaaa            >  2:592368/1‑36 (MQ=255)
                       aaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAg           <  2:1700059/36‑1 (MQ=255)
                        aCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc          <  2:1306227/36‑1 (MQ=255)
                        aCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc          <  1:53793/36‑1 (MQ=255)
                        aCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGc          >  2:2875780/1‑36 (MQ=255)
                         ccccGCCGGAGCGAGGTTTTTTGTTACATCAAAGCg         >  1:2298546/1‑36 (MQ=255)
                           ccgccgGAGCGAGGTTTCTTGATACATCAAAGTgag       >  1:229388/1‑36 (MQ=18)
                               cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatt   <  2:560316/36‑1 (MQ=255)
                               cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAatt   >  2:1178958/1‑36 (MQ=255)
                                ggAGCGAGGTTTTTTGTTACATCAAAGCGAGAatta  >  2:743421/1‑36 (MQ=255)
                                |                                   
ATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTA  >  REL606/3749540‑3749607

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG