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

AAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATA  >  NC_000913/3888382‑3888477
                                                 |                                              
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                >  1:4342075‑M1/1‑50 (MQ=255)
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                >  1:267564‑M1/1‑50 (MQ=255)
 aaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                               >  1:2371346‑M1/1‑49 (MQ=255)
  aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcg                                              <  1:3576389‑M1/50‑3 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                           <  1:921232‑M1/50‑6 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                           >  1:3054276‑M1/1‑45 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcacc                                      <  1:2104386‑M1/50‑11 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                       >  1:2334135‑M1/1‑40 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttc                                         <  1:392817‑M1/47‑8 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttc                                         <  1:1915858‑M1/47‑8 (MQ=255)
            tgtAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccag                                    >  1:4164232‑M1/1‑38 (MQ=255)
             gtAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccag                                    >  1:779122‑M1/1‑37 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                  <  1:800033‑M1/50‑15 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                  <  1:2446617‑M1/50‑15 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                  <  1:2469637‑M1/50‑15 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                  <  1:2848264‑M1/50‑15 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                >  1:279120‑M1/1‑34 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                >  1:310528‑M1/1‑34 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                <  1:626861‑M1/50‑17 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                >  1:1649616‑M1/1‑34 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                >  1:2417200‑M1/1‑34 (MQ=255)
                 taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                               >  1:4259662‑M1/1‑33 (MQ=255)
                   atGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctc                             >  1:1069148‑M1/1‑31 (MQ=255)
                    tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                            >  1:2266209‑M1/1‑30 (MQ=255)
                    tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                            <  1:3975253‑M1/50‑21 (MQ=255)
                      tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctgg                          <  1:1673542‑M1/50‑23 (MQ=255)
                      tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctg                           <  1:1700090‑M1/49‑22 (MQ=255)
                       aGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggt                         >  1:667934‑M1/1‑27 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggt                         <  1:2916004‑M1/49‑24 (MQ=255)
                         ccTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtta                       <  1:1199868‑M1/50‑26 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                      >  1:1480187‑M1/1‑24 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                      >  1:1069231‑M1/1‑24 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                      <  1:813769‑M1/50‑27 (MQ=255)
                             gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtc                   <  1:1201050‑M1/50‑30 (MQ=255)
                              tgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtca                  <  1:3868695‑M1/50‑31 (MQ=255)
                                tCTTGCGAGGATAAGTGCcggcttcaccagtattctctggttatgtcagg                <  1:1118466‑M1/50‑33 (MQ=255)
                                tCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcagg                <  1:842131‑M1/50‑33 (MQ=255)
                                tCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcag                 <  1:700635‑M1/49‑32 (MQ=255)
                                   tGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttt             >  1:1683819‑M1/1‑15 (MQ=255)
                                     cGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttg           <  1:2844249‑M1/50‑38 (MQ=255)
                                     cGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttg           >  1:1371644‑M1/1‑13 (MQ=255)
                                     cGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtttt            <  1:1629122‑M1/49‑37 (MQ=255)
                                      gAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgc          >  1:3770302‑M1/1‑12 (MQ=255)
                                            aaGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg     >  1:2128458‑M1/1‑6 (MQ=255)
                                            aaGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg     >  1:979301‑M1/1‑6 (MQ=255)
                                             aGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcga    <  1:3696471‑M1/49‑45 (MQ=255)
                                             aGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcga    <  1:2427849‑M1/49‑45 (MQ=255)
                                              gTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaat  <  1:4012815‑M1/50‑47 (MQ=255)
                                                 |                                              
AAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATA  >  NC_000913/3888382‑3888477

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