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%	106.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 (18/13);  total (18/13)

CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562038‑1562131                                                |                                               cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              >  1:100685/1‑50 (MQ=255) cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              >  1:2612104/1‑50 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                              <  1:1382437/49‑1 (MQ=255)  aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                             >  1:1769911/1‑50 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            >  1:397371/1‑50 (MQ=255)   ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                            <  1:1748190/50‑1 (MQ=255)    cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                           <  1:1930415/50‑1 (MQ=255)     gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                          >  1:107253/1‑50 (MQ=255)          ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                     <  1:2234957/50‑1 (MQ=255)               cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                >  1:1005327/1‑50 (MQ=255)                gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                               <  1:1148832/50‑1 (MQ=255)                 ggAGCATCAATCTTTAAATCTTTAGGAAATGTTTCTGCTGGCCCCTGGcc                              <  1:2094093/50‑1 (MQ=255)                  gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                             <  1:2532721/50‑1 (MQ=255)                   aGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGa                            >  1:1444316/1‑50 (MQ=255)                     cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                          >  1:2128068/1‑50 (MQ=255)                        cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:336348/1‑49 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttttt                       <  1:476642/49‑1 (MQ=255)                         aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                        >  1:870904/1‑48 (MQ=255)                          aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                     >  1:1614897/1‑50 (MQ=255)                          aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                     >  1:1540776/1‑50 (MQ=255)                             tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                  >  1:2474979/1‑50 (MQ=255)                              ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                 >  1:1301167/1‑50 (MQ=255)                              ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                 >  1:30495/1‑50 (MQ=255)                                aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                <  1:43232/49‑1 (MQ=255)                                aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg               >  1:2265993/1‑50 (MQ=255)                                 aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGc              <  1:2180587/50‑1 (MQ=255)                                  aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt             <  1:373732/50‑1 (MQ=255)                                  aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt             >  1:1254706/1‑50 (MQ=255)                                      ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:893583/1‑49 (MQ=255)                                        gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa          >  1:672178/1‑47 (MQ=255)                                         ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      <  1:2386862/50‑1 (MQ=255)                                         ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga      <  1:1974456/50‑1 (MQ=255)                                          gAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaa     <  1:1554634/50‑1 (MQ=255)                                             aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt  >  1:64793/1‑50 (MQ=255)                                                |                                               CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGT  >  NC_000913/1562038‑1562131

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