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 44 [0] [0] 43 [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)41 (1.410)
+GA
35/92 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

CTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGC  >  NC_000913/3136335‑3136430
                                               |                                                
cttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                 <  1:1038846‑M2/2‑1 (MQ=255)
cttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                 >  1:1485546‑M2/48‑49 (MQ=255)
 ttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACg                                                >  1:1557689‑M2/47‑49 (MQ=255)
  tctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                               >  1:814548‑M2/46‑49 (MQ=255)
   ctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACg                                                <  1:1839725‑M2/3‑1 (MQ=255)
     cagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAt                                            >  1:1566361‑M2/43‑49 (MQ=255)
      agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGcc                                              <  1:30105‑M2/5‑1 (MQ=255)
      agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCa                                             <  1:152712‑M2/6‑1 (MQ=255)
      agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtt                                           <  1:450964‑M2/8‑1 (MQ=255)
      agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAt                                            <  1:1602731‑M2/7‑1 (MQ=255)
       gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAttt                                          <  1:1754035‑M2/9‑1 (MQ=255)
       gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAttt                                          <  1:2295770‑M2/9‑1 (MQ=255)
         agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                         >  1:650615‑M2/39‑48 (MQ=255)
         agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                         >  1:1131867‑M2/39‑48 (MQ=255)
           cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                      >  1:239140‑M2/37‑49 (MQ=255)
                acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAttt                                   <  1:2163004‑M2/16‑1 (MQ=255)
                acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                  <  1:1766880‑M2/17‑1 (MQ=255)
                acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                                 <  1:2155276‑M2/18‑1 (MQ=255)
                acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                                 <  1:1235168‑M2/18‑1 (MQ=255)
                   tcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGaa                              <  1:1040444‑M2/21‑1 (MQ=255)
                     aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAAc                             <  1:79266‑M2/22‑1 (MQ=255)
                     aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAAc                             <  1:2017760‑M2/22‑1 (MQ=255)
                     aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                            >  1:319058‑M2/27‑49 (MQ=255)
                      atttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACgg                           <  1:1409848‑M2/24‑1 (MQ=255)
                      atttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                            >  1:1775948‑M2/26‑48 (MQ=255)
                         tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                         <  1:1844267‑M2/26‑1 (MQ=255)
                         tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                         <  1:1432650‑M2/26‑1 (MQ=255)
                         tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                        >  1:2213571‑M2/23‑49 (MQ=255)
                          gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                        <  1:318513‑M2/27‑1 (MQ=255)
                             tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAtt                    <  1:1974274‑M2/31‑1 (MQ=255)
                             tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAtt                    <  1:657244‑M2/31‑1 (MQ=255)
                             tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                     <  1:920364‑M2/30‑1 (MQ=255)
                               acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                  >  1:2099799‑M2/17‑49 (MQ=255)
                               acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                  <  1:1559003‑M2/33‑1 (MQ=255)
                                cccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATc                 <  1:35457‑M2/34‑1 (MQ=255)
                                 ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                  >  1:92083‑M2/15‑47 (MQ=255)
                                 ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCa                <  1:337282‑M2/35‑1 (MQ=255)
                                     agcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGa            <  1:837298‑M2/39‑1 (MQ=255)
                                         catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt        <  1:2076178‑M2/43‑1 (MQ=255)
                                          atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt        >  1:1342453‑M2/6‑48 (MQ=255)
                                           tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa      <  1:1722528‑M2/45‑1 (MQ=255)
                                             gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc     <  1:1236487‑M2/46‑1 (MQ=255)
                                               aCGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGc  <  1:1794155/49‑1 (MQ=255)
                                               |                                                
CTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGC  >  NC_000913/3136335‑3136430

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