New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 185468326 (0.560)5 (0.110) 3/54 NT 14.5% coding (79/255 nt) yoaF conserved outer membrane protein
?REL606 = 4266654 33 (0.550)intergenic (‑350/‑43) acs/nrfA acetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome

CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGA  >  REL606/4266619‑4266689
                                   |                                   
cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                     >  1:7925611/1‑36 (MQ=255)
cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                     >  1:2285525/1‑36 (MQ=255)
cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                     <  1:4828138/36‑1 (MQ=255)
cAGTGTAAGTGCGGTTACTGTCGTGAAAAACAATaa                                     <  1:5407866/36‑1 (MQ=255)
 agggtAGGTGCGGTTACTGTCGTGTAAAAAAATaaa                                    >  1:3269903/4‑36 (MQ=25)
 aGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAa                                    >  1:1397330‑M2/1‑35 (MQ=255)
  gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa                                   <  1:7278335‑M2/36‑3 (MQ=255)
  gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa                                   >  1:2262089‑M2/1‑34 (MQ=255)
   tgtAGGTGCGGTTACTGTCGTGAAAATCAATataaa                                  >  1:3913464/1‑36 (MQ=37)
   tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                                  >  1:3354270‑M2/1‑33 (MQ=255)
   tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                                  <  1:2036716‑M2/36‑4 (MQ=255)
   tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                                  <  1:2729121‑M2/36‑4 (MQ=255)
    gtAGGTGCGGTTACTGTTGTGAAAAACAATAAaaac                                 >  1:6458527‑M2/1‑32 (MQ=255)
    gtAGGCGCGCTTACTGTAGTGAAAAACAATAAAAAc                                 <  1:4256303/36‑1 (MQ=25)
     tAGGTGCGGTTACTGTCGTGAAAAACAATAAaaacc                                <  1:1786231‑M2/36‑6 (MQ=255)
     tAGGTGCGGTTACTGTCGTGAAAAACAATAAaaacc                                >  1:5475810‑M2/1‑31 (MQ=255)
        gggCGGTTACTGCAGTGAAAAACAATAAAAACCGcc                             <  1:6051002/34‑1 (MQ=37)
         tGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCa                            >  1:4912596/1‑36 (MQ=255)
           cGGTTACTGTCGTGAAAAACAATAAAAACCGCCAtt                          >  1:868948/1‑36 (MQ=255)
           cGGTTACTGTCGTGAAAAAAAATAAAAACCGCCAtt                          >  1:61833/1‑36 (MQ=255)
            ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                         >  1:1215817/1‑36 (MQ=255)
            ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                         >  1:6793296/1‑36 (MQ=255)
             gTTACTGTCGTAAAAAACAATAAAAACCGCCATTGc                        >  1:3563058/1‑36 (MQ=255)
              ttACTGTCGTGAAAAACAATAAAAACCGCCATTGca                       >  1:1563161/1‑36 (MQ=255)
                    tCGTGAAAAACAAGAAAAACCGCCATTGCAACAATg                 >  1:2285148/1‑36 (MQ=255)
                      gTGAAAAACAATAAAAACCGCCATTGCAACAACGgc               >  1:4217963/1‑36 (MQ=255)
                      aaaaaaaaAAAAAAAACCCGCCATTGCAACAATGgc               <  1:5309197/33‑1 (MQ=16)
                       tGAAAAACAATAAAAACCGCCATTGCAACAATGgcg              >  1:7850141/1‑36 (MQ=255)
                       aaaaaaaAAAAAAAACCCGCCATTGCAACAATGgcg              <  1:122419/34‑1 (MQ=21)
                        gAAAAACAATAAAAACCGCCATTGCAACAATGgcgc             >  1:7584869/1‑36 (MQ=255)
                        aaaaaaAAAAAAAAACCGCCATTGCAACAATGgcgc             <  1:3108806/35‑1 (MQ=25)
                         aaaaaCAATAAAAACCGCCATTGCAACAATGGCGCa            >  1:2059766/1‑36 (MQ=255)
                          aaaaCAATAAAAACCGCCATTGCAACAATGGCGCaa           >  1:1845099/1‑36 (MQ=255)
                          aaaaCAAAAAaaaccgccgcagattggtatggcaaa           >  1:4908588‑M2/1‑10 (MQ=255)
                           aaaCAATAAAAACCGCCCTTGCAACAATGCCGCaac          >  1:6462075/1‑35 (MQ=25)
                           aaaAAAAAAAAACCGCCATTGCAACAATGGCGCAAt          <  1:3095111/36‑1 (MQ=37)
                            aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt         >  1:5047552/1‑36 (MQ=255)
                            aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt         >  1:1741131/1‑36 (MQ=255)
                            aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt         >  1:1596554/1‑36 (MQ=255)
                             aCAATAAAAACCGCCATTGCAACAATGGCGCAATTc        >  1:819376/1‑36 (MQ=255)
                               aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg      >  1:6620795/1‑36 (MQ=255)
                               aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg      >  1:1433812/1‑36 (MQ=255)
                               aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg      >  1:6525060/1‑36 (MQ=255)
                                aaaaaaaCCGCCATTGCATCAATAGCGCAATTCGGa     <  1:5497865/34‑1 (MQ=25)
                                aTAAAAACCGCCATTGCAACAATGGCGCAATTTGGa     >  1:4116640/1‑36 (MQ=255)
                                 aaaaaaCCGCCATTGCAACAATGGCGCAATTCGGAt    <  1:6208594/35‑1 (MQ=255)
                                  caaacccgccgcagattggtatggcaaatccggcgt   >  1:4617256‑M2/2‑2 (MQ=255)
                                  caaaaccgccgcagattggtatggcaaatccggcgt   >  1:5053036‑M2/2‑2 (MQ=255)
                                  aaaacCCGCCATTGCAACAATGGCGCAATTCGGATg   <  1:3134978/36‑1 (MQ=255)
                                  aaaaaccgccgcagattggtatggcaaatccggcgt   >  1:6202243‑M2/1‑2 (MQ=255)
                                   aaaaccgccgcagattggtatggcaaatccggcgtc  >  1:1384665‑M2/1‑1 (MQ=255)
                                   |                                   
CAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGA  >  REL606/4266619‑4266689

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG