Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,110,680 25702 bp→GA [yghD]pitB 20 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3110680 3136381 25702 26 [0] [0] 27 [yghD]–pitB [yghD],yghE,yghF,yghG,pppA,yghJ,glcA,glcB,glcG,glcF,glcE,glcD,glcC,yghO,insH‑9,yghQ,yghR,yghS,yghT,pitB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 31106790 (0.000)26 (0.890)
+GA
21/94 0.1 100% coding (448/537 nt) yghD putative type II secretion system M‑type protein
?NC_000913 3136382 = 0 (0.000)intergenic (‑11/+281) pitB/gss metal phosphate:H(+) symporter PitB/fused glutathionylspermidine amidase/glutathionylspermidine synthetase

GCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110634‑3110728
                                             |                                                 
gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                               <  1:440080‑M1/50‑5 (MQ=255)
gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                               <  1:2484749‑M1/50‑5 (MQ=255)
  ttCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:2220407‑M1/1‑44 (MQ=255)
    ctcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgcca                                           >  1:1080972‑M1/1‑42 (MQ=255)
      cagcagCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatt                                         >  1:772916‑M1/1‑40 (MQ=255)
      cagcagCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatt                                         >  1:1481762‑M1/1‑40 (MQ=255)
        gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       <  1:156646‑M1/50‑13 (MQ=255)
         cagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       >  1:1151905‑M1/1‑37 (MQ=255)
                   aTCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaa                            >  1:950152‑M1/1‑27 (MQ=255)
                     cAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacg                          >  1:853372‑M1/1‑25 (MQ=255)
                       aTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgga                        >  1:1092894‑M1/1‑23 (MQ=255)
                               cACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:279934‑M1/1‑15 (MQ=255)
                                  ccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcag             <  1:1590515‑M1/50‑39 (MQ=255)
                                   cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            <  1:1432582‑M1/50‑40 (MQ=255)
                                    gCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtg           <  1:1400651‑M1/50‑41 (MQ=255)
                                     cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          <  1:1631806‑M1/50‑42 (MQ=255)
                                     cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          >  1:2044654‑M1/1‑9 (MQ=255)
                                     cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          <  1:669975‑M1/50‑42 (MQ=255)
                                      aGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtgac         <  1:2661872‑M1/50‑43 (MQ=255)
                                      aGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtgac         <  1:1329886‑M1/50‑43 (MQ=255)
                                         gcATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      <  1:738033‑M1/50‑46 (MQ=255)
                                         gcATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      <  1:1250581‑M1/50‑46 (MQ=255)
                                          cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttta     <  1:1577697‑M1/50‑47 (MQ=255)
                                          cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      >  1:1879001‑M1/1‑4 (MQ=255)
                                           aTAgaacgccattttgaatttacgaacggacgcattatcagtgactttaa    >  1:1123080‑M1/1‑3 (MQ=255)
                                             agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:1898602‑M1/1‑1 (MQ=255)
                                             |                                                 
GCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110634‑3110728

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