Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,940,584 C→A R351L (CGT→CTT)  mepM ← peptidoglycan DD‑endopeptidase MepM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,940,5840CA91.1% 122.2 / NA 47R351L (CGT→CTT) mepMpeptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (21/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGATT  >  NC_000913/1940445‑1940723
                                                                                                                                            |                                                                                                                                           
ttCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGCTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGCTGTAGCTGCGACCATGCCGaa                                                                                                                                          >  1:212425/1‑143 (MQ=255)
 tCAGAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAAt                                                                                                                                         <  1:724397/143‑1 (MQ=255)
       tGCAGATGCGGCCCGGTTGAACGTCCGGTATTATCGGAAAGCGCGATACGGTCGCCACGTTTCACCATCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAGGAATAGCcac                                                                                                                                   >  1:677916/1‑143 (MQ=255)
         cAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACat                                                                                                                                 <  1:788666/143‑1 (MQ=255)
           gATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGGGCGGTACGGTAGGGACGAAACAAAGAATGCAATCTGCTGACGATGAGCTTGCGCAAGTGCTTATCACGCGTGGTCTAGCTGCGACCATGCAGAATAGCCACataa                                                                                                                               >  1:241052/1‑143 (MQ=255)
           gATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACataa                                                                                                                               <  1:178177/143‑1 (MQ=255)
            aTGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGCATCTTGCGCAAGTGCATCTAACGCGTGGTGGAGCAGCGCCCATGGACGATACCCACataag                                                                                                                              >  2:761399/1‑142 (MQ=255)
              gCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGGCCCGGTTTGACCAGAATCTTGCGCAAGGGCATATAACGCGTGGTGTAGCTGCGCCCATGAAGAATAGCCACAtaataa                                                                                                                            >  1:730551/1‑143 (MQ=255)
                ggCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAAcc                                                                                                                          <  1:785579/143‑1 (MQ=255)
                 ccccGTTTGACCGTCCGGTCTTACCGGAACGCGCGTGACGGTCGCCACGTTTCACCTTCTGTCCCGATTTCACCAGATTCTTGCGCAAGTGCATATAACTCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACct                                                                                                                         <  2:871289/142‑1 (MQ=255)
                    cGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGGGTGGGGTAGCTGCGGCCATGAAGAATAGCCACAAGATAACctgct                                                                                                                      >  1:194975/1‑143 (MQ=255)
                              tCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGc                                                                                                            >  1:728354/1‑143 (MQ=255)
                                  ggtttACCGGAAGGCGCGGTACGGTCGCCGCGTTTCACCGTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGttt                                                                                                        <  2:321093/140‑1 (MQ=255)
                                    aTTACCGGACAGCGCCATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGCATCTCGCGCAAGTCCATATAACGCGTGGTGTAGATGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTgg                                                                                                      <  1:73495/143‑1 (MQ=255)
                                        ccGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAac                                                                                                  <  2:871715/143‑1 (MQ=255)
                                         cGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCCCCAGAATCTTGCGCAAGTGCATATACCCCGTGGTGTCGCCGCGACCACGCACCCTAGCCACCTAATCACCTGCTGCCCCACTCCGTCTCGCGacc                                                                                                 >  1:886218/1‑143 (MQ=255)
                                         cGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCACATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCGGCTGCGCCGCTGCGGCTCGCAacc                                                                                                 >  2:424803/1‑143 (MQ=255)
                                         cGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTATCCCGGTTTCACCAGAATCTTGCGGAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTCAGGTTGGCCacc                                                                                                 >  1:271094/1‑143 (MQ=255)
                                                 gcgATACGGTCGCCACGTTTCACCTTGTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTAGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTc                                                                                         <  2:668644/143‑1 (MQ=255)
                                                                  tttCACCTTCTGTCCCTGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGaaa                                                                        <  1:542557/143‑1 (MQ=255)
                                                                  tttCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAACTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTCGCAACCCCCACTCCACCGTCCCCCCCTGaca                                                                        >  2:448564/1‑141 (MQ=255)
                                                                  tttCACCTTCTGTCCCGGTTTCACCACACTCTTGCGCAACTGCCTATAACGCCTGGTGTAGCTGCCACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACCGCGTCTGGCACCCACCACTTCCCCCTCACCCCCTgcca                                                                        >  2:85079/1‑140 (MQ=255)
                                                                  tttCACCGTCTGTCCCGGTTTCACCAGAAGCGTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGaaa                                                                        <  1:526636/143‑1 (MQ=255)
                                                                     cACGTTCTGCCCCGGTTTCACCAGAACTTTCCCCAAATGCATCTAACGTGTGGTGTAGCTGCGCCCATGAGGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCa                                                                     <  2:675066/143‑1 (MQ=255)
                                                                      aCCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACGGCGGTTGGCAACCACCACTTCACCGTCACCGCCTGAACGCAc                                                                    >  1:225788/1‑143 (MQ=255)
                                                                        cgtCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTg                                                                  <  1:196628/141‑1 (MQ=255)
                                                                              tCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACATGCTGCGACACGACGCCTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTAc                                                            >  1:235903/1‑143 (MQ=255)
                                                                                       cACCAGAGTGGTGCGGAAGTGTATAGGACGCGTGGTGTAGCTGCGACCATGAAGAATAGGCACATAATAACCTGTTGCGCGACTGCGTGTGGCGACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCa                                                   <  2:682989/143‑1 (MQ=255)
                                                                                        aCCAGAATCTTGCGCAGGTGCATATAAGGCCTTGTTTAGCTGGGACCATGAAGAATAGCCACATAATAACCTGCTGCGTCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCAt                                                  <  2:212425/143‑1 (MQ=255)
                                                                                        aCCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCAt                                                  >  1:681131/1‑143 (MQ=255)
                                                                                        aCCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCAt                                                  >  1:211029/1‑143 (MQ=255)
                                                                                        aCCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCCGCGGCAt                                                  >  1:380088/1‑143 (MQ=255)
                                                                                          cAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTAACCCACGGAAAGCACTGGCGTACCCTGCGGCATgg                                                >  2:777134/1‑143 (MQ=255)
                                                                                          cAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCAAATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCAACAACTGAAAGCACTGGCGTACCCTGCGGCATgg                                                >  1:107802/1‑143 (MQ=255)
                                                                                            gAATCGTGCGCAAGTGCAGATCACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACCTAATAACCTGCTGCGCCACTGCGTGTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCg                                              <  2:701313/143‑1 (MQ=255)
                                                                                               tCTTGCGCAAGTGCATATAACACGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCCTGGCGaaa                                           >  1:898510/1‑143 (MQ=255)
                                                                                                  ggcgcAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCa                                        <  1:386035/142‑1 (MQ=255)
                                                                                                     gcAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCTCTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAaca                                     <  1:595000/143‑1 (MQ=255)
                                                                                                      cAAGTGCATATAACGCGTGGTGCAGCTGCGACCATGAAGAATAGCCACATACTAACCTGCTGCGCCACTGCGTTTGGCAACCACCACGTCACCGTCACCCCCTGAAAGCACAGCCGTACCCTGCTGCATGGCGAAATCAacac                                    <  1:150300/143‑1 (MQ=255)
                                                                                                          tGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACctct                                <  1:490106/143‑1 (MQ=255)
                                                                                                            cATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACACAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGTGTACCCTGCGGCATGGCGAAACCATCACCTCcgt                              >  2:550066/1‑143 (MQ=255)
                                                                                                                    gcgTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCa                      >  1:849981/1‑143 (MQ=255)
                                                                                                                         gtgtAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGCACCCTGCGGCATGGCCACATCAACACCGCTGCGTGGAGCAGcgcg                 >  2:466981/1‑143 (MQ=255)
                                                                                                                                   gACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGCACCCTGCGCCCCGGCGCAATCAGCACCACCGACTCGTGCAACGCCACCggacccc       >  2:253901/1‑138 (MQ=255)
                                                                                                                                      cATGAAGGATAGCCACATAATACCCTGCTGCGCCACTGCGTTTGGCAACCACCGCTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGa    <  2:235903/143‑1 (MQ=255)
                                                                                                                                       aTGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGAt   >  1:374110/1‑143 (MQ=255)
                                                                                                                                        tGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCGGTGTGGTGCAACGCGCCCGGTCACCGGAtt  >  1:743201/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                           
TTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGATT  >  NC_000913/1940445‑1940723

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: