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% 53.7 / NA 20D104N (GAT→AAT) mepMpeptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/7);  total (13/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGC  >  NC_000913/1941188‑1941440
                                                                                                                                          |                                                                                                                  
gCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATAGAGAATGCTGCTTAGCGTATTGCCa                                                                                                                >  2:44708/1‑143 (MQ=255)
   gTTCGCCGTCCGCGGTTAATGTACAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACACATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCCGTg                                                                                                             >  2:175319/1‑143 (MQ=255)
        ccGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCAATATCAATAACAAATGGAGTGAGAATGCTGCTTAGCGTATTGCCAGTGGAaac                                                                                                        >  1:483471/1‑143 (MQ=255)
             cgcgGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGGGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCGTAGCGTATTGGCAGTGGAAAGAACat                                                                                                   >  2:402893/1‑143 (MQ=255)
               cgGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCCTATCAATACCGTACTGATTGAGAATGCGGCTTAGCGTATTGCCAGTGGAAACTACatag                                                                                                 >  2:26602/1‑142 (MQ=255)
                     aTGTCCAGGCGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGTATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGc                                                                                           <  2:17821/143‑1 (MQ=255)
                                                    gTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACGATACTGATCGAGAGTGCTGCTTAGCGTATTGCCAGTGGAAACACCATATTCATGCACGCCCGCTTCACCGGCGACTTTGTCATCCa                                                            >  2:131770/1‑143 (MQ=255)
                                                          cAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTTAGAAGGCTGCGTAGCGTATTGCCAGTGGCAACAACATATTCATGCACGCCCGCTTCACCTGCGCCTTTGTAATCGAGTTctt                                                      >  2:109354/1‑141 (MQ=255)
                                                               cTTCGTCGGCCGCAGCCAGTTGGGTGATAACAATCATATCTGTAGCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTgg                                                 <  1:329841/143‑1 (MQ=255)
                                                                      ggCCGCAGCCAGTTGGGTGATATCACCTATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCGCGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATgg                                          <  2:41711/143‑1 (MQ=255)
                                                                           cAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTCGTCATCCAGTTCGTCCTGGGGAATGGCTTcc                                     >  2:463012/1‑142 (MQ=255)
                                                                            aGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCTCGTTCGTCCGGGGGAATGGCTCCAt                                    >  2:484885/1‑143 (MQ=255)
                                                                                      gTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACCCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGc                          >  1:244186/1‑143 (MQ=255)
                                                                                       tGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCa                         <  1:309429/143‑1 (MQ=255)
                                                                                            tCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTTTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGTAGCTTg                    <  1:126809/143‑1 (MQ=255)
                                                                                                ccATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCa                >  2:413900/1‑143 (MQ=255)
                                                                                                     tCAATACCATACGGATTGAGAATGCTGCTGAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCgg           <  1:154857/143‑1 (MQ=255)
                                                                                                           ccATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAACCAACAGATTCATGCACGCCCGCTGCACCGGCGATTTTGTCATCCAGTTCAGCCTGGGGAATGGCTGCATCTTCTTGTGCAGCTTAATCAAGCGGCTCACt     >  2:271432/1‑143 (MQ=255)
                                                                                                             aTACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCACTTTACCGGCGAGTTTGTCATACAGTTCGTCCTGGGGACTGGCTTCATCTTCTTGTGTAGCTTGATACACCGGCTGAATgg   >  2:7153/1‑143 (MQ=255)
                                                                                                              tACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGc  <  2:286883/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                  
GCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGC  >  NC_000913/1941188‑1941440

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: