Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,941,326 C→T 100% D104N (GAT→AAT)  mepM ← peptidoglycan DD‑endopeptidase MepM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,941,3260CT100.0% 69.5 / NA 26D104N (GAT→AAT) mepMpeptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/11);  total (15/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATTTC  >  NC_000913/1941192‑1941464
                                                                                                                                      |                                                                                                                                          
ttCGCCGTCCGCGGATAATGACCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCAGCCGCAGCCAGTTGGGTGATATCACCCATAGCAATACCATACTGATTGCGAATGCTGCTTCGCGTATTGCCAGtcg                                                                                                                                    >  2:192488/1‑141 (MQ=255)
    ggggccGCGGGTAATTTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCACTACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAaac                                                                                                                                <  1:709484/139‑1 (MQ=255)
           cgGTTAATGTCCAGGAGAGTTGTTGACCGATCTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTTGGGGAGATCACCCATATCAATACCATACTCAATGACAATGACACTTAGCGTATTACCAGTGGAaacaaactac                                                                                                                         >  2:796001/1‑138 (MQ=255)
             gTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTc                                                                                                                       <  1:196267/143‑1 (MQ=255)
                                    aCCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGa                                                                                                >  2:17804/1‑143 (MQ=255)
                                     ccGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGAt                                                                                               >  1:527591/1‑143 (MQ=255)
                                      cGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCCCCGGCGAtt                                                                                              >  1:129070/1‑143 (MQ=255)
                                      cGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGACTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCAgcggcgggt                                                                                              >  2:230457/1‑140 (MQ=255)
                                         ttttCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCCTCACCGGCGATTTTg                                                                                           >  1:623556/1‑143 (MQ=255)
                                          tctCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCACTACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGt                                                                                          <  2:784435/141‑1 (MQ=255)
                                                 ttACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCACTACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATGTTGTCATcccg                                                                                   >  1:459596/1‑141 (MQ=255)
                                                   aCGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATGCAGtt                                                                                 >  1:465661/1‑143 (MQ=255)
                                                   aCGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATACCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGGCATTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGtt                                                                                 <  1:758488/143‑1 (MQ=255)
                                                      cAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATCCTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGt                                                                              >  2:353927/1‑143 (MQ=255)
                                                            tttGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAACCAACATATTCTTGCACGCCCGCTTCACCGGCGATTTTGTAATCCAGTATGTCATggg                                                                        >  1:31551/1‑143 (MQ=255)
                                                                              gtgtgTGATATCACCCATCTCAATCCTCTACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATcttctt                                                      <  2:662493/139‑1 (MQ=255)
                                                                                     atatCACCCATTGCGATTCCAGACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGc                                               <  2:200256/143‑1 (MQ=255)
                                                                                          actcATATCAATTCCAGCCTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGAt                                          <  2:130674/140‑1 (MQ=255)
                                                                                              atatCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCCTCTTCTTGTGCAGCTGGCTCAAt                                      >  2:717511/1‑143 (MQ=255)
                                                                                                    aTACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTc                                >  1:378228/1‑143 (MQ=255)
                                                                                                     tACCATACTGAGTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCa                               <  2:137740/143‑1 (MQ=255)
                                                                                                             tGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTc                       >  1:558916/1‑143 (MQ=255)
                                                                                                                 tGAGAATGCTGCTTAGCGTATTGCCAGTGGGAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGt                   <  2:93256/143‑1 (MQ=255)
                                                                                                                       tgctgcGTAGCGTATTGCCAGTGGAAACAACATATTCGTGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGc             <  2:159741/143‑1 (MQ=255)
                                                                                                                          tggttAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGCTTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGcgaa          <  2:185216/140‑1 (MQ=255)
                                                                                                                                  gTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTCCGTCGTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTATGAGCGAACGAATTTc  >  1:757344/1‑143 (MQ=255)
                                                                                                                                      |                                                                                                                                          
TTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATTTC  >  NC_000913/1941192‑1941464

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