Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,349,538 T→C H91R (CAC→CGC)  mtgA ← peptidoglycan glycosyltransferase MtgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,349,5380TC97.4% 129.0 / ‑4.1 38H91R (CAC→CGC) mtgApeptidoglycan glycosyltransferase MtgA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/12);  total (26/12)
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.50e-01

AACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCTGACCTGTCGCTCGA  >  NC_000913/3349402‑3349673
                                                                                                                                        |                                                                                                                                       
aaCCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCt                                                                                                                                      <  2:464661/140‑1 (MQ=255)
    cAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAgg                                                                                                                                  <  1:1252837/140‑1 (MQ=255)
     aGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGa                                                                                                                                 <  1:30574/140‑1 (MQ=255)
      gCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGaa                                                                                                                                >  2:427859/1‑140 (MQ=255)
              catcccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGa                                                                                                                        <  2:1367315/140‑1 (MQ=255)
               atcccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCTACATCAAAGCCCCAGCGCTCAGGAAAttt                                                                                                                            <  1:1140133/135‑1 (MQ=255)
               atcccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCTACATCAAAGCCCCAGCGCTCAGGAAAttt                                                                                                                            >  2:1140133/1‑135 (MQ=255)
                 cccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGAtct                                                                                                                     >  1:695215/1‑140 (MQ=255)
                 cccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGAtct                                                                                                                     >  1:545174/1‑140 (MQ=255)
                 cccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGAtct                                                                                                                     >  2:927526/1‑140 (MQ=255)
                        aaaGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCgg                                                                                                              <  1:427859/140‑1 (MQ=255)
                                 tttGGCTGTTTGTGGAGAAATCTTTGAAGCCCCGCGAATGCGGTTTTCATTGCGCTCGTTGTGCGCGGGGGTTTTCTCAATGGAAGCGACATCAAACCCCCGGCGCTCAGAAATTTTCTGATCTTCTGCGGCAATCACCg                                                                                                     <  2:1400700/140‑1 (MQ=255)
                                 tttGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCAcca                                                                                                     <  1:1057465/140‑2 (MQ=255)
                                 tttGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCg                                                                                                     >  1:406488/1‑140 (MQ=255)
                                    ggCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATc                                                                                                         <  2:1275665/133‑1 (MQ=255)
                                    ggCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATc                                                                                                         >  1:1275665/1‑133 (MQ=255)
                                        gtctgtTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTcc                                                                                              >  1:514590/1‑140 (MQ=255)
                                         tctgtTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTccc                                                                                             >  2:612728/1‑140 (MQ=255)
                                                   aaTCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGc                                                                                                    >  2:69556/1‑123 (MQ=255)
                                                   aaTCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGc                                                                                                    <  1:69556/123‑1 (MQ=255)
                                                     tCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCTTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTACGGCAATCACCGCCAGTCCCATCCACGGCgag                                                                                 >  2:1255288/1‑140 (MQ=255)
                                                     tCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCgag                                                                                 >  1:427062/1‑140 (MQ=255)
                                                       gtggAAGCCCCGCGAATACGGTTTTCATGGCGTTGGTTGCGCGCCAGGGCTTTCTCAAGGAAAGCAACATAAAACCCCCAGCGCTCAGGAATTTTCTGTTCTTCTGGGCCAATCCCCGCCAGTCCCATCCACGGCGAGAt                                                                               <  2:406488/137‑1 (MQ=255)
                                                       gTTGAAGCCCCGGGAAATGGGTTTTCATTGGCCTCGTTGTGGCCCAGGGTTTTCCCAAGGGAAGGCAAATAAAAGCCCCAACCCCAAGGAAATTTCGGAACTTCTGGGGGAATCCCCGCCAGTCCCATCCACGGCGAGAt                                                                               <  2:882351/140‑1 (MQ=255)
                                                         taaacccccGCAAATACGTTTTTCATTGCGCTCGTTGTGCGCCGGGGTTTTCTAAATGGAAGCGACATCAAACCCCGAGCGCTCAGGAAATTTCTGTTCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGAttt                                                                             <  2:545174/133‑1 (MQ=255)
                                                          gAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTg                                                                            >  2:796489/1‑140 (MQ=255)
                                                          gAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTg                                                                            >  1:937267/1‑140 (MQ=255)
                                                            aGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGAt                                                                          <  1:1406999/140‑1 (MQ=255)
                                                               cccGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCa                                                                       <  1:796489/139‑1 (MQ=255)
                                                                     aaTACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGTTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGa                                                                 <  1:41829/140‑1 (MQ=255)
                                                                       tACGTTTTTCATTGCGCTCGTTGTGCGCCAGGGTTTTCTAAATGGAAGCGACATAAAACCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGAcc                                                               <  2:427062/140‑1 (MQ=255)
                                                                           tttttCATTGCGCTCGTTGTGCGCCGGGTTTTTCTAAATGGAAGCGACAAAAAACCCCCAGCGCTCAGAAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGt                                                           <  2:308917/139‑1 (MQ=255)
                                                                           gTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGt                                                           >  2:1394722/1‑140 (MQ=255)
                                                                           gTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGt                                                           >  1:308917/1‑140 (MQ=255)
                                                                               tCATTGCGCTCGTTGTGCGCCGGGGTTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTAAGAAATTTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGa                                                       <  1:1255288/140‑1 (MQ=255)
                                                                               tCATGGCGCTCGTTGTGCGCCGGGGTTTTCTCAAGGGAAGCGACACAAAACCCCCAGCGCTCAGAAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGa                                                       <  2:514590/140‑1 (MQ=255)
                                                                                    gcgcTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAAtgtg                                                  >  1:595466/1‑140 (MQ=255)
                                                                                    gcgcTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAACCCCCAGCGCTCAGAAAATTTCTGATCTTCTGCGGCAACCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAAtgtg                                                  <  2:15745/140‑1 (MQ=255)
                                                                                     ccgtCGTTGTGCGCGAGGGCTTTCTCAATGGAAGCTACATAAAAGCCCCCGCGCGGAGGAAATTTCTGACCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGc                                                 <  2:1377869/137‑1 (MQ=255)
                                                                                           ttGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTa                                           >  2:1153670/1‑140 (MQ=255)
                                                                                                     gggCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTg                                 >  1:751800/1‑140 (MQ=255)
                                                                                                      ggCTTTCTCAATGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGc                                >  2:1247331/1‑140 (MQ=255)
                                                                                                       gCTTTCTCAAGGGAAGCAACATAAAACCCCCAGCGCTCAGAAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGcc                               <  2:937267/140‑1 (MQ=255)
                                                                                                               aagggAAGCGACATCAAACCCCCAGCGCTCAGAAAATTTCTGATCTTCGGCGGCATCCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGcc                       <  1:927526/137‑1 (MQ=255)
                                                                                                               aaTGGAAGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGcc                       >  2:1029998/1‑140 (MQ=255)
                                                                                                                     aGCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGgcgc                 >  1:413643/1‑140 (MQ=255)
                                                                                                                      gCGACATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCt                >  1:487768/1‑140 (MQ=255)
                                                                                                                         aCATCAAAGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCTGAc             >  1:544964/1‑140 (MQ=255)
                                                                                                                               aaGCCCCAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCTGACCTGTCg       >  1:1032166/1‑140 (MQ=255)
                                                                                                                                 ccccaAGCGCTAAGGAAATTTCTGATCTTCGGCGGCAATCGCCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAAGCGCTGACCTGTCGCt     <  1:276530/139‑1 (MQ=255)
                                                                                                                                  ccccAGCGCTAAGGAAATTTCTGATTTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCTGACCTGTCGCTc    <  1:639814/140‑1 (MQ=255)
                                                                                                                                    ccAGCGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCTGACCTGTCGCTCGa  >  1:822926/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACCCAGTCAGAATGTGCCACGTAACGAAAATTGCCATGCAGCCAGGCGCTGACCTGTCGCTCGA  >  NC_000913/3349402‑3349673

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