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%	110.7 / NA	28	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/14);  total (14/14)

ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGA  >  NC_000913/1562039‑1562127                                               |                                           aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                          >  1:452994/1‑49 (MQ=255) aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                          <  1:1880507/49‑1 (MQ=255)  ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                         <  1:1094016/49‑1 (MQ=255)    gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                       >  1:775355/1‑49 (MQ=255)    gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                       >  1:248550/1‑49 (MQ=255)      aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                     >  1:1356172/1‑49 (MQ=255)       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                    <  1:2191370/49‑1 (MQ=255)         ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                  >  1:2444248/1‑49 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                            <  1:532244/49‑1 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                            >  1:1655452/1‑49 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                            >  1:2492110/1‑49 (MQ=255)               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                             <  1:1587365/48‑1 (MQ=255)                 gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                          <  1:2432230/49‑1 (MQ=255)                     aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                      >  1:455365/1‑49 (MQ=255)                     aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                      >  1:1481505/1‑49 (MQ=255)                      tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                     <  1:2102760/49‑1 (MQ=255)                      tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                     >  1:805236/1‑49 (MQ=255)                      tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                     <  1:951377/49‑1 (MQ=255)                       cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                    >  1:355934/1‑49 (MQ=255)                         aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                    >  1:2142773/1‑47 (MQ=255)                         aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                    >  1:1369816/1‑47 (MQ=255)                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta               <  1:1868413/47‑1 (MQ=255)                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagt                <  1:1062655/46‑1 (MQ=255)                                aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc            <  1:2395020/48‑1 (MQ=255)                                aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc            <  1:1500412/48‑1 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc        <  1:2067724/46‑1 (MQ=255)                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc        <  1:1088582/46‑1 (MQ=255)                                        ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa    <  1:1672886/48‑1 (MQ=255)                                         gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaag   >  1:1259698/1‑48 (MQ=255)                                           aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga  >  1:813362/1‑47 (MQ=255)                                               |                                           ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGA  >  NC_000913/1562039‑1562127

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