Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	801,422	G→A	A222T (GCA→ACA)	ybhH →	putative PrpF family isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	801,422	0	G	A	100.0%	208.0 / NA	59	A222T (GCA→ACA)	ybhH	putative PrpF family isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/36);  total (23/36)

ACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTT  >  NC_000913/801326‑801517                                                                                                 |                                                                                                aCCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACTGGCGGAACTGGATGCCGACAAAACat                                                                                              <  1:797037/100‑1 (MQ=255) aCCTGTATCGATATGGCGAGGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACat                                                                                              <  2:556846/100‑1 (MQ=255)   cTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACatta                                                                                            <  2:307836/100‑1 (MQ=255)    tGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACattat                                                                                           <  2:689713/100‑1 (MQ=255)    tGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACattat                                                                                           >  1:543025/1‑100 (MQ=255)    tGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACattat                                                                                           <  2:54358/100‑1 (MQ=255)     gTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACattatt                                                                                          >  2:554937/1‑100 (MQ=255)      tATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACAttatta                                                                                         >  1:608742/1‑100 (MQ=255)       aTCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGTCGACAAAACATTATTAg                                                                                        >  2:943738/1‑100 (MQ=255)       aTCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAg                                                                                        <  2:525726/100‑1 (MQ=255)          gATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGccc                                                                                     <  1:599313/100‑1 (MQ=255)            tatGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGc                                                                                   <  1:1043206/100‑1 (MQ=255)               ggCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGAATGCCGACAAAACATTATTAGCCCGCAtt                                                                                >  1:1101521/1‑100 (MQ=255)                gCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTg                                                                               <  2:296919/100‑1 (MQ=255)                gCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTg                                                                               <  2:1262542/100‑1 (MQ=255)                  gATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGa                                                                              <  2:1141825/99‑1 (MQ=255)                   aTGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGCTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAAt                                                                            <  2:590303/100‑1 (MQ=255)                       cAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAAGAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTAt                                                                        >  1:917548/1‑100 (MQ=255)                        aGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATc                                                                       <  1:368349/100‑1 (MQ=255)                           cgtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGt                                                                    >  1:1160802/1‑100 (MQ=255)                             tcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCt                                                                  <  1:169204/100‑1 (MQ=255)                              cATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTa                                                                 <  2:985106/100‑1 (MQ=255)                              cATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTa                                                                 >  2:1157815/1‑100 (MQ=255)                                 tattCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACaa                                                              >  1:780365/1‑100 (MQ=255)                                    tCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCg                                                           >  2:368195/1‑100 (MQ=255)                                     ccGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCgg                                                          <  1:364052/100‑1 (MQ=255)                                      cGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAACCCGCATTGAATATATACGTATACAAGCgg                                                          >  2:814435/1‑99 (MQ=255)                                             atatCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCa                                                  >  1:490944/1‑100 (MQ=255)                                                tCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGTCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATg                                               <  1:943746/100‑1 (MQ=255)                                                 cTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATgg                                              <  2:543021/100‑1 (MQ=255)                                                 cTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATgg                                              <  2:105954/100‑1 (MQ=255)                                                   gggAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGc                                            <  1:1148104/100‑1 (MQ=255)                                                   gggAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGc                                            >  1:798757/1‑100 (MQ=255)                                                     gAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCtt                                          >  1:546446/1‑100 (MQ=255)                                                            aGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGAt                                   >  1:1263452/1‑100 (MQ=255)                                                              gTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGt                                 <  1:939635/100‑1 (MQ=255)                                                              gTTATGAATTACCGGCGGAACTGAATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGt                                 <  2:1101511/100‑1 (MQ=255)                                                               ttATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTc                                <  1:554941/100‑1 (MQ=255)                                                                tATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCa                               <  1:78556/100‑1 (MQ=255)                                                                tATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCa                               <  1:262954/100‑1 (MQ=255)                                                                 aTGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAg                              <  2:546442/100‑1 (MQ=255)                                                                    aaTTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACGAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTaa                           >  1:914826/1‑100 (MQ=255)                                                                       tACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAatat                        <  1:436744/100‑1 (MQ=255)                                                                          cggcggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGt                     >  1:930548/1‑100 (MQ=255)                                                                          cggcggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGt                     <  1:1226583/100‑1 (MQ=255)                                                                          cggcggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAAGCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGt                     <  1:1018163/100‑1 (MQ=255)                                                                            gcggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTa                   <  1:334348/100‑1 (MQ=255)                                                                              ggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATc                 >  1:1054302/1‑100 (MQ=255)                                                                               gAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATcc                <  1:1213681/100‑1 (MQ=255)                                                                                aaCTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATccc               <  1:470414/100‑1 (MQ=255)                                                                                aaCTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATccc               <  2:608738/100‑1 (MQ=255)                                                                                   tGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTaa            <  1:747682/100‑1 (MQ=255)                                                                                   tGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTaa            >  1:68730/1‑100 (MQ=255)                                                                                    ggATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTaaa           >  2:418810/1‑100 (MQ=255)                                                                                       tGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCt        >  1:996141/1‑100 (MQ=255)                                                                                       tGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCt        >  1:556330/1‑100 (MQ=255)                                                                                         ccGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgt      >  1:680201/1‑100 (MQ=255)                                                                                          cGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgtg     <  2:917540/100‑1 (MQ=255)                                                                                             cAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCtt  <  2:675614/100‑1 (MQ=255)                                                                                                 |                                                                                                ACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTT  >  NC_000913/801326‑801517

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