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 58 [0] [0] 53 [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)52 (0.560)
+GA
40/94 0.3 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

CAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCAC  >  NC_000913/3110631‑3110727
                                                |                                                
cAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                 <  1:7012821‑M1/50‑2 (MQ=255)
 aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAga                                                <  1:7389745‑M1/50‑3 (MQ=255)
 aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAga                                                <  1:5272723‑M1/50‑3 (MQ=255)
 aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                 >  1:1228650‑M1/1‑48 (MQ=255)
 aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                 >  1:3851721‑M1/1‑48 (MQ=255)
 aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                 >  1:6228768‑M1/1‑48 (MQ=255)
  ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                               >  1:7170707‑M1/1‑47 (MQ=255)
   gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                              <  1:6567206‑M1/50‑5 (MQ=255)
   gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                              <  1:2274679‑M1/50‑5 (MQ=255)
   gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                              <  1:4169679‑M1/50‑5 (MQ=255)
   gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                              <  1:6350781‑M1/50‑5 (MQ=255)
      tctcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                            >  1:5052981‑M1/1‑43 (MQ=255)
         cagcagCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatt                                        >  1:177479‑M1/1‑40 (MQ=255)
          agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattt                                       >  1:5511504‑M1/1‑39 (MQ=255)
          agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccat                                         <  1:4847144‑M1/48‑10 (MQ=255)
           gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                      >  1:6173640‑M1/1‑38 (MQ=255)
           gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                      <  1:4580590‑M1/50‑13 (MQ=255)
            cagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                      >  1:2109856‑M1/1‑37 (MQ=255)
             agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                     >  1:4180869‑M1/1‑36 (MQ=255)
              gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                   <  1:6362382‑M1/50‑16 (MQ=255)
               caCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaat                                  >  1:5820561‑M1/1‑34 (MQ=255)
                aCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                 >  1:4227211‑M1/1‑33 (MQ=255)
                aCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                 >  1:4928561‑M1/1‑33 (MQ=255)
                 cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                 >  1:798434‑M1/1‑32 (MQ=255)
                 cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                 >  1:2200317‑M1/1‑32 (MQ=255)
                  tGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattta                               >  1:6785073‑M1/1‑31 (MQ=255)
                  tGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattta                               >  1:4540480‑M1/1‑31 (MQ=255)
                   gACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                              >  1:1981548‑M1/1‑30 (MQ=255)
                   gACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                              >  1:1664367‑M1/1‑30 (MQ=255)
                    aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                             <  1:3101157‑M1/50‑22 (MQ=255)
                    aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                             <  1:1525403‑M1/50‑22 (MQ=255)
                       tCAATTTGTGTCTCCCGCAGCGCATAgaacgccattttgaatttacgaac                          >  1:6444694‑M1/1‑26 (MQ=255)
                       tCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaac                          <  1:2968449‑M1/50‑25 (MQ=255)
                           tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                      <  1:2707646‑M1/50‑29 (MQ=255)
                            ttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                     >  1:7410459‑M1/1‑21 (MQ=255)
                               tgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcat                  <  1:7150999‑M1/50‑33 (MQ=255)
                               tgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcat                  <  1:2145221‑M1/50‑33 (MQ=255)
                                 tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatta                <  1:6869232‑M1/50‑35 (MQ=255)
                                 tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatta                >  1:6525512‑M1/1‑16 (MQ=255)
                                 tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatta                >  1:7308394‑M1/1‑16 (MQ=255)
                                 tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatta                <  1:5134136‑M1/50‑35 (MQ=255)
                                 tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatt                 <  1:6259342‑M1/49‑34 (MQ=255)
                                   aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat               <  1:5899367‑M1/49‑36 (MQ=255)
                                     ccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcag            >  1:205763‑M1/1‑12 (MQ=255)
                                      cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt           <  1:1256468‑M1/50‑40 (MQ=255)
                                      cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt           >  1:3240581‑M1/1‑11 (MQ=255)
                                        cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga         <  1:6239733‑M1/50‑42 (MQ=255)
                                         aGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtgac        >  1:3228434‑M1/1‑8 (MQ=255)
                                          gcgcATAgaacgccattttgaatttacgaacggacgcattatcagtgact       <  1:6840485‑M1/50‑44 (MQ=255)
                                          gcgcATAgaacgccattttgaatttacgaacggacgcattatcagtgact       >  1:3321882‑M1/1‑7 (MQ=255)
                                           cgcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt      <  1:6313132‑M1/50‑45 (MQ=255)
                                            gcATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt     <  1:46166‑M1/50‑46 (MQ=255)
                                            gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt      >  1:5247135‑M1/1‑5 (MQ=255)
                                            gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt      >  1:7223396‑M1/1‑5 (MQ=255)
                                             cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttta    <  1:7437890‑M1/50‑47 (MQ=255)
                                             cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt     >  1:2487175‑M1/1‑4 (MQ=255)
                                              aTAgaacgccattttgaatttacgaacggacgcattatcagtgactttaa   >  1:4245033‑M1/1‑3 (MQ=255)
                                               tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac  <  1:6816933‑M1/50‑49 (MQ=255)
                                                |                                                
CAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCAC  >  NC_000913/3110631‑3110727

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