Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,224,823 A→G A16A (GCT→GCC)  pbpG ← D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,224,8230AG97.7% 157.2 / NA 44A16A (GCT→GCC) pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/32);  minor base C (0/1);  total (11/33)
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.99e-01

ATCGGACGCACCAGATCCGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGC  >  NC_000913/2224677‑2224968
                                                                                                                                                  |                                                                                                                                                 
aTCGGACGCACCAGATCCGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCAcggcc                                                                                                                                                >  1:45947/1‑150 (MQ=255)
   ggACGCACCAGATCCGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCa                                                                                                                                            <  2:250480/151‑1 (MQ=255)
    gACGCACCAGATCCGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCAt                                                                                                                                           >  1:211497/1‑151 (MQ=255)
            aGATCCGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCa                                                                                                                                   <  2:273986/151‑1 (MQ=255)
                ccGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGc                                                                                                                                <  2:411316/150‑1 (MQ=255)
                   ggtggtTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCtaaa                                                                                                                            <  1:137061/151‑1 (MQ=255)
                      ggtTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCtaaataa                                                                                                                         <  2:297902/151‑1 (MQ=255)
                       gtTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCtaattaaa                                                                                                                        <  2:486082/151‑1 (MQ=255)
                         tCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAg                                                                                                                       <  2:703000/150‑1 (MQ=255)
                          cGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGaa                                                                                                                     <  1:558961/151‑1 (MQ=255)
                          cGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGaa                                                                                                                     <  2:636438/151‑1 (MQ=255)
                             aTAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAAc                                                                                                                   >  2:581593/1‑150 (MQ=255)
                                       ttgttgGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTc                                                                                                         <  1:706169/150‑1 (MQ=255)
                                           tgGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCa                                                                                                     <  2:283431/150‑1 (MQ=255)
                                                ttCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCatgatg                                                                                                >  1:228783/1‑150 (MQ=255)
                                                 tCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCAtgatga                                                                                               <  2:149964/150‑1 (MQ=255)
                                                         acaaTCATCGCGCTACCGGAGGCAAGTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGa                                                                                      <  2:538855/151‑1 (MQ=255)
                                                          caaTCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGa                                                                                      <  1:587657/150‑1 (MQ=255)
                                                              catcGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTg                                                                                 >  2:291238/1‑151 (MQ=255)
                                                                  gcgcTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTc                                                                             <  1:651619/151‑1 (MQ=255)
                                                                       aCCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCAcgcg                                                                        >  1:731250/1‑151 (MQ=255)
                                                                              gCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGt                                                                 <  2:731206/151‑1 (MQ=255)
                                                                               cAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGtt                                                                <  1:259107/151‑1 (MQ=255)
                                                                                aaTTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTg                                                               <  2:394692/151‑1 (MQ=255)
                                                                                    tCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTgcgc                                                            >  2:459029/1‑150 (MQ=255)
                                                                                     ccGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCa                                                           <  1:13248/150‑1 (MQ=255)
                                                                                     ccGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCAc                                                          <  1:48233/151‑1 (MQ=255)
                                                                                        gTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCAccgc                                                       >  1:220386/1‑151 (MQ=255)
                                                                                         ttGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCAccgcc                                                      >  2:517220/1‑151 (MQ=255)
                                                                                                gCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTa                                                >  1:441598/1‑150 (MQ=255)
                                                                                                 cGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACGCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAg                                              <  2:738538/151‑1 (MQ=255)
                                                                                                    tggtAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGAt                                           >  2:731317/1‑151 (MQ=255)
                                                                                                           ggCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGa                                    <  1:612584/151‑1 (MQ=255)
                                                                                                                 cGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACCGCCAGCATCAGGGCCAGGCTACATAAAGAAACTCGAAATTTCCGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTCCTGAGGCTAg                              <  2:231863/151‑1 (MQ=255)
                                                                                                                    tttCGCAACGGCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCg                            >  2:59602/1‑150 (MQ=255)
                                                                                                                              gCCTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATa                 <  2:72501/151‑1 (MQ=255)
                                                                                                                               ccTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATaa                <  2:658700/151‑1 (MQ=255)
                                                                                                                               ccTGCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATaa                <  1:554784/151‑1 (MQ=255)
                                                                                                                                  gCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACg              <  1:724923/150‑1 (MQ=255)
                                                                                                                                  gCGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGa             <  2:45946/151‑1 (MQ=255)
                                                                                                                                   cGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGa             <  1:763728/150‑1 (MQ=255)
                                                                                                                                   cGGTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAg            <  1:315344/151‑1 (MQ=255)
                                                                                                                                    ggTGCAAAAGGCAAGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCGTCATAACGAGc           <  2:501911/151‑1 (MQ=255)
                                                                                                                                     gTGCAAAAGGCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCa          <  2:675434/151‑1 (MQ=255)
                                                                                                                                             ggCACGGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTTCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGc  <  2:657537/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCGGACGCACCAGATCCGGGTGGTTCGAATAGATCACTTTGTTGGTATTCAGATCAACAATCATCGCGCTACCGGAGGCAATTTCCGGTTGTGAAGCGGTGGTAGCGGCTGCCGTTTTCGCAACGGCCTGCGGTGCAAAAGGCACAGCCAGCATCAGGGCCAGGCTAAATAAAGAAACTCGAAATTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGC  >  NC_000913/2224677‑2224968

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