Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,562,085	C→T	A133T (GCA→ACA)	ddpA ←	putative D,D‑dipeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,562,085	0	C	T	100.0%	112.6 / NA	29	A133T (GCA→ACA)	ddpA	putative D,D‑dipeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/15);  total (14/15)

ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562039‑1562131                                               |                                               aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:1942605/49‑1 (MQ=255) aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             <  1:2192104/50‑1 (MQ=255)  ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            <  1:1138124/50‑1 (MQ=255)   cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                           <  1:1982523/50‑1 (MQ=255)     tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                         <  1:914476/50‑1 (MQ=255)     tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                         <  1:491346/50‑1 (MQ=255)     tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                         <  1:1857702/50‑1 (MQ=255)       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                       >  1:2002397/1‑50 (MQ=255)        gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                      >  1:117467/1‑50 (MQ=255)         ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                      <  1:580828/49‑1 (MQ=255)            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                   <  1:360754/49‑1 (MQ=255)            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                   <  1:2186438/49‑1 (MQ=255)              cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                >  1:741352/1‑50 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                               <  1:1702606/50‑1 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                               >  1:715134/1‑50 (MQ=255)                   gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                           <  1:1878077/50‑1 (MQ=255)                     aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                         >  1:972408/1‑50 (MQ=255)                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:471580/1‑48 (MQ=255)                         aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                     >  1:1786495/1‑50 (MQ=255)                                 aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt             <  1:791001/50‑1 (MQ=255)                                   cTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc            >  1:1669111/1‑49 (MQ=255)                                    tttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa           >  1:1703371/1‑49 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc            <  1:272923/46‑1 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa           <  1:193455/47‑1 (MQ=255)                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:1110080/1‑47 (MQ=255)                                        ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaag       <  1:1240053/49‑1 (MQ=255)                                          aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:1723282/1‑48 (MQ=255)                                          aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      >  1:1700538/1‑48 (MQ=255)                                           aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag   >  1:819696/1‑50 (MQ=255)                                            aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt  >  1:1477270/1‑50 (MQ=255)                                               |                                               ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562039‑1562131

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