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

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269                                                |                                               gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:1193719/50‑1 (MQ=255) gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:2826892/50‑1 (MQ=255)  ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:1304076/1‑50 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            <  1:2406140/50‑1 (MQ=255)    ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:2592347/1‑49 (MQ=255)    ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:851232/50‑1 (MQ=255)       gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         <  1:5026791/49‑1 (MQ=255)       gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:2401263/50‑1 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:3681269/1‑50 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:2547988/1‑50 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:2116911/1‑50 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     >  1:3215594/1‑50 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:3710624/50‑1 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:4912574/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:177852/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:2046355/1‑50 (MQ=255)            cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:3030818/1‑50 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:4433479/1‑49 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:2226888/1‑49 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:1768719/50‑1 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:1245731/50‑1 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:4342964/50‑1 (MQ=255)              tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  >  1:4798727/1‑49 (MQ=255)                aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:2575645/50‑1 (MQ=255)                aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               >  1:5075456/1‑50 (MQ=255)                 cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              >  1:4950284/1‑50 (MQ=255)                   tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:1716736/1‑49 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:488917/50‑1 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:1565611/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:2991047/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:1080523/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:1342521/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:340095/1‑50 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:892184/1‑50 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:3641778/50‑1 (MQ=255)                       tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:630556/50‑1 (MQ=255)                        cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:4761516/1‑49 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:4452841/49‑1 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:1824417/1‑48 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:1537508/1‑48 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:1473255/50‑1 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:553238/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:2889278/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     >  1:588881/1‑50 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:4845810/50‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCTACCGCATTCCAGCg                   <  1:3152876/50‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:3285674/49‑1 (MQ=255)                             atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:4402663/1‑50 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCGGCCGCATTCCAGCGtt                 <  1:4686964/50‑1 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 <  1:4907641/50‑1 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:5010740/1‑50 (MQ=255)                               ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                >  1:1726745/1‑50 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:4193584/49‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:2193424/50‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:16428/1‑49 (MQ=255)                                   accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4895455/1‑48 (MQ=255)                                   accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4932225/1‑48 (MQ=255)                                   accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:349593/50‑1 (MQ=255)                                   accaACGTTGGCAGACATCACTGCGGGCAACCGCATTCCAGCGTTCTTCa            <  1:511523/50‑1 (MQ=255)                                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:3488995/1‑47 (MQ=255)                                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg           <  1:4096285/50‑1 (MQ=255)                                     caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:1364138/1‑46 (MQ=255)                                     caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:730975/50‑1 (MQ=255)                                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:1524485/1‑50 (MQ=255)                                      aaCGTCGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:3719909/1‑50 (MQ=39)                                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:4533236/1‑50 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:4460954/1‑49 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:3224394/50‑1 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:4219279/1‑49 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:2645871/50‑1 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:204960/50‑1 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:4476904/1‑50 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:3279309/1‑50 (MQ=255)                                           tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:3736302/1‑50 (MQ=255)                                           tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:705154/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   >  1:3235245/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:232572/49‑1 (MQ=255)                                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:3616794/50‑1 (MQ=255)                                                |                                               GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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