Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,468,078 Δ1 bp 9.1% coding (814/1047 nt) mltB ← membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,468,0780T. 9.1% 166.5 / 11.7 55coding (814/1047 nt)mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base . (2/3);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGC  >  NC_002947/5467938‑5468224
                                                                                                                                             |                                                                                                                                                  
gAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtga                                                                                                                                            <  2:43024/149‑2 (MQ=255)
 aaGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGACACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATCAACCGACCATCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtcgt                                                                                                                                           >  3:182391/1‑149 (MQ=255)
   gTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCggtggt                                                                                                                                           >  2:69798/1‑147 (MQ=255)
    ttcttcAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAg                                                                                                                                        <  2:23034/149‑1 (MQ=255)
         cAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGc                                                                                                                                      <  7:73806/146‑1 (MQ=255)
          aGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtc                                                                                                                                  <  5:285075/149‑1 (MQ=255)
            gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtc                                                                                                                                  >  4:8299/1‑147 (MQ=255)
            gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtc                                                                                                                                  <  3:8299/147‑1 (MQ=255)
            gCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGAACCATGCCCGTGGACACGCCTtcgt                                                                                                                                >  8:146770/1‑149 (MQ=255)
              ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtcg                                                                                                                              <  4:317295/149‑1 (MQ=255)
              ccATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTtcgtcg                                                                                                                              >  8:171554/1‑149 (MQ=255)
                aTCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgc                                                                                                                            >  6:40174/1‑149 (MQ=255)
                 tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                         >  4:89930/1‑131 (MQ=255)
                 tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                         >  4:218717/1‑131 (MQ=255)
                 tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                         <  3:89930/135‑5 (MQ=255)
                 tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCATCGAGCCggtgaaaa                                                                                                                                         <  3:218717/135‑5 (MQ=255)
                 tCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcg                                                                                                                           <  6:277058/149‑1 (MQ=255)
                   cAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcgcg                                                                                                                         >  6:218327/1‑149 (MQ=255)
                    aGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTAACGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGgcgcgc                                                                                                                        >  5:299334/1‑149 (MQ=255)
                            gggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGc                                                                                                                                                  <  1:139626/115‑1 (MQ=255)
                            gggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGc                                                                                                                                                  >  2:139626/1‑115 (MQ=255)
                             ggCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACTGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTcc                                                                                                               <  1:69798/149‑1 (MQ=255)
                               ccGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCAc                                                                                                             <  6:289366/149‑1 (MQ=255)
                                         ccTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGa                                                                                                   >  6:46584/1‑149 (MQ=255)
                                             gagaCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACt                                                                                                 <  8:86036/147‑1 (MQ=255)
                                              agaCGGAAGGCGGTAACCGGCAGATCATCGGGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGAc                                                                                              <  7:146770/149‑1 (MQ=255)
                                                aCGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACca                                                                                            >  2:100534/1‑149 (MQ=255)
                                                   gAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCacac                                                                                         >  1:196755/1‑149 (MQ=255)
                                                      ggCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCt                                                                                      >  2:258261/1‑149 (MQ=255)
                                                               cGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCTCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCa                                                                             <  2:85165/149‑1 (MQ=255)
                                                                 gCAGATCATCTCGCATCGAATCATGAACCGACCAGCCAATCCTACGCAACTCCCCTACCGTTTTCACCGGCTCGATCCCGGTGGTCAGGCCTTCCTCGGAGCGCTCACCTTACACCAATGCTCGAATAACCACACCCTCGCCGGCCAcc                                                                           >  8:305881/1‑149 (MQ=255)
                                                                   aGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCCCCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCa                                                                         >  6:256762/1‑149 (MQ=255)
                                                                         tcGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTg                                                                   >  2:77407/1‑149 (MQ=255)
                                                                                  gAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGaa                                                          <  1:100534/149‑1 (MQ=255)
                                                                                  gAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGaa                                                          <  7:171554/149‑1 (MQ=255)
                                                                                     tCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCATGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTa                                                       <  4:182391/149‑1 (MQ=255)
                                                                                       aTGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAg                                                      <  5:40174/148‑1 (MQ=255)
                                                                                          aaCCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTg                                                   <  6:299334/148‑1 (MQ=255)
                                                                                                   gCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                         <  3:180507/149‑1 (MQ=255)
                                                                                                    ccAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGcc                                        >  7:243435/1‑149 (MQ=255)
                                                                                                      aaGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGa                                      >  8:77047/1‑149 (MQ=255)
                                                                                                       aGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGAt                                     <  7:77047/149‑1 (MQ=255)
                                                                                                         cGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCGGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATgg                                   <  5:218327/149‑1 (MQ=255)
                                                                                                               cAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGt                             <  8:97489/149‑1 (MQ=255)
                                                                                                                aaCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCCGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg                                    >  7:314921/1‑141 (MQ=255)
                                                                                                                aaCTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCCGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATg                                    <  8:314921/141‑1 (MQ=255)
                                                                                                                     cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                         >  8:102580/1‑131 (MQ=255)
                                                                                                                     cccTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGc                                         <  7:102580/131‑1 (MQ=255)
                                                                                                                      ccTACCGTTTTCACCGGCTCGATGACGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTCACCACACCCTCGCCGCCCACCCAGCCGTGGCGAATGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGCgtt                      >  2:30880/1‑149 (MQ=255)
                                                                                                                       cTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGgttg                     >  1:339631/1‑149 (MQ=255)
                                                                                                                            gTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCa                <  2:196755/149‑1 (MQ=255)
                                                                                                                                  aCCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATAtt          <  2:339631/149‑1 (MQ=255)
                                                                                                                                   ccGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTg         <  1:77407/149‑1 (MQ=255)
                                                                                                                                    cGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACTCCCTCGCCGCCCACCCAGCCGTGGCGGTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGa        <  1:30880/149‑1 (MQ=255)
                                                                                                                                           gATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGc  >  6:120247/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GAAGTTCTTCAGGCCCATCCAGTACTCGGGGCCGCTCTCGCCCTCGAGACGGAAGGCGGTAACCGGCAGATCATCGCGCAGCGAATCATGAACCGACCAGCCAAGCGTACGCAACTCCCCTACCGTTTTCACCGGCTCGATGCCGGTGGTCAGGCCTTCGTCGGCGCGCTCACCTTCCACCCAGGCGCGACTGACCACACCCTCGCCGGCCACCCAGCCGTGGCGCTTGAAGTAGCTGGCCACGCTGCCGATGGCATCGTCAGGGTTGTTCCAGATATTGATATGGC  >  NC_002947/5467938‑5468224

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