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

AAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCA  >  NC_000913/3136348‑3136416
                                  |                                  
ctgacatcaatttgtgtcacccgcagcgcatagaa                                    >  1:6497918‑M2/35‑35 (MQ=255)
ctgacatcaatttgtgtcacccgcagcgcatagaa                                    >  1:4613192‑M2/35‑35 (MQ=255)
ctgacatcaatttgtgtcacccgcagcgcatagaa                                    >  1:4448670‑M2/35‑35 (MQ=255)
ctgacatcaatttgtgtcacccgcagcgcatagaa                                    <  1:2154487‑M2/1‑1 (MQ=255)
 tgacatcaatttgtgtcacccgcagcgcatagaAc                                   >  1:981445‑M2/34‑35 (MQ=255)
 tgacatcaatttgtgtcacccgcagcgcatagaAc                                   >  1:1194208‑M2/34‑35 (MQ=255)
   acatcaatttgtgtcacccgcagcgcatagaACGc                                 <  1:7371379‑M2/4‑1 (MQ=255)
   acatcaatttgtgtcacccgcagcgcatagaACGc                                 <  1:5830421‑M2/4‑1 (MQ=255)
    catcaatttgtgtcacccgcagcgcatagaACGcc                                <  1:416532‑M2/5‑1 (MQ=255)
     atcaatttgtgtcacccgcagcgcatagaACGCCa                               <  1:4840402‑M2/6‑1 (MQ=255)
      tcaatttgtgtcacccgcagcgcatagaACGCCAt                              <  1:5598225‑M2/7‑1 (MQ=255)
       caatttgtgtcacccgcagcgcatagaACGCCAtt                             <  1:4680584‑M2/8‑1 (MQ=255)
       caatttgtgtcacccgcagcgcatagaACGCCAtt                             <  1:2214497‑M2/8‑1 (MQ=255)
       caatttgtgtcacccgcagcgcatagaACGCCAtt                             <  1:1414940‑M2/8‑1 (MQ=255)
       caatttgtgtcacccgcagcgcatagaACGCCAtt                             <  1:1444465‑M2/8‑1 (MQ=255)
        aatttgtgtcacccgcagcgcatagaACGCCAttt                            >  1:7685000‑M2/27‑35 (MQ=255)
        aatttgtgtcacccgcagcgcatagaACGCCAttt                            >  1:5754990‑M2/27‑35 (MQ=255)
        aatttgtgtcacccgcagcgcatagaACGCCAttt                            <  1:4225943‑M2/9‑1 (MQ=255)
        aatttgtgtcacccgcagcgcatagaACGCCAttt                            <  1:3871202‑M2/9‑1 (MQ=255)
        aatttgtgtcacccgcagcgcatagaACGCCAttt                            >  1:3674400‑M2/27‑35 (MQ=255)
        aatttgtgtcacccgcagcgcatagaACGCCAttt                            <  1:684355‑M2/9‑1 (MQ=255)
         atttgtgtcacccgcagcgcatagaACGCCCtttt                           >  1:1257106‑M2/26‑35 (MQ=255)
         atttgtgtcacccgcagcgcatagaACGCCAtttt                           <  1:813417‑M2/10‑1 (MQ=255)
         atttgtgtcacccgcagcgcatagaACGCCAtttt                           >  1:1395986‑M2/26‑35 (MQ=255)
         atttgtgtcacccgcagcgcatagaACGCCAtttt                           >  1:379213‑M2/26‑35 (MQ=255)
            tgtgtcacccgcagcgcatagaACGCCATTTTGaa                        >  1:6946695‑M2/23‑35 (MQ=255)
              tgtcacccgcagcgcatagaACGCCATTTTGAAtt                      <  1:4998667‑M2/15‑1 (MQ=255)
              tgtcacccgcagcgcatagaACGCCATTTTGAAtt                      >  1:3874998‑M2/21‑35 (MQ=255)
              tgtcacccgcagcgcatagaACGCCATTTTGAAtt                      >  1:4743449‑M2/21‑35 (MQ=255)
              tgtcacccgcagcgcatagaACGCCATTTTGAAtt                      <  1:1563780‑M2/15‑1 (MQ=255)
              tgtcacccgcagcgcatagaACGCCATTTTGAAtt                      <  1:7367554‑M2/15‑1 (MQ=255)
              tgtcacccgcagcgcatagaACGCCATTTTGAAtt                      >  1:7533047‑M2/21‑35 (MQ=255)
               gtcacccgcagcgcatagaACGCCATTTTGAAttt                     <  1:5156975‑M2/16‑1 (MQ=255)
                tcacccgcagcgcatagaACGCCATTTTGAATTTa                    <  1:7281376‑M2/17‑1 (MQ=255)
                tcacccgcagcgcatagaACGCCATTTTGAATTTa                    >  1:4952373‑M2/19‑35 (MQ=255)
                tcacccgcagcgcatagaACGCCATTTTGAATTTa                    <  1:5256369‑M2/17‑1 (MQ=255)
                tcacccgcagcgcatagaACGCCATTTTGAATTTa                    <  1:7699678‑M2/17‑1 (MQ=255)
                tcacccgcagcgcatagaACGCCATTTTGAATTTa                    <  1:2200463‑M2/17‑1 (MQ=255)
                 cacccgcagcgcatagaACGCCATTTTGAATTTAc                   <  1:2761367‑M2/18‑1 (MQ=255)
                 cacccgcagcgcatagaACGCCATTTTGAATTTAc                   <  1:4566036‑M2/18‑1 (MQ=255)
                  ccccgcagcgcatagaACGCCATTTTGAATTTACg                  <  1:4006310‑M2/19‑1 (MQ=255)
                  acccgcagcgcatagaACGCCATTTTGAATTTACg                  >  1:5194722‑M2/17‑35 (MQ=255)
                    ccgcagcgcatagaACGCCATTTTGAATTTACGaa                <  1:2340608‑M2/21‑1 (MQ=255)
                    ccgcagcgcatagaACGCCATTTTGAATTTACGaa                <  1:7261797‑M2/21‑1 (MQ=255)
                     cgcagcgcatagaACGCCATTTTGAATTTACGAAc               >  1:3612794‑M2/14‑35 (MQ=255)
                        agcgcatagaACGCCATTTTGAATTTACGAACGGa            <  1:4514116‑M2/25‑1 (MQ=255)
                        agcgcatagaACGCCATTTTGAATTTACGAACGGa            <  1:4319170‑M2/25‑1 (MQ=255)
                        agcgcatagaACGCCATTTTGAATTTACGAACGGa            >  1:3871106‑M2/11‑35 (MQ=255)
                         gcgcatagaACGCCATTTTGAATTTACGAACGGAc           >  1:6326519‑M2/10‑35 (MQ=255)
                         gcgcatagaACGCCATTTTGAATTTACGAACGGAc           >  1:1926519‑M2/10‑35 (MQ=255)
                         gcgcatagaACGCCATTTTGAATTTACGAACGGAc           <  1:1871050‑M2/26‑1 (MQ=255)
                          cgcatagaACGCCATTTTGAATTTACGAACGGACg          <  1:540342‑M2/27‑1 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        >  1:2522335‑M2/7‑35 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        >  1:1886811‑M2/7‑35 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        >  1:7849053‑M2/7‑35 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        <  1:5646010‑M2/29‑1 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        >  1:5769323‑M2/7‑35 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        >  1:188514‑M2/7‑35 (MQ=255)
                            catagaACGCCATTTTGAATTTACGAACGGACGCa        >  1:213423‑M2/7‑35 (MQ=255)
                             atagaACGCCATTTTGAATTTACGAACGGACGCAt       >  1:2814725‑M2/6‑35 (MQ=255)
                             atagaACGCCATTTTGAATTTACGAACGGACGCAt       <  1:6457166‑M2/30‑1 (MQ=255)
                             atagaACGCCATTTTGAATTTACGAACGGACGCAt       <  1:4146081‑M2/30‑1 (MQ=255)
                              tagaACGCCATTTTGAATTTACGAACGGACGCAtt      >  1:1954148‑M2/5‑35 (MQ=255)
                              tagaACGCCATTTTGAATTTACGAACGGACGCAtt      >  1:1343440‑M2/5‑35 (MQ=255)
                              tagaACGCCATTTTGAATTTACGAACGGACGCAtt      >  1:1106552‑M2/5‑35 (MQ=255)
                              tagaACGCCATTTTGAATTTACGAACGGACGCAtt      <  1:5388647‑M2/31‑1 (MQ=255)
                               agaACGCCATTTTGAATTTACGAACGGACGCATTa     >  1:831398‑M2/4‑35 (MQ=255)
                               agaACGCCATTTTGAATTTACGAACGGACGCATTa     >  1:1191668‑M2/4‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:3075590‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    <  1:7786730‑M2/33‑1 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:1194021‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:2354905‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:5025140‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:6228451‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:5859725‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    >  1:4801310‑M2/3‑35 (MQ=255)
                                gaACGCCATTTTGAATTTACGAACGGACGCATTAt    <  1:4953441‑M2/33‑1 (MQ=255)
                                 aACGCCATTTTGAATTTACGAACGGACGCATTATc   <  1:7238218‑M2/34‑1 (MQ=255)
                                 aACGCCATTTTGAATTTACGAACGGACGCATTATc   >  1:7084439‑M2/2‑35 (MQ=255)
                                 aACGCCATTTTGAATTTACGAACGGACGCATTATc   <  1:1547804‑M2/34‑1 (MQ=255)
                                  aCGCCATTTTGAATTTACGAACGGACGCATTATCa  >  1:287750/1‑35 (MQ=255)
                                  aCGCCATTTTGAATTTACGAACGGACGCATTATCa  >  1:5775665/1‑35 (MQ=255)
                                  aCGCCATTTTGAATTTACGAACGGACGCATTATCa  >  1:562434/1‑35 (MQ=255)
                                  aCGCCATTTTGAATTTACGAACGGACGCATTATCa  <  1:5351904/35‑1 (MQ=255)
                                  |                                  
AAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCA  >  NC_000913/3136348‑3136416

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