Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,645,776	C→T	G517D (GGC→GAC)	pbpC ←	peptidoglycan glycosyltransferase PbpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,645,776	0	C	T	100.0%	168.0 / NA	47	G517D (GGC→GAC)	pbpC	peptidoglycan glycosyltransferase PbpC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/17);  total (30/17)

GTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATC  >  NC_000913/2645701‑2645852                                                                            |                                                                             gTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAg                                                                         >  1:422448/1‑81 (MQ=255) gTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAg                                                                         >  2:422457/1‑81 (MQ=255)   gCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCt                                                                       >  2:422444/1‑81 (MQ=255)     cGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTgg                                                                     <  1:422401/81‑1 (MQ=255)      gTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGt                                                                    >  2:422442/1‑81 (MQ=255)       tCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTg                                                                   >  1:422440/1‑81 (MQ=255)       tCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTg                                                                   >  1:422445/1‑81 (MQ=255)        ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGc                                                                  >  2:422458/1‑81 (MQ=255)            tCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGt                                                              >  2:422473/1‑81 (MQ=255)              tGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGttt                                                            >  2:422462/1‑81 (MQ=255)               gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGtttt                                                           >  2:422463/1‑81 (MQ=255)                 cAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTcc                                                         <  2:422402/81‑1 (MQ=255)                       aGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGc                                                     >  1:422452/1‑79 (MQ=255)                       aGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCa                                                   <  2:422403/81‑1 (MQ=255)                        gATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAg                                                  >  1:422464/1‑81 (MQ=255)                          tCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTg                                                <  2:422404/81‑1 (MQ=255)                              aaTGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGc                                            <  2:422405/81‑1 (MQ=255)                                  aCATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcg                                        <  2:422406/81‑1 (MQ=255)                                  aCATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcg                                        >  1:422447/1‑81 (MQ=255)                                     tAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcgc                                     <  2:422407/81‑1 (MQ=255)                                         gagcgTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGc                                  >  2:422450/1‑80 (MQ=255)                                         gagcgTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                 <  2:422409/81‑1 (MQ=255)                                         gagcgTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                 <  1:422408/81‑1 (MQ=255)                                          agcgTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGc                                  >  2:422474/1‑79 (MQ=255)                                          agcgTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGc                                  >  1:422454/1‑79 (MQ=255)                                          agcgTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGc                                  >  2:422477/1‑79 (MQ=255)                                             gTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAgg                              >  2:422461/1‑80 (MQ=255)                                             gTTAACCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTCTCCATGCCAGTGGGGCGACGCGCGGCAAgg                              >  2:422480/1‑80 (MQ=255)                                                aaCCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                           <  1:422410/80‑1 (MQ=255)                                                aaCCCCAATCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                           >  2:422476/1‑80 (MQ=255)                                                      aaTCGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCgg                    <  2:422411/81‑1 (MQ=255)                                                         cGCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGc                   >  1:422465/1‑79 (MQ=255)                                                          gCCCAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAg                >  2:422467/1‑81 (MQ=255)                                                           cccAGGCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg               >  1:422460/1‑81 (MQ=255)                                                               ggCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg           >  2:422475/1‑81 (MQ=255)                                                                gCGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgc          <  2:422412/81‑1 (MQ=255)                                                                 cGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg         >  2:422459/1‑81 (MQ=255)                                                                 cGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg         <  1:422413/81‑1 (MQ=255)                                                                 cGTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg         >  1:422441/1‑81 (MQ=255)                                                                  gTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcgc        >  2:422488/1‑81 (MQ=255)                                                                  gTCACGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcgc        >  2:422455/1‑81 (MQ=255)                                                                     aCGATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCtt      >  2:422489/1‑80 (MQ=255)                                                                       gATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcat   <  1:422414/81‑1 (MQ=255)                                                                       gATAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcat   <  1:422415/81‑1 (MQ=255)                                                                        aTAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcatc  <  1:422416/81‑1 (MQ=255)                                                                        aTAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcatc  <  1:422417/81‑1 (MQ=255)                                                                        aTAGTCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcat   >  2:422470/1‑80 (MQ=255)                                                                            |                                                                             GTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATC  >  NC_000913/2645701‑2645852

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑