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

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269                                               |                                               ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:1308650/1‑50 (MQ=255) ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:6885210/1‑50 (MQ=255)   ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:6885359/50‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:6045878/49‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:3000565/49‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:6764980/49‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:6307366/49‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:1383727/1‑50 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:5606646/50‑1 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:6884606/50‑1 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:6549171/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:4516781/50‑1 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:1532402/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:7419574/1‑50 (MQ=255)         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:6632211/50‑1 (MQ=255)         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     >  1:2870546/1‑50 (MQ=255)          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:2795225/1‑50 (MQ=255)           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   <  1:3994742/50‑1 (MQ=255)           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   <  1:918281/50‑1 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:6305073/50‑1 (MQ=255)             tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  >  1:1223395/1‑49 (MQ=255)              aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:4624696/1‑50 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:1846942/50‑1 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:6442415/50‑1 (MQ=255)                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:8102828/1‑50 (MQ=255)                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             <  1:1871757/50‑1 (MQ=255)                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:2146646/1‑50 (MQ=255)                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:325014/1‑50 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:5286948/1‑50 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:4277078/50‑1 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:2299520/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:1604179/1‑50 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:3915747/50‑1 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:7325266/1‑49 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5665746/1‑48 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:3085237/1‑48 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:6180746/1‑48 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:7614487/1‑48 (MQ=255)                         tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:1538830/50‑1 (MQ=255)                          atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:4115245/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:525693/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:6875437/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:7349994/50‑1 (MQ=255)                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:7736664/1‑50 (MQ=255)                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:2039982/1‑50 (MQ=255)                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:6405876/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:583216/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:7334927/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  <  1:7173175/49‑1 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  <  1:3263633/49‑1 (MQ=255)                              ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                >  1:4954747/1‑50 (MQ=255)                               tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:2770423/49‑1 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4848055/1‑50 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:317700/1‑50 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:3225226/1‑50 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:2802151/1‑49 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:2503592/1‑48 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:3138883/50‑1 (MQ=255)                                   ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg           <  1:7800393/50‑1 (MQ=255)                                   ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg           <  1:7534289/50‑1 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:1548317/1‑46 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:8180024/1‑46 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:3984896/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:1975988/1‑50 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:4959214/1‑50 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:1806678/1‑50 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:1356720/1‑50 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:3361933/1‑50 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:8175736/50‑1 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:2019190/1‑49 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:2278848/1‑49 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:496519/1‑49 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:1777408/50‑1 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:5762338/50‑1 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:2465072/50‑1 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:1034533/1‑49 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:341557/1‑49 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:4029321/1‑49 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:7268421/1‑49 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:5185216/50‑1 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:2358647/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:1750238/50‑1 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:5039321/1‑50 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:6483886/1‑50 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:5934463/50‑1 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:3581843/50‑1 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  >  1:6486177/1‑50 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:5227535/49‑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: ‑