Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,735,088 A→G G218G (GGA→GGG)  mepH → peptidoglycan DD‑endopeptidase MepH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,735,0880AG98.0% 173.8 / ‑5.0 51G218G (GGA→GGG) mepHpeptidoglycan DD‑endopeptidase MepH
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (31/19);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTTATAAAGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTC  >  NC_000913/1734961‑1735218
                                                                                                                               |                                                                                                                                  
cTTATAAAGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGTGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCgg                                                                                                                               <  2:741354/133‑1 (MQ=255)
cTTATAAAGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGTGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCgg                                                                                                                               >  1:741354/1‑133 (MQ=255)
cTTATAAAGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGcc                                                                                                                        >  2:821116/1‑140 (MQ=255)
cTTATAAAGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTCAAAAACGGCGACATGGCCTTTTTCCATACTCAGGGGCGCGGCACAGcc                                                                                                                        >  2:296351/1‑140 (MQ=255)
     aaaGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCa                                                                                                                   >  1:760293/1‑140 (MQ=255)
     aaaGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCa                                                                                                                   >  1:442672/1‑140 (MQ=255)
           ttGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCgg                                                                                                             >  1:110666/1‑140 (MQ=255)
             ggTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCg                                                                                                           >  1:664379/1‑140 (MQ=255)
                        gTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTcggc                                                                                                <  2:417122/140‑1 (MQ=255)
                             ccGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCTGCACAGCCGATCATGTCGGCGGGTATGTcggcaacgg                                                                                           >  1:396218/1‑140 (MQ=255)
                                cgTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTcggcaacggcaa                                                                                        >  1:415748/1‑140 (MQ=255)
                                cgTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTcggcaacggcaa                                                                                        >  1:522213/1‑140 (MQ=255)
                                           tGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGt                                                                             <  2:405597/140‑1 (MQ=255)
                                           tGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGt                                                                             <  2:176718/140‑1 (MQ=255)
                                                    tCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGca                                                                    <  2:107232/140‑1 (MQ=255)
                                                    tCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGca                                                                    <  2:760293/140‑1 (MQ=255)
                                                    tCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGca                                                                    <  2:730113/140‑1 (MQ=255)
                                                       ccTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGAATGGCGGCCACGGCAAATTTATTCAGTCACCGCGCACAg                                                                 >  2:378278/1‑140 (MQ=255)
                                                          gCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAgg                                                                >  1:450374/1‑138 (MQ=255)
                                                          gCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAgg                                                                <  2:450374/138‑1 (MQ=255)
                                                          gCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTc                                                              >  2:153117/1‑140 (MQ=255)
                                                           cGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCa                                                             >  1:573227/1‑140 (MQ=255)
                                                           cGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGACGGCAACTGCAAATTTATTCATTCACCGCGCACAGGTCa                                                             >  2:320893/1‑140 (MQ=255)
                                                            gTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGGGTGTATGTAGGCAACGGCAAATTTATTCAGGCACCGCGCACAGGTCAg                                                            >  2:269899/1‑140 (MQ=255)
                                                            gTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTGGGCAACGGCAAATTTAGTCAGCCACCGCGCACAGGTCAg                                                            >  2:262992/1‑140 (MQ=255)
                                                            gTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAg                                                            >  1:588426/1‑140 (MQ=255)
                                                            gTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAg                                                            >  1:573242/1‑140 (MQ=255)
                                                            gTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGAGTGTATGTAGGCAACGGCAAATTTACTCAGCCACCGCGCACAGGTCAg                                                            >  2:211601/1‑140 (MQ=255)
                                                                tGCAGCGCCAATCGACCGTAGTTAAATGACAAACGGCTACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGaaa                                                        <  1:262992/140‑1 (MQ=255)
                                                                    gcgcCAATCGAACGTATAGAACTGAAAAACCTCTACCTGGTCTTTTTCAGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCa                                                    <  2:396218/140‑1 (MQ=255)
                                                                    gcgcCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCa                                                    <  1:75258/140‑1 (MQ=255)
                                                                         aaTCGAACGTAGTGACCTGAAAGACGGCGTCCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCa                                               <  1:269899/140‑1 (MQ=255)
                                                                           tCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACt                                             <  1:821116/140‑1 (MQ=255)
                                                                                cGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTtctct                                        >  1:406034/1‑140 (MQ=255)
                                                                                      gAACTGAAAAACGGCAACCTGGTCTGTTTACGTACTCTGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGa                                  <  1:296351/140‑1 (MQ=255)
                                                                                        aCTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAg                                <  1:446000/140‑1 (MQ=255)
                                                                                           gAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTAACCAAATTCACATCACCTCTCTCAACGAAGact                             >  2:390966/1‑140 (MQ=255)
                                                                                            aaaaaCGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCa                                      <  1:533345/130‑1 (MQ=255)
                                                                                            aaaaaCGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCa                                      >  2:533345/1‑130 (MQ=255)
                                                                                            aaaaaCGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGacta                            >  1:166394/1‑140 (MQ=255)
                                                                                              aaaCGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGactact                          <  2:166394/140‑1 (MQ=255)
                                                                                               aaCGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTGAGGAAGTGCAAATCGCTTCTCGCAGCGAAGACTACTg                         >  1:255810/1‑140 (MQ=255)
                                                                                               aaCGGCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATGCAAGTCCCTTCTCCCAGCGAAGACTACTg                         >  1:403903/1‑140 (MQ=255)
                                                                                                   gCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAg                     <  2:406034/140‑1 (MQ=255)
                                                                                                   gCGACCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAg                     <  1:797626/140‑1 (MQ=255)
                                                                                                      aCCTGGTCTTTTTCCGTACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAgcgc                  >  1:731679/1‑140 (MQ=255)
                                                                                                              tttttCCGTACTCAGAGGCGCGACAAAGCTGATCATGCCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAGCGCCACTATGt          <  1:390966/140‑1 (MQ=255)
                                                                                                                      tACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCATAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAGCGCCACTATGTTGGCGCtc  >  2:790507/1‑140 (MQ=255)
                                                                                                                      tACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAGCGCCACTATGTTGGCGCtc  >  2:602982/1‑140 (MQ=255)
                                                                                                                      tACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGCGAAGACTACTGGCAGCGCCACTATGTTGGCGCtc  >  2:822512/1‑140 (MQ=255)
                                                                                                                      tACTCAGGGGCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTATCAGCGAAGACTACTGCCAGCGCCACTATGTTGGCGCtc  >  2:372617/1‑140 (MQ=255)
                                                                                                                               |                                                                                                                                  
CTTATAAAGATTTGGTGAAAATTCGTATTCCGCGTACGGCGAATGAAATGTATCACCTGCGTGATGCAGCGCCAATCGAACGTAGTGAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTC  >  NC_000913/1734961‑1735218

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