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

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269                                               |                                               ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:3151758/1‑50 (MQ=255)  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:12235/1‑50 (MQ=255)   ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:1733446/50‑1 (MQ=255)   ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:5695841/50‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:985039/49‑1 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:2121765/1‑50 (MQ=255)     aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:1668005/1‑50 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         <  1:4199413/49‑1 (MQ=255)      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:5052867/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:5385962/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:4692061/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:3622977/50‑1 (MQ=255)       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  1:3008854/1‑50 (MQ=255)         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:4868023/50‑1 (MQ=255)         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:3508938/50‑1 (MQ=255)          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:5229750/50‑1 (MQ=255)          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:6437723/1‑50 (MQ=255)          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:950276/1‑50 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:1240756/50‑1 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:1330700/50‑1 (MQ=255)            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:6417598/50‑1 (MQ=255)             tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  >  1:4576518/1‑49 (MQ=255)             tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                 <  1:1055745/50‑1 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               >  1:497696/1‑50 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               >  1:4437930/1‑50 (MQ=255)               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:6608500/50‑1 (MQ=255)                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:120740/50‑1 (MQ=255)                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             <  1:1013134/50‑1 (MQ=255)                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:6125197/1‑50 (MQ=255)                  tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGc                            <  1:2544990/50‑1 (MQ=255)                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:2469354/1‑50 (MQ=255)                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:1961444/50‑1 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:3664110/1‑50 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:4852837/50‑1 (MQ=255)                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:635745/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCGACCGCATTc                        <  1:5879548/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:5902097/50‑1 (MQ=255)                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        <  1:2011761/50‑1 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:924653/1‑49 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:3465308/1‑49 (MQ=255)                       cATATATTTCCACCAACGTTGGCAGACATCAGCGCTGGCAACCGCATTc                        >  1:1493795/1‑49 (MQ=39)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:5695406/49‑1 (MQ=255)                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:3050239/50‑1 (MQ=255)                         tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:2742303/50‑1 (MQ=255)                         tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:3697050/50‑1 (MQ=255)                         tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     >  1:5252322/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:1646798/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:2230338/49‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:1096194/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:2008701/50‑1 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:4135746/1‑50 (MQ=255)                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:5838657/50‑1 (MQ=255)                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:5839064/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:5861346/1‑50 (MQ=255)                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  <  1:6089749/49‑1 (MQ=255)                               tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:2945169/49‑1 (MQ=255)                               tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:4528612/49‑1 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCCACCGCATTCCAGCGttct               <  1:4935633/49‑1 (MQ=255)                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:288855/1‑50 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAGCCGCATTCCAGCGttctt              <  1:4688458/49‑1 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:1884991/50‑1 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:900050/1‑49 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4192857/1‑49 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:4163272/49‑1 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:791827/1‑49 (MQ=255)                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:701485/49‑1 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:4188337/1‑48 (MQ=255)                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:4785109/50‑1 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:6616944/1‑46 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:3188392/50‑1 (MQ=255)                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:1976172/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:2052515/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:1361657/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:266077/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:4106004/50‑1 (MQ=255)                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:4339656/1‑50 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCGGCCGCATTCCAGCGTTCTTCAGAtt        <  1:1907741/50‑1 (MQ=255)                                      aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:1093495/50‑1 (MQ=255)                                       ctttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:2037898/3‑49 (MQ=39)                                       cGTTGGCAGACATCAGTGCTGGCCGCCGCATTCCAGCGTTCTTCAGATTc       <  1:1313217/50‑1 (MQ=38)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:772010/1‑49 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:2640855/50‑1 (MQ=255)                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:5854514/50‑1 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:3660102/1‑49 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:4574839/50‑1 (MQ=255)                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:467410/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAGCCGCATTCCAGCGTTCTTCAGATTCaa     <  1:3173039/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:2855540/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:1152972/50‑1 (MQ=255)                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:1661781/1‑50 (MQ=255)                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:2362462/1‑50 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:4750656/49‑1 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:912035/49‑1 (MQ=255)                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGACTCAAtc   >  1:2391157/1‑50 (MQ=39)                                            gCAGACATCAGTGCTGGCGACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:1703478/49‑1 (MQ=39)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  >  1:5286796/1‑50 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:3859499/49‑1 (MQ=255)                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:2143535/49‑1 (MQ=255)                                               |                                               GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269

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