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

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269                                                |                                               gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:4086175/50‑1 (MQ=255) gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:6266651/50‑1 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:2875682/1‑50 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:11298/1‑50 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            <  1:3319717/50‑1 (MQ=255)   gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            <  1:6166828/50‑1 (MQ=255)    ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:4436940/1‑49 (MQ=255)    ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:4589239/1‑49 (MQ=255)    ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:2283977/50‑1 (MQ=255)      aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:899548/1‑50 (MQ=255)      aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:1799553/1‑50 (MQ=255)       gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        >  1:2895800/1‑50 (MQ=255)        cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:3823150/1‑50 (MQ=255)        cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:6444811/50‑1 (MQ=255)        cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:6349179/1‑50 (MQ=255)         gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:110365/1‑50 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:2492117/50‑1 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:4446631/50‑1 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:6479547/50‑1 (MQ=255)          ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:3960516/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:6254832/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:6558343/50‑1 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:1331162/1‑50 (MQ=255)           gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:3909394/50‑1 (MQ=255)            cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:2732528/1‑50 (MQ=255)            cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   <  1:5804879/50‑1 (MQ=255)            cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   <  1:21142/50‑1 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:1813194/1‑49 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:1991430/1‑49 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:6536269/50‑1 (MQ=255)             aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:5056904/50‑1 (MQ=255)              tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  >  1:3423370/1‑49 (MQ=255)               aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:707099/1‑50 (MQ=255)               aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:747314/1‑50 (MQ=255)               aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:1146697/1‑50 (MQ=255)                aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:2118388/50‑1 (MQ=255)                aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:2405154/50‑1 (MQ=255)                aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:1549418/50‑1 (MQ=255)                 cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              >  1:5337476/1‑50 (MQ=255)                 cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:4802326/50‑1 (MQ=255)                  aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:4447081/1‑50 (MQ=255)                   tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:1036362/1‑49 (MQ=255)                    cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           >  1:1858227/1‑50 (MQ=255)                    cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           >  1:4333731/1‑50 (MQ=255)                    cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           >  1:6114632/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:2874861/1‑50 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:2345289/50‑1 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:691311/50‑1 (MQ=255)                     ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:1069482/50‑1 (MQ=255)                      gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:5788170/1‑50 (MQ=255)                       tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:5973066/50‑1 (MQ=255)                       tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:6388634/1‑50 (MQ=255)                       tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5481052/1‑50 (MQ=255)                       tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:4537035/1‑50 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCTACCGCATTCCa                      <  1:6617407/50‑1 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:4256093/1‑48 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:6642959/1‑48 (MQ=255)                         atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:1795099/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:5477027/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:938253/50‑1 (MQ=255)                          tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:107465/50‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAGCCGCATTCCAGc                    <  1:4844949/49‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:3980406/49‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:4975664/1‑50 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:4696118/1‑50 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:2221554/50‑1 (MQ=255)                            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:2432385/1‑50 (MQ=255)                              tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:6001217/1‑50 (MQ=255)                               ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                >  1:2235651/1‑50 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               >  1:114769/1‑50 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               >  1:4096003/1‑50 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               >  1:3640797/1‑50 (MQ=255)                                tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:4889443/49‑1 (MQ=255)                                 ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               <  1:4416550/49‑1 (MQ=255)                                  caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:80117/50‑1 (MQ=255)                                   accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:74012/50‑1 (MQ=255)                                     caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:6587019/1‑46 (MQ=255)                                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:6387414/50‑1 (MQ=255)                                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:2899315/1‑50 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:1714862/1‑49 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:3348937/50‑1 (MQ=255)                                        cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:2222623/50‑1 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:4694004/1‑49 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:3762184/1‑49 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:5733653/1‑49 (MQ=255)                                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:6259217/1‑49 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:4408210/50‑1 (MQ=255)                                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:1855426/1‑50 (MQ=255)                                           tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:1852460/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   >  1:5858492/1‑50 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:2665087/49‑1 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:3546311/49‑1 (MQ=255)                                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:3944628/49‑1 (MQ=255)                                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  >  1:1637343/1‑50 (MQ=255)                                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:3019840/49‑1 (MQ=255)                                                |                                               GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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