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 27 [0] [0] 25 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)24 (0.820) 20/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

GTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATG  >  NC_000913/3906551‑3906645
                                                |                                              
aaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                                >  1:1285790‑M2/49‑49 (MQ=255)
aaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                                >  1:653951‑M2/49‑49 (MQ=255)
 aactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGa                                              <  1:1348536‑M2/3‑1 (MQ=255)
  actctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGAc                                             <  1:229712‑M2/4‑1 (MQ=255)
        agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCa                                         <  1:190099‑M2/8‑1 (MQ=255)
               taatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtc                                >  1:115069‑M2/34‑50 (MQ=255)
               taatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtc                                >  1:201300‑M2/34‑50 (MQ=255)
                aatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtct                               >  1:1312146‑M2/33‑50 (MQ=255)
                aatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtct                               >  1:1721237‑M2/33‑50 (MQ=255)
                aatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtct                               >  1:2255104‑M2/33‑50 (MQ=255)
                  tgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTctct                             <  1:1105667‑M2/20‑1 (MQ=255)
                    tagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTgg                           >  1:1888872‑M2/29‑50 (MQ=255)
                    tagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTg                            <  1:1381338‑M2/21‑1 (MQ=255)
                     agcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGt                          >  1:946951‑M2/28‑50 (MQ=255)
                      gcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGtt                         <  1:546668‑M2/24‑1 (MQ=255)
                           gtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTc                    <  1:951828‑M2/29‑1 (MQ=255)
                            tgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCa                   <  1:591709‑M2/30‑1 (MQ=255)
                              tcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAg                  <  1:354006‑M2/31‑1 (MQ=255)
                               cttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGt                >  1:1238849‑M2/18‑50 (MQ=255)
                                   cgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGtttt             <  1:487691‑M2/36‑1 (MQ=255)
                                    gaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGc           >  1:140873‑M2/13‑50 (MQ=255)
                                       gataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTg        >  1:1662197‑M2/10‑50 (MQ=255)
                                          aagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg      >  1:1169679‑M2/7‑49 (MQ=255)
                                             tgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATg  <  1:1814686‑M2/47‑1 (MQ=255)
                                               cCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATg  >  1:1772527‑M2/2‑48 (MQ=255)
                                                |                                              
GTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATG  >  NC_000913/3906551‑3906645

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.