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 76 [0] [0] 84 [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)68 (1.000)
+GA
34/64 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

GCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGG  >  NC_000913/3110645‑3110713
                                  |                                  
gcaCTGACATCAATTTGTGTCACCCGCAGCGCata                                    <  1:6348890/35‑1 (MQ=255)
gcaCTGACATCAATTTGTGTCACCCGCAGCGCata                                    <  1:2119275/35‑1 (MQ=255)
 cACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                   <  1:2376418‑M1/35‑2 (MQ=255)
 cACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                   >  1:5836306‑M1/1‑34 (MQ=255)
 cACTGACATCAATTTGTGTCACCCACAGCGCATAg                                   <  1:4168272‑M1/35‑2 (MQ=255)
  aCTGACATCAATTTGTGTCACCCGCAGCGCATAga                                  <  1:3951322‑M1/35‑3 (MQ=255)
  aCTGACATCAATTTGTGTCACCCGCAGCGCATAga                                  >  1:96087‑M1/1‑33 (MQ=255)
  aCTGACATCAATTTGTGTCACCCGCAGCGCATAga                                  >  1:6845459‑M1/1‑33 (MQ=255)
   cTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                 >  1:4448670‑M1/1‑32 (MQ=255)
   cTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                 <  1:2154487‑M1/35‑4 (MQ=255)
   cTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                 >  1:4613192‑M1/1‑32 (MQ=255)
   cTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                 >  1:6497918‑M1/1‑32 (MQ=255)
    tGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                >  1:981445‑M1/1‑31 (MQ=255)
    tGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                >  1:1194208‑M1/1‑31 (MQ=255)
      aCATCAATTTGTGTCACCCGCAGCGCATAgaacgc                              <  1:7371379‑M1/35‑7 (MQ=255)
      aCATCAATTTGTGTCACCCGCAGCGCATAgaacgc                              <  1:5830421‑M1/35‑7 (MQ=255)
       cATCAATTTGTGTCACCCGCAGCGCATAgaacgcc                             <  1:416532‑M1/35‑8 (MQ=255)
        aTCAATTTGTGTCACCCGCAGCGCATAgaacgcca                            <  1:4840402‑M1/35‑9 (MQ=255)
         tCAATTTGTGTCACCCGCAGCGCATAgaacgccat                           <  1:5598225‑M1/35‑10 (MQ=255)
          cAATTTGTGTCACCCGCAGCGCATAgaacgccatt                          <  1:4680584‑M1/35‑11 (MQ=255)
          cAATTTGTGTCACCCGCAGCGCATAgaacgccatt                          <  1:2214497‑M1/35‑11 (MQ=255)
          cAATTTGTGTCACCCGCAGCGCATAgaacgccatt                          <  1:1444465‑M1/35‑11 (MQ=255)
          cAATTTGTGTCACCCGCAGCGCATAgaacgccatt                          <  1:1414940‑M1/35‑11 (MQ=255)
           aaTTTGTGTCACCCGCAGCGCATAgaacgccattt                         >  1:5754990‑M1/1‑24 (MQ=255)
           aaTTTGTGTCACCCGCAGCGCATAgaacgccattt                         <  1:684355‑M1/35‑12 (MQ=255)
           aaTTTGTGTCACCCGCAGCGCATAgaacgccattt                         <  1:4225943‑M1/35‑12 (MQ=255)
           aaTTTGTGTCACCCGCAGCGCATAgaacgccattt                         <  1:3871202‑M1/35‑12 (MQ=255)
           aaTTTGTGTCACCCGCAGCGCATAgaacgccattt                         >  1:3674400‑M1/1‑24 (MQ=255)
           aaTTTGTGTCACCCGCAGCGCATAgaacgccattt                         >  1:7685000‑M1/1‑24 (MQ=255)
            aTTTGTGTCACCCGCAGCGCATAgaacgccctttt                        >  1:1257106‑M1/1‑23 (MQ=255)
            aTTTGTGTCACCCGCAGCGCATAgaacgccatttt                        >  1:379213‑M1/1‑23 (MQ=255)
            aTTTGTGTCACCCGCAGCGCATAgaacgccatttt                        >  1:1395986‑M1/1‑23 (MQ=255)
            aTTTGTGTCACCCGCAGCGCATAgaacgccatttt                        <  1:813417‑M1/35‑13 (MQ=255)
               tgtgTCACCCGCAGCGCATAgaacgccattttgaa                     >  1:6946695‑M1/1‑20 (MQ=255)
                 tgtCACCCGCAGCGCATAgaacgccattttgaatt                   >  1:7533047‑M1/1‑18 (MQ=255)
                 tgtCACCCGCAGCGCATAgaacgccattttgaatt                   <  1:1563780‑M1/35‑18 (MQ=255)
                 tgtCACCCGCAGCGCATAgaacgccattttgaatt                   <  1:4998667‑M1/35‑18 (MQ=255)
                 tgtCACCCGCAGCGCATAgaacgccattttgaatt                   >  1:4743449‑M1/1‑18 (MQ=255)
                 tgtCACCCGCAGCGCATAgaacgccattttgaatt                   <  1:7367554‑M1/35‑18 (MQ=255)
                 tgtCACCCGCAGCGCATAgaacgccattttgaatt                   >  1:3874998‑M1/1‑18 (MQ=255)
                  gtCACCCGCAGCGCATAgaacgccattttgaattt                  <  1:5156975‑M1/35‑19 (MQ=255)
                   tCACCCGCAGCGCATAgaacgccattttgaattta                 <  1:2200463‑M1/35‑20 (MQ=255)
                   tCACCCGCAGCGCATAgaacgccattttgaattta                 <  1:7281376‑M1/35‑20 (MQ=255)
                   tCACCCGCAGCGCATAgaacgccattttgaattta                 <  1:7699678‑M1/35‑20 (MQ=255)
                   tCACCCGCAGCGCATAgaacgccattttgaattta                 >  1:4952373‑M1/1‑16 (MQ=255)
                   tCACCCGCAGCGCATAgaacgccattttgaattta                 <  1:5256369‑M1/35‑20 (MQ=255)
                    cACCCGCAGCGCATAgaacgccattttgaatttac                <  1:4566036‑M1/35‑21 (MQ=255)
                    cACCCGCAGCGCATAgaacgccattttgaatttac                <  1:2761367‑M1/35‑21 (MQ=255)
                     ccccgCAGCGCATAgaacgccattttgaatttacg               <  1:4006310‑M1/34‑22 (MQ=255)
                     aCCCGCAGCGCATAgaacgccattttgaatttacg               >  1:5194722‑M1/1‑14 (MQ=255)
                       ccGCAGCGCATAgaacgccattttgaatttacgaa             <  1:2340608‑M1/35‑24 (MQ=255)
                       ccGCAGCGCATAgaacgccattttgaatttacgaa             <  1:7261797‑M1/35‑24 (MQ=255)
                        cGCAGCGCATAgaacgccattttgaatttacgaac            >  1:3612794‑M1/1‑11 (MQ=255)
                           aGCGCATAgaacgccattttgaatttacgaacgga         <  1:4514116‑M1/35‑28 (MQ=255)
                           aGCGCATAgaacgccattttgaatttacgaacgga         <  1:4319170‑M1/35‑28 (MQ=255)
                           aGCGCATAgaacgccattttgaatttacgaacgga         >  1:3871106‑M1/1‑8 (MQ=255)
                            gcgcATAgaacgccattttgaatttacgaacggac        <  1:1871050‑M1/35‑29 (MQ=255)
                            gcgcATAgaacgccattttgaatttacgaacggac        >  1:6326519‑M1/1‑7 (MQ=255)
                            gcgcATAgaacgccattttgaatttacgaacggac        >  1:1926519‑M1/1‑7 (MQ=255)
                             cgcATAgaacgccattttgaatttacgaacggacg       <  1:540342‑M1/35‑30 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     <  1:5646010‑M1/35‑32 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     >  1:7849053‑M1/1‑4 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     >  1:188514‑M1/1‑4 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     >  1:1886811‑M1/1‑4 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     >  1:213423‑M1/1‑4 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     >  1:2522335‑M1/1‑4 (MQ=255)
                               cATAgaacgccattttgaatttacgaacggacgca     >  1:5769323‑M1/1‑4 (MQ=255)
                                aTAgaacgccattttgaatttacgaacggacgcat    >  1:2814725‑M1/1‑3 (MQ=255)
                                aTAgaacgccattttgaatttacgaacggacgcat    <  1:4146081‑M1/35‑33 (MQ=255)
                                aTAgaacgccattttgaatttacgaacggacgcat    <  1:6457166‑M1/35‑33 (MQ=255)
                                 tAgaacgccattttgaatttacgaacggacgcatt   >  1:1954148‑M1/1‑2 (MQ=255)
                                 tAgaacgccattttgaatttacgaacggacgcatt   >  1:1106552‑M1/1‑2 (MQ=255)
                                 tAgaacgccattttgaatttacgaacggacgcatt   <  1:5388647‑M1/35‑34 (MQ=255)
                                 tAgaacgccattttgaatttacgaacggacgcatt   >  1:1343440‑M1/1‑2 (MQ=255)
                                  agaacgccattttgaatttacgaacggacgcatta  >  1:831398‑M1/1‑1 (MQ=255)
                                  agaacgccattttgaatttacgaacggacgcatta  >  1:1191668‑M1/1‑1 (MQ=255)
                                  |                                  
GCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGG  >  NC_000913/3110645‑3110713

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