Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,940,584 C→A 100% 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,5840CA100.0% 44.4 / NA 19R351L (CGT→CTT) mepMpeptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (9/10);  total (9/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGT  >  NC_000913/1940451‑1940714
                                                                                                                                     |                                                                                                                                  
gTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGCATCTTGCGCAAGTGCATATAACGCGCGGTGTAGCTGCGACCATGCAGAATCGCca                                                                                                                           >  2:350905/1‑143 (MQ=255)
    aGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATCCGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAAAATAATAAGCAAGTGCATATAACGCGTGGTGTAGCTGCGACAATGAAGAATAGCCACata                                                                                                                       >  1:86250/1‑143 (MQ=255)
      ctggggCCCGGTTGAACGTCGGGTATTACCGGGAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACataat                                                                                                                     <  2:446124/139‑1 (MQ=255)
       tGCGGCCCGGTTGAACGTCCGGTATTAACGGAAAGCGCGATACGGTCACCACGTTTCACCATCTGTCCCGGTTTGACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGATGCGACAATGAAGAATAGGCACataata                                                                                                                    >  2:462469/1‑143 (MQ=255)
         cGGCCCGGGTGAACGTCCGTTATTACCGGAAAGCGCGATACTGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATTACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAAc                                                                                                                  <  2:144292/143‑1 (MQ=255)
              cGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCTCCTCGTTTCACCTTCTGTGCCGGTTGCACCAGAAACTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACctgct                                                                                                             <  1:295292/143‑1 (MQ=255)
               ggTTGAGCGTCCGGGAGTGGCGGACAGCGCGTGACGGTCGCCTCGTTTCCGCGTGTGTCGCGGTTTCACCAGAATCTTGCGCACGTTCAGATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACctgctg                                                                                                            <  2:109864/143‑1 (MQ=255)
                                 aCCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCaa                                                                                          >  1:363808/1‑143 (MQ=255)
                                          cgcgTTACGATCGCAACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCAGATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACtt                                                                                 <  2:473278/143‑1 (MQ=255)
                                             gATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCAc                                                                              >  2:458087/1‑143 (MQ=255)
                                                      gCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTTCATCTAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCAACCa                                                                     >  1:161846/1‑143 (MQ=255)
                                                        taCGTTTCTCCTTCGGTTCGGGTTTCACCAGAATCTTGCGCATGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACt                                                                   <  2:72998/142‑1 (MQ=255)
                                                                          ccGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAACCACATAATAACCTGCTGCGCCACTGCGTTTGACCACCACCACTACACAATCACACACTGAAAGCACTGGCGTAccc                                                 >  2:283546/1‑143 (MQ=255)
                                                                           cGGTTTCACCAGAATCTTGAACAAGAGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCACCCACCACCTCACCGTCACCCACTGACACCACTGGCGTCCCCt                                                >  2:53571/1‑143 (MQ=255)
                                                                            ggTTTCACCAGAATGTTGCGCAATTTCATATAACGCGTGGTGTAGCTGCGACTATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTg                                               <  1:25465/143‑1 (MQ=255)
                                                                               ttCACCAGAATCTCGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCgg                                            <  1:45076/143‑1 (MQ=255)
                                                                                                         ataACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCtgtgtg                  <  1:458087/143‑1 (MQ=255)
                                                                                                                       aGCTGCGACCATGAAGAATAGCCACCTCATAACCTGCTGCGCCACTGCGTGTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCg    <  2:398518/143‑1 (MQ=255)
                                                                                                                         cTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGCCACCCACTGAAACAAACGCCCTACCCCGCGCCATGGCGAAATCACCACCTCAGCACGATGCACCGCGACCgat  >  2:155957/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGT  >  NC_000913/1940451‑1940714

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