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

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093176‑4093268                                                |                                              gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                             <  1:7057651/50‑1 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           <  1:1456069/50‑1 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           <  1:4768784/50‑1 (MQ=255)    ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                          <  1:442186/50‑1 (MQ=255)     tAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           >  1:6466362/1‑48 (MQ=255)      aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                        >  1:4352388/1‑50 (MQ=255)      aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                        >  1:2679962/1‑50 (MQ=255)       gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                        <  1:5858228/49‑1 (MQ=255)       gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                       <  1:4117462/50‑1 (MQ=255)        cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                      >  1:2441500/1‑50 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     <  1:2201159/50‑1 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     >  1:14990/1‑50 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                    <  1:3411138/50‑1 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                    >  1:1887202/1‑50 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   >  1:585703/1‑50 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   >  1:827089/1‑50 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:1816572/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:1212119/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:7066133/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   >  1:2555962/1‑50 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:2077820/50‑1 (MQ=255)            cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  <  1:5115453/50‑1 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  >  1:2523361/1‑49 (MQ=255)              tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:5785503/1‑49 (MQ=255)              tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:5877248/1‑49 (MQ=255)              tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:4904995/1‑49 (MQ=255)              tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:3162124/1‑49 (MQ=255)               aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                               >  1:5870426/1‑50 (MQ=255)               aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                               >  1:6329946/1‑50 (MQ=255)                aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                              <  1:5523864/50‑1 (MQ=255)                 cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCGAcc                             <  1:1778026/50‑1 (MQ=255)                 cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                             >  1:6654590/1‑50 (MQ=255)                  aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCATCCg                            <  1:5701956/50‑1 (MQ=255)                  aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:4165454/1‑50 (MQ=255)                   tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGc                           <  1:3298345/50‑1 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         >  1:3128729/1‑50 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        <  1:4692059/50‑1 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:2851995/1‑50 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:3945608/1‑50 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:4920505/1‑50 (MQ=255)                      gTCATATATTTCCACCAACGCTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:5343054/1‑50 (MQ=255)                       tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       <  1:3907753/50‑1 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                      <  1:791686/49‑1 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       >  1:6371306/1‑48 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                     <  1:5487751/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                    <  1:3809861/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                    <  1:2487644/50‑1 (MQ=255)                           atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   >  1:2519833/1‑50 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   <  1:3183663/49‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                  >  1:6801965/1‑50 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                >  1:5753479/1‑50 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                >  1:2773454/1‑50 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                >  1:6774986/1‑50 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                 <  1:5546579/49‑1 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                 <  1:4875299/49‑1 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct              >  1:5777193/1‑50 (MQ=255)                                 ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct              <  1:475615/49‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCTACCGCATTCCAGCGttctt             <  1:4974244/49‑1 (MQ=39)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc            <  1:5738678/50‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc            <  1:5535655/50‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             <  1:7059988/49‑1 (MQ=255)                                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:5607919/1‑47 (MQ=255)                                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg          <  1:5915022/50‑1 (MQ=255)                                     caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa         <  1:5240821/50‑1 (MQ=255)                                     caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa         <  1:6557516/50‑1 (MQ=255)                                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt        <  1:818787/50‑1 (MQ=255)                                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt        >  1:2963970/1‑50 (MQ=255)                                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       <  1:5423927/50‑1 (MQ=255)                                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       <  1:5925913/50‑1 (MQ=255)                                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       <  1:4449541/50‑1 (MQ=255)                                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       <  1:5772110/50‑1 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       >  1:2642339/1‑49 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      <  1:284384/50‑1 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:6892163/1‑49 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:4295797/1‑49 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:362872/1‑49 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa    >  1:6764631/1‑50 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa    >  1:4154352/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc  >  1:2069429/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt   <  1:839928/49‑1 (MQ=255)                                                |                                              GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093176‑4093268

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