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

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269                                               |                                               ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             <  1:5506516/50‑1 (MQ=255) ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:790459/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:2989538/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:4695352/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:3386847/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:3932861/1‑50 (MQ=255)    tAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:5955742/1‑48 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:2566078/49‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:1098509/1‑50 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:1835275/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:2087621/1‑50 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:2092307/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:1305441/1‑50 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:2832224/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:342510/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:4118257/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:3884511/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:2730185/1‑50 (MQ=255)        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      <  1:2729333/50‑1 (MQ=255)         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:5106365/50‑1 (MQ=255)          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:3916187/1‑50 (MQ=255)           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:5529548/1‑50 (MQ=255)           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:878231/1‑50 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:5043116/1‑49 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:5716132/50‑1 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:843602/50‑1 (MQ=255)             tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                 <  1:1903870/50‑1 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               >  1:2970646/1‑50 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:1567156/50‑1 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:1622639/50‑1 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:1826315/50‑1 (MQ=255)                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:3775118/1‑50 (MQ=255)                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGCGCTGGCAACCGCAt                          >  1:6205991/1‑50 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:5664728/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:2792130/1‑50 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:2314546/1‑49 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5139571/1‑49 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:1279152/49‑1 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:6341279/1‑48 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:4711149/1‑48 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5019842/1‑48 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:5039538/1‑48 (MQ=255)                         tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     >  1:5663126/1‑50 (MQ=255)                          atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:2109814/1‑50 (MQ=255)                          atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:5449370/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:1416385/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:3282896/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:3818158/50‑1 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 <  1:6018360/50‑1 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:5947412/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  <  1:3705381/49‑1 (MQ=255)                              ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                >  1:4545825/1‑50 (MQ=255)                               tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:1940568/49‑1 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4967834/1‑50 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               <  1:5565271/49‑1 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:4761470/49‑1 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:97084/1‑48 (MQ=255)                                  accaACGTTGGCAGACAACAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:459394/50‑1 (MQ=39)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:890727/1‑46 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:4148730/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:2041145/50‑1 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:412015/1‑50 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:5665165/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:3413920/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:652267/1‑50 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:2997353/1‑50 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:979702/1‑50 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:5486002/1‑50 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:3928708/1‑50 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   >  1:5099394/1‑50 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   >  1:3426484/1‑50 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:5184252/49‑1 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:3234881/50‑1 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  >  1:1273103/1‑50 (MQ=255)                                            gCAGACATCAGAGCTGGCAACCGCATTCTAGCGTTCTTCAGATTCAAtct  >  1:5045503/1‑50 (MQ=255)                                               |                                               GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑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: ‑