Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,093,223	A→G	V72V (GTT→GTC)	rhaM ←	L‑rhamnose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,223	0	A	G	100.0%	131.1 / NA	36	V72V (GTT→GTC)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/18);  total (18/18)

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269                                                |                                               gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:1713264/50‑1 (MQ=255)  ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:1985569/1‑50 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:164596/1‑50 (MQ=255)     tAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:2672058/1‑48 (MQ=255)       gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        >  1:1203867/1‑50 (MQ=255)        cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:1114169/1‑50 (MQ=255)        cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:42776/1‑50 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:1094958/1‑50 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:248123/1‑50 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:1098274/50‑1 (MQ=255)               aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:1377502/1‑50 (MQ=255)                 cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:1912666/50‑1 (MQ=255)                  aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             <  1:2275864/50‑1 (MQ=255)                    cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           >  1:688896/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           <  1:1493524/49‑1 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           <  1:1614848/49‑1 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:104562/50‑1 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:602123/1‑48 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:1343827/49‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:115412/49‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:119577/50‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:1155322/1‑50 (MQ=255)                             atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:1705411/1‑50 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:898639/49‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:1139665/50‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:92682/49‑1 (MQ=255)                                   accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:118774/50‑1 (MQ=255)                                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:2695338/1‑47 (MQ=255)                                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:507269/1‑50 (MQ=255)                                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:2322227/1‑50 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:2628442/50‑1 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:1089472/50‑1 (MQ=255)                                           tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:2536155/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   >  1:863308/1‑50 (MQ=255)                                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:2640859/50‑1 (MQ=255)                                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:205892/49‑1 (MQ=255)                                                |                                               GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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: ‑