Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,681,810 A→G 9.8% S34G (AGC→GGC)  mltD → membrane‑bound lytic murein transglycosylase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,681,8100AG 9.8% 137.9 / 4.9 51S34G (AGC→GGC) mltDmembrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base G (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGCC  >  NC_002947/4681680‑4681953
                                                                                                                                   |                                                                                                                                               
aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                                >  1:329806/1‑148 (MQ=255)
         gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATa                                                                                                                      <  1:150731/149‑1 (MQ=255)
          cccGGATACCTGTATCAGCCGATGTCTTACCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCAAGAGCACCCGTCAGCTCGACGAATCCGATAg                                                                                                                     >  1:235017/1‑149 (MQ=255)
                   cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                            >  2:28304/1‑149 (MQ=255)
                      tATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGca                                                                                                         >  2:199933/1‑149 (MQ=255)
                                      cccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCC‑AATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                                           <  4:20533/146‑1 (MQ=255)
                                           aGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATc                                                                                      >  2:281816/1‑147 (MQ=255)
                                               gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGca                                                                                                         >  8:65464/1‑124 (MQ=255)
                                               gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGca                                                                                                         <  7:65464/124‑1 (MQ=255)
                                                   aaCCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAc                                                                            <  6:82876/149‑1 (MQ=255)
                                                    aCCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAc                                                                            <  2:329806/148‑1 (MQ=255)
                                                     ccTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCt                                                                          >  3:51525/1‑149 (MQ=255)
                                                      ctctCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGT‑GGCTGCCAGGGc                                                                                                                                               <  8:85736/79‑1 (MQ=18)
                                                      ctctCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGT‑GGCTGCCAGGGc                                                                                                                                               >  7:85736/1‑79 (MQ=18)
                                                      ctctCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                                       >  3:85175/1‑135 (MQ=255)
                                                      ctctCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                                       <  4:85175/135‑1 (MQ=255)
                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                                       <  3:52929/146‑1 (MQ=255)
                                                           aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCg                                                                    >  3:190695/1‑149 (MQ=255)
                                                            ttCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                                    >  3:123123/1‑132 (MQ=255)
                                                            ttCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                                    <  4:123123/132‑1 (MQ=255)
                                                               cGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCGCAAACCTTCCCCgctgc                                                                >  4:59003/1‑149 (MQ=255)
                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAcc                                                                                                 >  6:213099/1‑115 (MQ=255)
                                                                gTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAcc                                                                                                 <  5:213099/115‑1 (MQ=255)
                                                                            cgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                                   <  1:281816/149‑1 (MQ=255)
                                                                             gcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggcc                                                  <  4:190695/149‑1 (MQ=255)
                                                                               cTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggccg                                                 >  5:266974/1‑148 (MQ=255)
                                                                                   cccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAg                                            >  1:20689/1‑149 (MQ=255)
                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgcc                                                      >  1:167531/1‑138 (MQ=255)
                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgcc                                                      <  2:167531/138‑1 (MQ=255)
                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgcc                                                      <  2:172929/138‑1 (MQ=255)
                                                                                    ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgcc                                                      >  1:172929/1‑138 (MQ=255)
                                                                                      aaaTCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCaa                                                         <  1:88513/133‑1 (MQ=255)
                                                                                      aaaTCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCaa                                                         >  2:88513/1‑133 (MQ=255)
                                                                                                gcTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccgc                                     <  6:321904/143‑1 (MQ=255)
                                                                                                gcTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccgc                                     >  5:321904/1‑143 (MQ=255)
                                                                                                   tgctcTGGCAGCCCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt                            >  7:172966/5‑149 (MQ=255)
                                                                                                     ctcTGGCAGCCCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGCATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAGCCAGCACAAACCTTCCCAGCTGCTGGTCACGCCGGCCGAACAGGCGCCGCAGGATGTCTgg                          >  5:6202/3‑149 (MQ=255)
                                                                                                     cccTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCGCAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg                          <  3:59003/149‑1 (MQ=255)
                                                                                                       cTGGCAGCCCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa                        <  8:172966/149‑1 (MQ=255)
                                                                                                          gCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgcg                     <  2:20689/149‑1 (MQ=255)
                                                                                                            aG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCa                   <  1:28304/149‑1 (MQ=255)
                                                                                                                   cTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCa                                    >  1:75906/1‑126 (MQ=255)
                                                                                                                   cTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCa                                    <  2:75906/126‑1 (MQ=255)
                                                                                                                   cTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc              >  5:298048/1‑148 (MQ=255)
                                                                                                                    ttGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc              >  7:105175/1‑147 (MQ=255)
                                                                                                                      ggtggGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCa            <  2:235017/147‑1 (MQ=255)
                                                                                                                           gCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTc     <  6:266974/149‑1 (MQ=255)
                                                                                                                               ccAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAg                                                                                   <  4:230991/67‑1 (MQ=255)
                                                                                                                               ccAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAg           <  1:81940/139‑1 (MQ=255)
                                                                                                                               ccAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAg           >  2:81940/1‑139 (MQ=255)
                                                                                                                               ccAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGcc  <  4:51525/148‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                               
AACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGCC  >  NC_002947/4681680‑4681953

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: ‑