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 48 [0] [0] 46 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)46 (0.810) 35/96 0.1 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

GGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAAT  >  NC_000913/3906550‑3906644
                                                 |                                             
aaaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcc                                               >  1:2371346‑M2/50‑50 (MQ=255)
 aaactctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCg                                              <  1:3576389‑M2/2‑1 (MQ=255)
    ctctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACt                                           >  1:3054276‑M2/46‑50 (MQ=255)
    ctctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACt                                           <  1:921232‑M2/5‑1 (MQ=255)
         agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTc                                         <  1:1915858‑M2/7‑1 (MQ=255)
         agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTc                                         <  1:392817‑M2/7‑1 (MQ=255)
         agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAcc                                      <  1:2104386‑M2/10‑1 (MQ=255)
         agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAc                                       >  1:2334135‑M2/41‑49 (MQ=255)
           tgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAg                                    >  1:4164232‑M2/39‑50 (MQ=255)
            gtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAg                                    >  1:779122‑M2/38‑49 (MQ=255)
             taataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTa                                  <  1:800033‑M2/14‑1 (MQ=255)
             taataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTa                                  <  1:2446617‑M2/14‑1 (MQ=255)
             taataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTa                                  <  1:2469637‑M2/14‑1 (MQ=255)
             taataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTa                                  <  1:2848264‑M2/14‑1 (MQ=255)
               ataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAtt                                >  1:1649616‑M2/35‑50 (MQ=255)
               ataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAtt                                >  1:310528‑M2/35‑50 (MQ=255)
               ataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAtt                                <  1:626861‑M2/16‑1 (MQ=255)
               ataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAtt                                >  1:2417200‑M2/35‑50 (MQ=255)
               ataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAtt                                >  1:279120‑M2/35‑50 (MQ=255)
                taatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtc                               >  1:4259662‑M2/34‑50 (MQ=255)
                  atgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtctc                             >  1:1069148‑M2/32‑50 (MQ=255)
                   tgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTctct                            <  1:3975253‑M2/20‑1 (MQ=255)
                   tgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTctct                            >  1:2266209‑M2/31‑50 (MQ=255)
                     tagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTgg                          <  1:1673542‑M2/22‑1 (MQ=255)
                     tagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTg                           <  1:1700090‑M2/21‑1 (MQ=255)
                      agcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGt                         >  1:667934‑M2/28‑50 (MQ=255)
                       gcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGt                         <  1:2916004‑M2/23‑1 (MQ=255)
                        cctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTa                       <  1:1199868‑M2/25‑1 (MQ=255)
                         ctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTAt                      <  1:813769‑M2/26‑1 (MQ=255)
                         ctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTAt                      >  1:1480187‑M2/25‑50 (MQ=255)
                         ctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTAt                      >  1:1069231‑M2/25‑50 (MQ=255)
                            gtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTc                   <  1:1201050‑M2/29‑1 (MQ=255)
                             tgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCa                  <  1:3868695‑M2/30‑1 (MQ=255)
                               tcttgcgaggataagtgcCGGCTTCACCAGTATTCTCTGGTTATGTCAgg                <  1:1118466‑M2/32‑1 (MQ=255)
                               tcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAgg                <  1:842131‑M2/32‑1 (MQ=255)
                               tcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAg                 <  1:700635‑M2/31‑1 (MQ=255)
                                  tgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGttt             >  1:1683819‑M2/16‑50 (MQ=255)
                                    cgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGtttt            <  1:1629122‑M2/36‑1 (MQ=255)
                                    cgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTg           <  1:2844249‑M2/37‑1 (MQ=255)
                                    cgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTg           >  1:1371644‑M2/14‑50 (MQ=255)
                                     gaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGc          >  1:3770302‑M2/13‑50 (MQ=255)
                                           aagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg     >  1:2128458‑M2/7‑49 (MQ=255)
                                           aagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg     >  1:979301‑M2/7‑49 (MQ=255)
                                            agtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGa    <  1:3696471‑M2/44‑1 (MQ=255)
                                            agtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGa    <  1:2427849‑M2/44‑1 (MQ=255)
                                             gtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAAt  <  1:4012815‑M2/46‑1 (MQ=255)
                                                 |                                             
GGTGATTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAAT  >  NC_000913/3906550‑3906644

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