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 73 [0] [0] 70 [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)67 (0.960)
+GA
47/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

CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGG  >  NC_000913/3136333‑3136429
                                                 |                                               
ggcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaa                                                 >  1:5214048‑M2/50‑50 (MQ=255)
ggcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaa                                                 <  1:1079727‑M2/1‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                <  1:1256885‑M2/2‑1 (MQ=255)
   ttctcagcaacactgacatcaatttgtgtcacccgcagcgcatagaACGc                                              >  1:3571503‑M2/47‑50 (MQ=255)
    tctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                              >  1:1495326‑M2/46‑49 (MQ=255)
    tctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                              >  1:701592‑M2/46‑49 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAttt                                         <  1:4024431‑M2/9‑1 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                       <  1:2604080‑M2/11‑1 (MQ=255)
           agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                       >  1:3423153‑M2/39‑49 (MQ=255)
           agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                       >  1:2992601‑M2/39‑49 (MQ=255)
            gcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                     >  1:1885183‑M2/38‑50 (MQ=255)
            gcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                     >  1:730414‑M2/38‑50 (MQ=255)
            gcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                     <  1:2926216‑M2/13‑1 (MQ=255)
            gcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                     >  1:2452031‑M2/38‑50 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                    >  1:3751208‑M2/37‑50 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAttt                                  <  1:4383385‑M2/16‑1 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                   >  1:3064552‑M2/35‑49 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                   >  1:3191874‑M2/35‑49 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                 <  1:2466864‑M2/17‑1 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                                >  1:3990504‑M2/32‑49 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                                <  1:5024461‑M2/18‑1 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                               <  1:540241‑M2/19‑1 (MQ=255)
                   catcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                               >  1:2658558‑M2/31‑49 (MQ=255)
                   catcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                               >  1:2678506‑M2/31‑49 (MQ=255)
                    atcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGaa                             <  1:2312868‑M2/21‑1 (MQ=255)
                       aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                           >  1:4545385‑M2/27‑49 (MQ=255)
                       aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                           >  1:566319‑M2/27‑49 (MQ=255)
                        atttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGa                         <  1:2998150‑M2/25‑1 (MQ=255)
                        atttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGa                         >  1:5013220‑M2/26‑50 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                        <  1:1953116‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                        <  1:1806314‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                        <  1:378618‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                        >  1:3151031‑M2/25‑50 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                       >  1:693006‑M2/23‑49 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                       >  1:695988‑M2/23‑49 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCa                     >  1:2166625‑M2/22‑50 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCa                     >  1:4554299‑M2/22‑50 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCa                     >  1:1780405‑M2/22‑50 (MQ=255)
                             tgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                    >  1:5055269‑M2/21‑50 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAtt                   <  1:4194806‑M2/31‑1 (MQ=255)
                                cacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                 >  1:743289‑M2/18‑50 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                 <  1:5513341‑M2/33‑1 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                 >  1:1941216‑M2/17‑49 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                 <  1:4426230‑M2/33‑1 (MQ=255)
                                  cccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCa               <  1:3561058‑M2/35‑1 (MQ=255)
                                   ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAg              <  1:942953‑M2/36‑1 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt             <  1:477892‑M2/37‑1 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt             >  1:3109838‑M2/14‑50 (MQ=255)
                                      cagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGa           <  1:4823540‑M2/39‑1 (MQ=255)
                                       agcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGAc          >  1:3177520‑M2/11‑50 (MQ=255)
                                        gcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACt         <  1:4878667‑M2/41‑1 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt        >  1:4906506‑M2/9‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt        >  1:3952468‑M2/9‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt        >  1:2537736‑M2/9‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt        >  1:5226974‑M2/9‑50 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt       <  1:3356822‑M2/43‑1 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt        >  1:1318424‑M2/8‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa      <  1:2277648‑M2/44‑1 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa      <  1:4216709‑M2/44‑1 (MQ=255)
                                            atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa     <  1:4406949‑M2/45‑1 (MQ=255)
                                            atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa     >  1:3482907‑M2/6‑50 (MQ=255)
                                            atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa     <  1:2280738‑M2/45‑1 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc    <  1:514654‑M2/46‑1 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc    >  1:5411700‑M2/5‑50 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc    <  1:5024352‑M2/46‑1 (MQ=255)
                                              agaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg   >  1:271762‑M2/4‑50 (MQ=255)
                                              agaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg   >  1:1525793‑M2/4‑50 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg  >  1:1293237‑M2/3‑50 (MQ=255)
                                                aACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg  >  1:3286583‑M2/2‑49 (MQ=255)
                                                aACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg  >  1:431086‑M2/2‑49 (MQ=255)
                                                 |                                               
CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGG  >  NC_000913/3136333‑3136429

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