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%	185.2 / NA	47	V72V (GTT→GTC)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/17);  total (30/17)

GTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAA  >  NC_000913/4093178‑4093266                                              |                                            gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                          <  1:1539679/49‑1 (MQ=255) gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                          <  1:1464041/49‑1 (MQ=255)     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                      <  1:2009032/49‑1 (MQ=255)       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                    <  1:1926132/49‑1 (MQ=255)       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                    >  1:306406/1‑49 (MQ=255)        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                   >  1:1631505/1‑49 (MQ=255)          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                 >  1:2274748/1‑49 (MQ=255)          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                 >  1:1002985/1‑49 (MQ=255)          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                 >  1:93099/1‑49 (MQ=255)            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                               <  1:219491/49‑1 (MQ=255)             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                               >  1:919216/1‑48 (MQ=255)               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                            >  1:2034273/1‑49 (MQ=255)                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                          >  1:2276177/1‑49 (MQ=255)                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                          <  1:1595195/49‑1 (MQ=255)                  cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                          >  1:315602/1‑48 (MQ=255)                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                       >  1:269714/1‑49 (MQ=255)                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                       >  1:628227/1‑49 (MQ=255)                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                      >  1:1189609/1‑49 (MQ=255)                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                      >  1:919856/1‑49 (MQ=255)                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                     >  1:2028741/1‑49 (MQ=255)                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                    <  1:1104115/49‑1 (MQ=255)                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                     >  1:2158203/1‑48 (MQ=255)                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                   <  1:382693/49‑1 (MQ=255)                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                  >  1:1485305/1‑49 (MQ=255)                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                  >  1:1415617/1‑49 (MQ=255)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                  <  1:1360339/48‑1 (MQ=255)                             ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt              >  1:472389/1‑49 (MQ=255)                             ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt              >  1:1494746/1‑49 (MQ=255)                             ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt              >  1:1205155/1‑49 (MQ=255)                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt              <  1:208837/48‑1 (MQ=255)                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt              <  1:289553/48‑1 (MQ=255)                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct            >  1:237126/1‑49 (MQ=255)                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc             <  1:233060/48‑1 (MQ=255)                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc          <  1:1170245/49‑1 (MQ=255)                                 accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt           >  1:1628389/1‑48 (MQ=255)                                  ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt           >  1:1580694/1‑47 (MQ=255)                                  ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt           >  1:202126/1‑47 (MQ=255)                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt           >  1:2334942/1‑45 (MQ=255)                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt           >  1:1317126/1‑45 (MQ=255)                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt           >  1:1156021/1‑45 (MQ=255)                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa       <  1:293328/49‑1 (MQ=255)                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt      >  1:2272943/1‑49 (MQ=255)                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt     <  1:1376300/49‑1 (MQ=255)                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt     <  1:2189811/49‑1 (MQ=255)                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt     >  1:2318023/1‑48 (MQ=255)                                        ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa   <  1:2006948/49‑1 (MQ=255)                                         tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa  >  1:1103497/1‑49 (MQ=255)                                              |                                            GTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAA  >  NC_000913/4093178‑4093266

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