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

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269                                               |                                               ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:4105438/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:1779810/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:1711827/1‑50 (MQ=255)   ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:279919/50‑1 (MQ=255)   ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:1875989/50‑1 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         <  1:3653933/49‑1 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:789585/50‑1 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:795316/50‑1 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:5362105/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:3328677/1‑50 (MQ=255)         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:2251050/50‑1 (MQ=255)          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:3557939/1‑50 (MQ=255)           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:1778458/1‑50 (MQ=255)             tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                 <  1:3934606/50‑1 (MQ=255)              aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:4222442/1‑50 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:4293928/50‑1 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:126400/50‑1 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              >  1:1019970/1‑50 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:1558662/50‑1 (MQ=255)                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCTACCg                             <  1:4356728/50‑1 (MQ=255)                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:5369786/1‑50 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:5136539/1‑50 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:2573289/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:5549495/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:1516288/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:3811440/1‑50 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5339701/1‑49 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:2570145/49‑1 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:419040/50‑1 (MQ=255)                         tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     >  1:2040820/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:3471042/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:2306228/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:4502804/49‑1 (MQ=255)                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:2686249/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 <  1:1133033/50‑1 (MQ=255)                              ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                >  1:3672301/1‑50 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:2619765/1‑50 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               <  1:1779664/49‑1 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:5409667/50‑1 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4627060/1‑48 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4784218/1‑48 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:101356/50‑1 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:3271645/50‑1 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:5434939/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCTACCGCATTCCAGCGTTCTTCAGAt         <  1:4041064/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:1580119/50‑1 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:4331171/50‑1 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:1176375/1‑50 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:3963759/1‑49 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:987638/50‑1 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:67683/1‑49 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:9281/1‑49 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:2394849/1‑50 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:1303690/49‑1 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:496977/50‑1 (MQ=255)                                               |                                               GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269

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