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

GGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAG  >  NC_000913/3906550‑3906648
                                                 |                                                 
aaaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                                   >  1:2574940‑M2/50‑50 (MQ=255)
aaaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                                   <  1:1374481‑M2/1‑1 (MQ=255)
aaaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                                   >  1:1473786‑M2/50‑50 (MQ=255)
aaaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                                   <  1:1099828‑M2/1‑1 (MQ=255)
 aaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCg                                                  <  1:2461484‑M2/2‑1 (MQ=255)
  aactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGa                                                 >  1:1358031‑M2/48‑50 (MQ=255)
  aactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGa                                                 >  1:1932457‑M2/48‑50 (MQ=255)
    ctctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACt                                               <  1:1803505‑M2/5‑1 (MQ=255)
       tgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCa                                            >  1:2447567‑M2/43‑50 (MQ=255)
        gagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAc                                           >  1:1438516‑M2/42‑50 (MQ=255)
           tgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAg                                        <  1:1391269‑M2/12‑1 (MQ=255)
             taataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGt                                       <  1:654329‑M2/13‑1 (MQ=255)
                taatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtc                                   >  1:2312500‑M2/34‑50 (MQ=255)
                 aatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtct                                  >  1:1686180‑M2/33‑50 (MQ=255)
                   tgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTctct                                >  1:1747126‑M2/31‑50 (MQ=255)
                     tagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTgg                              >  1:742940‑M2/29‑50 (MQ=255)
                        cctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTa                           <  1:135657‑M2/25‑1 (MQ=255)
                        cctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTa                           <  1:830464‑M2/25‑1 (MQ=255)
                            gtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTc                       <  1:487298‑M2/29‑1 (MQ=255)
                                    cgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGtttt                <  1:16100‑M2/36‑1 (MQ=255)
                                       ggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCt            <  1:2019677‑M2/40‑1 (MQ=255)
                                          taagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg         <  1:1849566‑M2/43‑1 (MQ=255)
                                            agtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGaa       >  1:1191652‑M2/6‑50 (MQ=255)
                                                 cGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAg  <  1:720501/50‑1 (MQ=255)
                                                 |                                                 
GGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAG  >  NC_000913/3906550‑3906648

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: