Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,888,432 Δ18,167 bp tnaCbglG 17 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3888432 3906598 18167 24 [0] [0] 24 tnaC–bglG tnaC,tnaA,tnaB,mdtL,yidZ,yieE,yieF,adeP,yieH,cbrB,cbrC,yieK,yieL,bglH,bglB,bglF,bglG

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 38884310 (0.000)23 (0.770) 19/96 0.2 100% intergenic (+239/‑4) mnmE/tnaC 5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit/tnaAB operon leader peptide
?NC_000913 3906599 = 0 (0.000)intergenic (‑32/+254) bglG/phoU transcriptional antiterminator BglG/negative regulator of the pho regulon

AAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAAT  >  NC_000913/3888383‑3888476
                                                |                                             
aaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                              <  1:1099828‑M1/50‑2 (MQ=255)
aaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                              >  1:2574940‑M1/1‑49 (MQ=255)
aaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                              >  1:1473786‑M1/1‑49 (MQ=255)
aaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                              <  1:1374481‑M1/50‑2 (MQ=255)
 aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCt                                              >  1:1631413‑M1/1‑48 (MQ=255)
 aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcg                                             <  1:2461484‑M1/50‑3 (MQ=255)
  aaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcga                                            >  1:1358031‑M1/1‑47 (MQ=255)
  aaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcga                                            >  1:1932457‑M1/1‑47 (MQ=255)
    ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                          <  1:1803505‑M1/50‑6 (MQ=255)
       tGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttca                                       >  1:2447567‑M1/1‑42 (MQ=255)
        gAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                      >  1:1438516‑M1/1‑41 (MQ=255)
           tgtAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccag                                   <  1:1391269‑M1/50‑13 (MQ=255)
             taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagt                                  <  1:654329‑M1/49‑14 (MQ=255)
                taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                              >  1:2312500‑M1/1‑33 (MQ=255)
                 aatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattct                             >  1:1686180‑M1/1‑32 (MQ=255)
                   tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                           >  1:1747126‑M1/1‑30 (MQ=255)
                     tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctgg                         >  1:742940‑M1/1‑28 (MQ=255)
                        ccTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtta                      <  1:830464‑M1/50‑26 (MQ=255)
                        ccTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtta                      <  1:135657‑M1/50‑26 (MQ=255)
                            gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtc                  <  1:487298‑M1/50‑30 (MQ=255)
                                    cGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtttt           <  1:16100‑M1/49‑37 (MQ=255)
                                       ggATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcct       <  1:2019677‑M1/50‑41 (MQ=255)
                                          tAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg    <  1:1849566‑M1/50‑44 (MQ=255)
                                            aGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaa  >  1:1191652‑M1/1‑5 (MQ=255)
                                                |                                             
AAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAAT  >  NC_000913/3888383‑3888476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: