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

CATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTT  >  NC_000913/4093192‑4093254                                |                                cATCGGTCATATATTTCCCCCAACGTTGGCAGACa                              <  1:679832/35‑1 (MQ=25) cATCGGTCATATATTTCCACCAACGTTGGCAGACa                              <  1:6332946/35‑1 (MQ=38) cATCGGTCATATATTTCCACCAACGTTGGCAGACa                              <  1:4191350/35‑1 (MQ=38)  aTCGGTCATATATTTCCACCAACGTTGGCAGACAt                             >  1:5103280/1‑35 (MQ=255)  aTCGGTCATATATTTCCACCAACGTTGGCAGACAt                             >  1:881202/1‑35 (MQ=38)  aTCGGTCATATATTTCCACCAACGTTGGCAGACAt                             <  1:4719092/35‑1 (MQ=38)  aTCGGTCATATATTTCCACCAACGTTGGCAGACAt                             <  1:4719056/35‑1 (MQ=38)  aTCGGTCATATATTTCCACCAACGTTGGCAGACAt                             <  1:5490414/35‑1 (MQ=255)  aTCGGTCATATATTTCCACCAACGTTGGCAGACAt                             <  1:2746800/35‑1 (MQ=255)   tCGGTCATATATTTCCACCAACGTTGGCAGACATc                            >  1:1935452/1‑35 (MQ=38)   tCGGTCATATATTTCCACCAACGTTGGCAGACATc                            <  1:8234/35‑1 (MQ=255)   tCGGTCATATATTTCCACCAACGTTGGCAGACATc                            >  1:65784/1‑35 (MQ=38)     ggTCATATATTTCCACCAACGTTGGCAGACATCac                          >  1:4619654/1‑34 (MQ=38)      gTCATATATTTCCACCAACGTTGGCAGACATCAGt                         >  1:36161/1‑35 (MQ=255)      gTCATATATTTCCACCAACGTTGGCAGACATCAGt                         >  1:2826302/1‑35 (MQ=255)      gTCATATATTTCCACCAACGTTGGCAGACATCAGt                         >  1:2574603/1‑35 (MQ=255)      gTCATATATTTCCACCAACGTTGGCAGACATCAGt                         >  1:6280151/1‑35 (MQ=38)        cATATATTTCCACCAACGTTGGCAGACATCAGTGc                       >  1:4041420/1‑35 (MQ=38)         atatatTTCCCCCAACGTTGGCAGACATCAGTGCt                      <  1:5441926/35‑1 (MQ=25)         atatatTTCCACCAACGTTGGCAGACATCAGTGCt                      >  1:5260518/1‑35 (MQ=38)         atatatTTCCACCAACGTTGGCAGACATCAGTGCt                      <  1:7744880/35‑1 (MQ=38)          tatatTTCCACCAACGTTGGCAGACATCAGTGCTg                     <  1:2480289/35‑1 (MQ=38)          tatatTTCCACCAACGTTGGCAGACATCAGTGCTg                     >  1:867459/1‑35 (MQ=255)          tatatTTCCACCAACGTTGGCAGACATCAGTGCTg                     <  1:2159558/35‑1 (MQ=38)           atatTTCCACCAACGTTGGCAGACATCAGTGCTgg                    <  1:7125707/35‑1 (MQ=255)           atatTTCCACCAACGTTGGCAGACATCAGTGCTgg                    >  1:812266/1‑35 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   >  1:1348888/1‑35 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   >  1:5658222/1‑35 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   >  1:5719808/1‑35 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   >  1:6220628/1‑35 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   <  1:6534822/35‑1 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   >  1:6814470/1‑35 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   <  1:2349972/35‑1 (MQ=255)            tatTTCCACCAACGTTGGCAGACATCAGTGCTGGc                   >  1:2308109/1‑35 (MQ=38)             atTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                  >  1:55669/1‑35 (MQ=38)              tttCCACCAACGTTGGCAGACATCAGTGCTGGCaa                 >  1:1905010/1‑35 (MQ=255)              tttCCACCAACGTTGGCAGACATCAGTGCTGGCaa                 <  1:6121585/35‑1 (MQ=255)               ttCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                <  1:6626908/35‑1 (MQ=255)               ttCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                <  1:6871866/35‑1 (MQ=255)               ttCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                <  1:7158773/35‑1 (MQ=38)               ttCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                <  1:6149671/35‑1 (MQ=255)                tCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc               >  1:998290/1‑35 (MQ=255)                 ccaccaACGTTGGCAGACATCAGTGCTGGCAACCg              >  1:517540/1‑35 (MQ=255)                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCAt           <  1:3567838/35‑1 (MQ=255)                    ccaACGTTGGCAGACATCAGTGCTGGCAACCGCAt           <  1:7149962/35‑1 (MQ=255)                     caACGTTGGCAGACATCAGTGCTGGCAACCGCAtt          <  1:4442040/35‑1 (MQ=38)                     caACGTTGGCAGACATCAGTGCTGGCAACCGCAtt          <  1:6977585/35‑1 (MQ=255)                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTc         <  1:1709606/35‑1 (MQ=255)                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTc         <  1:2174890/35‑1 (MQ=255)                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTc         <  1:4291646/35‑1 (MQ=255)                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTc         <  1:3707098/35‑1 (MQ=255)                      aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTc         <  1:7722292/35‑1 (MQ=38)                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTcc        <  1:5933372/35‑1 (MQ=255)                       aCGTTGGCAGACATCAGTGCTGGCAACCGCATTcc        <  1:5502449/35‑1 (MQ=38)                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg      <  1:1388647/35‑1 (MQ=37)                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg      <  1:910453/35‑1 (MQ=255)                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg      <  1:2194326/35‑1 (MQ=37)                         gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg      <  1:6513370/35‑1 (MQ=37)                          ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGc     >  1:2351659/1‑35 (MQ=255)                           tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg    <  1:3975993/35‑1 (MQ=255)                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGt   <  1:3581834/35‑1 (MQ=255)                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGt   >  1:3143173/1‑35 (MQ=255)                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGt   >  1:4370815/1‑35 (MQ=37)                            ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGt   <  1:5583710/35‑1 (MQ=37)                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt  <  1:6965237/35‑1 (MQ=255)                             gCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt  <  1:6395789/35‑1 (MQ=255)                                |                                CATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTT  >  NC_000913/4093192‑4093254

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