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

GTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093178‑4093268                                              |                                              gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           <  1:668469/50‑1 (MQ=255) gTTAGCGGGCATAACATCGGTCATATATTGCCACCAACGTTGGCAGACAt                                           <  1:7003968/50‑1 (MQ=39)  ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           >  1:5452209/1‑49 (MQ=255)  ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                          <  1:4857432/50‑1 (MQ=255)   tAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           >  1:359544/1‑48 (MQ=255)    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                         <  1:1899768/49‑1 (MQ=255)    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                        >  1:6459976/1‑50 (MQ=255)     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                       <  1:367418/50‑1 (MQ=255)     gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                       <  1:2148906/50‑1 (MQ=255)      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                      >  1:1427026/1‑50 (MQ=255)      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                      >  1:4830420/1‑50 (MQ=255)      cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                      <  1:5083812/50‑1 (MQ=255)       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     >  1:2766249/1‑50 (MQ=255)       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     >  1:6511393/1‑50 (MQ=255)       gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                     <  1:811345/50‑1 (MQ=255)        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                    <  1:2982645/50‑1 (MQ=255)        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                    >  1:4260582/1‑50 (MQ=255)         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   <  1:3734514/50‑1 (MQ=255)         gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                   >  1:5147206/1‑50 (MQ=255)          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  <  1:437320/50‑1 (MQ=255)          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  <  1:1513264/50‑1 (MQ=255)          cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  >  1:2873813/1‑50 (MQ=255)           aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                  >  1:6069280/1‑49 (MQ=255)           aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 <  1:1997946/50‑1 (MQ=255)            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:5806823/1‑49 (MQ=255)            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                 >  1:3041367/1‑49 (MQ=255)            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                <  1:1893620/50‑1 (MQ=255)             aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                               >  1:2189856/1‑50 (MQ=255)              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                               <  1:3691082/49‑1 (MQ=255)              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                              <  1:420265/50‑1 (MQ=255)              aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                              <  1:1246328/50‑1 (MQ=255)               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                             >  1:706179/1‑50 (MQ=255)               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                             >  1:3107685/1‑50 (MQ=255)               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                             >  1:7484085/1‑50 (MQ=255)               cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                             >  1:2905591/1‑50 (MQ=255)                aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:1036682/1‑50 (MQ=255)                aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            <  1:899137/50‑1 (MQ=255)                aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:2539779/1‑50 (MQ=255)                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                            >  1:4809779/1‑49 (MQ=255)                 tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGc                           <  1:3353157/50‑1 (MQ=255)                  cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                          >  1:3509208/1‑50 (MQ=255)                  cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                          >  1:2560623/1‑50 (MQ=255)                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         >  1:1877653/1‑50 (MQ=255)                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         <  1:6115291/50‑1 (MQ=255)                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                         >  1:2137059/1‑50 (MQ=255)                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:624318/1‑50 (MQ=255)                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:2669286/1‑50 (MQ=255)                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        <  1:7313196/50‑1 (MQ=255)                    gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                        >  1:2292063/1‑50 (MQ=255)                     tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       >  1:442214/1‑50 (MQ=255)                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       >  1:7006763/1‑49 (MQ=255)                      cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       >  1:4781205/1‑49 (MQ=255)                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCACCCGCATTCCa                     <  1:6007523/50‑1 (MQ=255)                       atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                      <  1:5668759/49‑1 (MQ=255)                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTTGCAACCGCATTCCAg                    >  1:1584694/1‑50 (MQ=255)                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCGACCGCATTCCAg                    <  1:2216170/50‑1 (MQ=255)                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCATCCGCATTCCAg                    <  1:2212170/50‑1 (MQ=255)                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                    >  1:3437288/1‑50 (MQ=255)                        tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                    >  1:700467/1‑50 (MQ=255)                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   >  1:723853/1‑50 (MQ=255)                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   >  1:1348867/1‑50 (MQ=255)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCGACCGCATTCCAGc                   <  1:1173981/49‑1 (MQ=39)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   <  1:944068/49‑1 (MQ=255)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   <  1:7036882/49‑1 (MQ=255)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                   <  1:6718514/49‑1 (MQ=255)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                  >  1:4478160/1‑50 (MQ=255)                          tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                  <  1:5520055/50‑1 (MQ=255)                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                <  1:1298611/50‑1 (MQ=255)                            tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                >  1:7048573/1‑50 (MQ=255)                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct              >  1:1931621/1‑50 (MQ=255)                              tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc               <  1:3523695/49‑1 (MQ=255)                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:2016781/1‑50 (MQ=255)                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct              <  1:2774530/49‑1 (MQ=255)                               ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct              <  1:6219575/49‑1 (MQ=255)                                caccaACGTTGGCAGACATCAGTGCTGGCCACCGCATTCCAGCGttcttc            <  1:3425390/50‑1 (MQ=255)                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc            <  1:5287958/50‑1 (MQ=255)                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:1642698/1‑49 (MQ=255)                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:5918414/1‑46 (MQ=255)                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:6549428/1‑46 (MQ=255)                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt             >  1:2540948/1‑46 (MQ=255)                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa         <  1:6378027/50‑1 (MQ=255)                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa         <  1:4162077/50‑1 (MQ=255)                                   caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa         <  1:3181489/50‑1 (MQ=255)                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt        >  1:6642678/1‑50 (MQ=255)                                    aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt        <  1:738485/50‑1 (MQ=255)                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       >  1:1950007/1‑50 (MQ=255)                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       <  1:5086560/50‑1 (MQ=255)                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       >  1:1019395/1‑49 (MQ=255)                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt       >  1:888252/1‑49 (MQ=255)                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:6662276/1‑49 (MQ=255)                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:6753227/1‑49 (MQ=255)                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:2129967/1‑49 (MQ=255)                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc      >  1:3219963/1‑49 (MQ=255)                                        ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa    >  1:660161/1‑50 (MQ=255)                                         tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt   >  1:7289359/1‑50 (MQ=255)                                         tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt   >  1:3344997/1‑50 (MQ=255)                                          ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc  >  1:6403042/1‑50 (MQ=255)                                          ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt   <  1:3810364/49‑1 (MQ=255)                                          ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt   <  1:4215783/49‑1 (MQ=255)                                          ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt   <  1:3264184/49‑1 (MQ=255)                                              |                                              GTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093178‑4093268

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