Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 95,730 A→G 32.4% M361V (ATG→GTG)  murF → D‑alanyl‑D‑alanine‑adding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091395,7300AG32.4% 38.6 / 29.3 46M361V (ATG→GTG) murFD‑alanyl‑D‑alanine‑adding enzyme
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base G (8/7);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

CCAGGCCGTCTGTTCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTATCAGCACCGCCAGCGGCGTTGGC  >  NC_000913/95589‑95861
                                                                                                                                             |                                                                                                                                   
gcAGGCCGGTTGTTCCTCCGCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGCCGGTTCAATGCCTGTAGCAGTCCAGGTACTGTCTTAAATGCCGTGCTACCGCGTGCTGGTGGTGGGCGatat                                                                                                                                    <  2:128999/142‑1 (MQ=255)
    gCCGTCTGTCCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACGCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACGGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCg                                                                                                                                <  1:391562/143‑1 (MQ=255)
     ccGTCTGTTCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCgg                                                                                                                               >  1:394857/1‑143 (MQ=255)
     ccGTCTGTTCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGAGGTGGCgg                                                                                                                               >  1:818126/1‑143 (MQ=255)
          tGTTCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACACCGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGGTGTAGCGGAACta                                                                                                                          >  2:306308/1‑142 (MQ=255)
            ttCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACGGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACAGCGTCCGGGGGGTAGGCGCTAGGGCgggcccggg                                                                                                                        >  2:765259/1‑136 (MQ=255)
               cccATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGgcgc                                                                                                                     <  1:187806/143‑1 (MQ=255)
                ccATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTATAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGGGGTGGGCGATGTGGCGGAGCTGGGCGCt                                                                                                                    >  1:688279/1‑143 (MQ=255)
                  accaaaCTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATAACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGa                                                                                                                  >  2:170638/5‑143 (MQ=255)
                   tCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGTAGCTGGGCGCTGaa                                                                                                                 >  1:874685/1‑143 (MQ=255)
                    ccAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCCGCAGTCCAGGTACTGGCTGAAATGCCGGGCTCCCGCGGGCGGGTGGTGGGCGAGAGGGCGGAGATGGGCGCTGaaa                                                                                                                >  1:78683/1‑143 (MQ=255)
                           acAGAAAAGCCGTTGCTGCTCGACGACTCCTACAGCGCCAATGTCGGTTCAATGGCTGCAGCAGTCCAGGGACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGc                                                                                                         <  2:715188/142‑1 (MQ=255)
                             aGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGGTATGGCGGAACCGGGCGCTGAAGGCGAAGACt                                                                                                       >  2:452587/1‑143 (MQ=255)
                                 aaCCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGGTGGTGGTGTTCGATATGGCGGAACAGGGCGCTGAACGAGAAGCGTGccg                                                                                                   >  1:324841/1‑142 (MQ=255)
                                    gaGTGGCTGAGCGACGACTGCGACAGCGCCAGGGTCGGTTCAAGGNCTGCGGCGGTCCAGGTAGTGGCTGAAGTGCCGGGCTACCGCGTGCTGGTGGTGGGCGAGATGGCGGAACGGGGCGCTGAAAGCGAAGCCTGccccgc                                                                                                <  2:229400/142‑5 (MQ=255)
                                          ctgctCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTAAAATGCAGAGCAACCGAGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGt                                                                                          >  1:770717/1‑143 (MQ=255)
                                            tctCGACGGCTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTgg                                                                                        <  1:127473/142‑1 (MQ=255)
                                             ctCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACTGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCAGGTGCAGGAggg                                                                                       >  1:93331/1‑143 (MQ=255)
                                                gacgacTCCTACAACGCCAATGTCGGTACAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGAGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGAGGGCGa                                                                                    >  2:785094/1‑143 (MQ=255)
                                                      tCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCGGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGGGGGCGAggcggc                                                                              >  1:227288/1‑143 (MQ=255)
                                                        ggtccaCGCAAATGTCGAGTCAGTGATGGCAGCAAGCCATGTACTGGCTGAAAGGCCGGGCTAGCGCGTGGGGGGGGTGGGCGATATGGCGGGAGTGGGGGCTGAAAGCGAAGCCGGCCATGTACAGGTGGGCGAGGCGGCaa                                                                            <  2:940064/138‑1 (MQ=255)
                                                         tccaccGCCTATCTCGGTTCACTGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCTCGTGCTGGTGGTGGTCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCaaa                                                                           <  2:22619/138‑1 (MQ=255)
                                                         caCAACGCCACTGCCGGTGCAATGACTGCAGCAGTCCAGGTACTCGCTGAAATGCCGCGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCaaa                                                                           <  1:279458/142‑1 (MQ=255)
                                                         caCAACGCCAATGTCGGTTCAATGACTGCAGCACTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCaaa                                                                           <  1:806822/142‑1 (MQ=255)
                                                          gccaCGCCAAGGTCGCTTCAATGACTGCAGCAGTCCAGGTACTGTCTGAAATGCCAGGGGACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCaaaa                                                                          <  2:345634/140‑1 (MQ=255)
                                                            aaCGCCAAGGTCGGTTCAGTGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAgc                                                                        <  1:578046/143‑1 (MQ=255)
                                                              cGCCAATGTCGGTTCAAGGACGGCAGCAGAGCAGGGACGGGCTGAAAGGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAgctg                                                                      <  1:72269/143‑1 (MQ=255)
                                                                  aaTGCCGGTCCAATGACTGCAGCAGTCCAGGTACGGGCTGAAATGCCGGGCGACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTgg                                                                  <  1:58501/143‑1 (MQ=255)
                                                                      tCGGTTGAATGCCTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTAtt                                                              <  1:479253/143‑1 (MQ=255)
                                                                                 aCTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGGGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGtt                                                   >  1:296456/1‑143 (MQ=255)
                                                                                   tGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTaa                                                 >  1:597173/1‑143 (MQ=255)
                                                                                   ggcagcaGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTaa                                                 <  1:732684/142‑1 (MQ=255)
                                                                                        cagTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCGGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTg                                            <  1:849400/143‑1 (MQ=255)
                                                                                        cagTCCAGGTACTGGCTGAAATGCCGCGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTg                                            <  1:627688/143‑1 (MQ=255)
                                                                                           tCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGCATTGACCGCGTGTTAAGCGTGGGt                                         >  2:548237/1‑143 (MQ=255)
                                                                                              aGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACAGCGGGTTAAGCGTGGGTaga                                      >  2:192326/1‑141 (MQ=255)
                                                                                                    tGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGc                                >  2:277280/1‑143 (MQ=255)
                                                                                                      gCTGGAATGCCGGGCTGCCGCGTGCGGGGGGTGGGCGATATGGCGGAACGGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCa                              <  2:458076/143‑1 (MQ=255)
                                                                                                      gCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTCAAGCGCGCGTCAACAAAGCCa                              >  2:981933/1‑143 (MQ=255)
                                                                                                         gAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTATTCGTGGGTAAACAAAGCCATGc                           >  1:272229/1‑143 (MQ=255)
                                                                                                             tGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTCCAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGCGGGTAAACAAAGCCATGCTATc                       >  1:688462/1‑143 (MQ=255)
                                                                                                             tGCCGGGCTACCGCGTGCTGGTGGTGGGCGATGTGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTATc                       >  1:941636/1‑143 (MQ=255)
                                                                                                                cGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTATCAGc                    >  1:214938/1‑143 (MQ=255)
                                                                                                                  ggCCACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTATCAGCAc                  <  1:729771/143‑1 (MQ=255)
                                                                                                                   gCTGCCGCGTTGTGGTGGTGGGCGATATGGCGGAACTGGGCGCGGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTATCAGCAcc                 <  2:770717/143‑1 (MQ=255)
                                                                                                                     tACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCAAGCTATCAGCACCGc               >  2:347507/1‑143 (MQ=255)
                                                                                                                                  tggtggGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTAGAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTCTCACCACCGCGACCGGCGTTCGc  >  1:987926/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                   
CCAGGCCGTCTGTTCCCCATCCAACTGGCAGAAAACCAGTTGCTGCTCGACGACTCCTACAACGCCAATGTCGGTTCAATGACTGCAGCAGTCCAGGTACTGGCTGAAATGCCGGGCTACCGCGTGCTGGTGGTGGGCGATATGGCGGAACTGGGCGCTGAAAGCGAAGCCTGCCATGTACAGGTGGGCGAGGCGGCAAAAGCTGCTGGTATTGACCGCGTGTTAAGCGTGGGTAAACAAAGCCATGCTATCAGCACCGCCAGCGGCGTTGGC  >  NC_000913/95589‑95861

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