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 49 [0] [0] 52 [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)47 (0.850)
+GA
39/94 0.0 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

CCAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110630‑3110728
                                                 |                                                 
ccAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCata                                                   <  1:266862/50‑1 (MQ=255)
  aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  >  1:228746‑M1/1‑48 (MQ=255)
   ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                >  1:644628‑M1/1‑47 (MQ=255)
    gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                               <  1:1996824‑M1/50‑5 (MQ=255)
      ttCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:334426‑M1/1‑44 (MQ=255)
         tcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccat                                          >  1:690755‑M1/1‑41 (MQ=255)
           agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatt                                         <  1:3614348‑M1/49‑11 (MQ=255)
           agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccat                                          <  1:3578971‑M1/48‑10 (MQ=255)
           agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgcca                                           <  1:1168051‑M1/47‑9 (MQ=255)
            gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       <  1:788630‑M1/50‑13 (MQ=255)
             cagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                      <  1:3113161‑M1/50‑14 (MQ=255)
             cagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                      <  1:2698015‑M1/50‑14 (MQ=255)
              agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttga                                     <  1:1117214‑M1/50‑15 (MQ=255)
                caCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaat                                   >  1:3266070‑M1/1‑34 (MQ=255)
                caCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaat                                   <  1:1369427‑M1/50‑17 (MQ=255)
                 aCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                  <  1:2393988‑M1/50‑18 (MQ=255)
                  cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattt                                 <  1:4073728‑M1/50‑19 (MQ=255)
                   tGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattta                                >  1:495455‑M1/1‑31 (MQ=255)
                    gACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                               >  1:887372‑M1/1‑30 (MQ=255)
                     aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              <  1:2121550‑M1/50‑22 (MQ=255)
                      cATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              >  1:3114500‑M1/1‑28 (MQ=255)
                          aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgg                         <  1:1735542‑M1/50‑27 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:541517‑M1/50‑29 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:295902‑M1/50‑29 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:2481287‑M1/50‑29 (MQ=255)
                              tgtgTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                      >  1:2230135‑M1/1‑20 (MQ=255)
                              tgtgTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                      >  1:4077135‑M1/1‑20 (MQ=255)
                               gtgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgca                    >  1:3770932‑M1/1‑19 (MQ=255)
                               gtgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgca                    >  1:2305121‑M1/1‑19 (MQ=255)
                               gtgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgc                     <  1:3800215‑M1/49‑31 (MQ=255)
                                tgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcat                   >  1:589152‑M1/1‑18 (MQ=255)
                                tgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcat                   <  1:693985‑M1/50‑33 (MQ=255)
                                  tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatta                 <  1:1566721‑M1/50‑35 (MQ=255)
                                    aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatc               <  1:4313433‑M1/50‑37 (MQ=255)
                                    aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                <  1:1094931‑M1/49‑36 (MQ=255)
                                     cccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatca              <  1:4145141‑M1/50‑38 (MQ=255)
                                     cccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatca              <  1:904484‑M1/50‑38 (MQ=255)
                                      ccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcag             <  1:1453240‑M1/50‑39 (MQ=255)
                                       cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            <  1:3129610‑M1/50‑40 (MQ=255)
                                       cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            >  1:2881401‑M1/1‑11 (MQ=255)
                                            cgcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       <  1:939376‑M1/50‑45 (MQ=255)
                                             gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:566402‑M1/1‑5 (MQ=255)
                                             gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:1609205‑M1/1‑5 (MQ=255)
                                              cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttta     <  1:4394357‑M1/50‑47 (MQ=255)
                                              cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttta     <  1:2301576‑M1/50‑47 (MQ=255)
                                              cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      >  1:647885‑M1/1‑4 (MQ=255)
                                               aTAgaacgccattttgaatttacgaacggacgcattatcagtgactttaa    >  1:1291053‑M1/1‑3 (MQ=255)
                                                tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   <  1:2998215‑M1/50‑49 (MQ=255)
                                                 agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:1111635‑M1/1‑1 (MQ=255)
                                                 |                                                 
CCAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110630‑3110728

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