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 34 [0] [0] 33 [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)32 (0.890)
+GA
28/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

AGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110632‑3110728
                                               |                                                 
aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAga                                                 <  1:1795356‑M1/50‑3 (MQ=255)
aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  >  1:659632‑M1/1‑48 (MQ=255)
 ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                >  1:1141070‑M1/1‑47 (MQ=255)
 ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                <  1:52752‑M1/50‑4 (MQ=255)
     tctcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:2280728‑M1/1‑43 (MQ=255)
       tcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccat                                          >  1:290061‑M1/1‑41 (MQ=255)
         agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattt                                        >  1:2124164‑M1/1‑39 (MQ=255)
         agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccat                                          <  1:975614‑M1/48‑10 (MQ=255)
          gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       <  1:216896‑M1/50‑13 (MQ=255)
          gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       <  1:2137586‑M1/50‑13 (MQ=255)
            agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttga                                     <  1:1020641‑M1/50‑15 (MQ=255)
             gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                    <  1:2233771‑M1/50‑16 (MQ=255)
              caCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaat                                   >  1:986911‑M1/1‑34 (MQ=255)
                  gACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                               >  1:1786332‑M1/1‑30 (MQ=255)
                   aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              <  1:1212059‑M1/50‑22 (MQ=255)
                   aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                               <  1:2399943‑M1/49‑21 (MQ=255)
                   aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                               >  1:511842‑M1/1‑29 (MQ=255)
                    cATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacga                             <  1:1145217‑M1/50‑23 (MQ=255)
                      tCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaac                           <  1:601976‑M1/50‑25 (MQ=255)
                          tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       >  1:1828422‑M1/1‑22 (MQ=255)
                          tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       >  1:199679‑M1/1‑22 (MQ=255)
                           ttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                      >  1:2496603‑M1/1‑21 (MQ=255)
                                 cACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:1331414‑M1/1‑15 (MQ=255)
                                  aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatc               <  1:2412580‑M1/50‑37 (MQ=255)
                                  aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:2613087‑M1/1‑14 (MQ=255)
                                   cccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:501572‑M1/1‑13 (MQ=255)
                                     cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            >  1:474687‑M1/1‑11 (MQ=255)
                                       cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          >  1:1370347‑M1/1‑9 (MQ=255)
                                       cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          >  1:1851915‑M1/1‑9 (MQ=255)
                                           gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:2765407‑M1/1‑5 (MQ=255)
                                              tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   >  1:2697861‑M1/1‑2 (MQ=255)
                                              tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   <  1:983613‑M1/50‑49 (MQ=255)
                                               agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:926978‑M1/1‑1 (MQ=255)
                                               agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:969013‑M1/1‑1 (MQ=255)
                                               |                                                 
AGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110632‑3110728

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