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 52 [0] [0] 49 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)47 (0.940) 35/96 0.0 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

AAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTG  >  NC_000913/3888382‑3888480
                                                 |                                                 
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                   >  1:2573947‑M1/1‑50 (MQ=255)
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                   >  1:2826387‑M1/1‑50 (MQ=255)
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                   <  1:1779318‑M1/50‑1 (MQ=255)
  aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcg                                                 <  1:514782‑M1/50‑3 (MQ=255)
  aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                                  >  1:185148‑M1/1‑48 (MQ=255)
   aaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcga                                                >  1:1918948‑M1/1‑47 (MQ=255)
    aCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgac                                               <  1:2770529‑M1/50‑5 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                              <  1:1103059‑M1/50‑6 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                              <  1:1217502‑M1/50‑6 (MQ=255)
      tctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                             >  1:1874753‑M1/1‑44 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttca                                           <  1:1618665‑M1/48‑9 (MQ=255)
           gtgtAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcacca                                        <  1:2613483‑M1/50‑12 (MQ=255)
            tgtAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccag                                       <  1:1028021‑M1/50‑13 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                     >  1:2456365‑M1/1‑36 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                     >  1:742502‑M1/1‑36 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagt                                      <  1:3752785‑M1/49‑14 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagt                                      <  1:1764696‑M1/49‑14 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagt                                      <  1:1577598‑M1/49‑14 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                   >  1:2412345‑M1/1‑34 (MQ=255)
                   atGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctc                                <  1:3053069‑M1/50‑20 (MQ=255)
                   atGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattct                                 <  1:3025015‑M1/49‑19 (MQ=255)
                    tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                               >  1:1378845‑M1/1‑30 (MQ=255)
                    tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                               <  1:214877‑M1/50‑21 (MQ=255)
                    tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                               <  1:1207298‑M1/50‑21 (MQ=255)
                     gTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctg                              >  1:1800708‑M1/1‑29 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtt                           <  1:3543091‑M1/50‑25 (MQ=255)
                         ccTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtta                          >  1:1326467‑M1/1‑25 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                         <  1:367376‑M1/50‑27 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                         >  1:3932013‑M1/1‑24 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                         >  1:1576218‑M1/1‑24 (MQ=255)
                            cGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       >  1:631263‑M1/1‑22 (MQ=255)
                              tgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtca                     <  1:1529773‑M1/50‑31 (MQ=255)
                              tgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtca                     <  1:1546079‑M1/50‑31 (MQ=255)
                                tCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcagg                   >  1:2504397‑M1/1‑18 (MQ=255)
                                    gCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtttt               >  1:2071862‑M1/1‑14 (MQ=255)
                                    gCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttt                <  1:2026330‑M1/49‑36 (MQ=255)
                                     cGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtttt               <  1:2511132‑M1/49‑37 (MQ=255)
                                      gAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgc             >  1:2054382‑M1/1‑12 (MQ=255)
                                      gAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgc             >  1:1982797‑M1/1‑12 (MQ=255)
                                         gATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctg          <  1:2276093‑M1/50‑42 (MQ=255)
                                          aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc         <  1:3201806‑M1/50‑43 (MQ=255)
                                          aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc         >  1:1288896‑M1/1‑8 (MQ=255)
                                          aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc         <  1:10171‑M1/50‑43 (MQ=255)
                                           tAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc         >  1:2406274‑M1/1‑7 (MQ=255)
                                            aaGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg        >  1:839942‑M1/1‑6 (MQ=255)
                                            aaGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg        >  1:2364745‑M1/1‑6 (MQ=255)
                                             aGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaa      >  1:498410‑M1/1‑5 (MQ=255)
                                               tGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaatg    >  1:2213256‑M1/1‑3 (MQ=255)
                                               tGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaatg    >  1:187938‑M1/1‑3 (MQ=255)
                                               tGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaatg    >  1:1011173‑M1/1‑3 (MQ=255)
                                                 ccgacttcaccagtattctctggttatgtcaggttttgcctgcgaatgca  <  1:751698‑M1/50‑50 (MQ=255)
                                                 ccgacttcaccagtattctctggttatgtcaggttttgcctgcgaatg    >  1:2689358‑M1/1‑1 (MQ=255)
                                                 |                                                 
AAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTG  >  NC_000913/3888382‑3888480

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