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%	118.6 / NA	31	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/17);  total (14/17)

ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAA  >  NC_000913/1562039‑1562129                                               |                                             aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                            <  1:990394/49‑1 (MQ=255) aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                            <  1:1251994/49‑1 (MQ=255) aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                           <  1:437638/50‑1 (MQ=255) aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                           >  1:2063582/1‑50 (MQ=255)   cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                         <  1:2700022/50‑1 (MQ=255)    gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                        >  1:2195233/1‑50 (MQ=255)    gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                        <  1:54458/50‑1 (MQ=255)       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                      <  1:35518/49‑1 (MQ=255)       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                     <  1:778802/50‑1 (MQ=255)       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                     <  1:761185/50‑1 (MQ=255)         ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                   <  1:2672148/50‑1 (MQ=255)         ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                    <  1:768505/49‑1 (MQ=255)            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                >  1:2677473/1‑50 (MQ=255)            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                >  1:729779/1‑50 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                              <  1:1404512/49‑1 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                              <  1:1036188/49‑1 (MQ=255)                ggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                            >  1:1549110/1‑50 (MQ=255)                 gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                            >  1:1217808/1‑49 (MQ=255)                 gAGCATCAAGCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                           <  1:1268791/50‑1 (MQ=255)                       cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                      >  1:2371277/1‑49 (MQ=255)                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                      >  1:2579548/1‑48 (MQ=255)                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                    <  1:2052829/50‑1 (MQ=255)                          tCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                   >  1:766445/1‑49 (MQ=255)                            tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                >  1:390487/1‑50 (MQ=255)                             ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc               >  1:2723632/1‑50 (MQ=255)                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                 <  1:496423/47‑1 (MQ=255)                               aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg             >  1:2214408/1‑50 (MQ=255)                               aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg             >  1:1795299/1‑50 (MQ=255)                                aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGc            >  1:106491/1‑50 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa         <  1:182695/47‑1 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa         <  1:1693926/47‑1 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa         <  1:1488788/47‑1 (MQ=255)                                          aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaa  <  1:2074986/50‑1 (MQ=255)                                               |                                             ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAA  >  NC_000913/1562039‑1562129

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