Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,080,959	C→T	R226H (CGC→CAC)	tktA ←	transketolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,080,959	0	C	T	100.0%	159.2 / NA	47	R226H (CGC→CAC)	tktA	transketolase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/21);  total (26/21)

GCCCAGCGGCGCACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAA  >  NC_000913/3080838‑3081095                                                                                                                          |                                                                                                                                         gCCCAGCGGCGCACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGt                                                                                                                        >  2:32286/1‑140 (MQ=255)  cccAGCGGCGCACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTACACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTc                                                                                                                       >  2:425293/1‑140 (MQ=255)          cgcACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGt                                                                                                               >  2:428787/1‑140 (MQ=255)           gcACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGCTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTc                                                                                                              >  1:87878/1‑140 (MQ=255)            cACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCg                                                                                                             <  1:855909/140‑1 (MQ=255)                 tGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTc                                                                                                        >  1:730567/1‑140 (MQ=255)                 tGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTc                                                                                                        >  1:497968/1‑140 (MQ=255)                    gAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTcgcg                                                                                                     >  2:917704/1‑140 (MQ=255)                      gTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGaa                                                                                                   >  1:748099/1‑140 (MQ=255)                       tCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAAt                                                                                                  >  1:24816/1‑140 (MQ=255)                                      ttGTTCGGGGAACCGAAACCGATGAGGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTa                                                                                   <  1:261457/140‑1 (MQ=255)                                       tGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTaa                                                                                  >  2:363122/1‑140 (MQ=255)                                              ggAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGa                                                                           >  1:332224/1‑140 (MQ=255)                                               gAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGaa                                                                          <  2:300296/140‑1 (MQ=255)                                                  ccgaaaccgaTGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACg                                                                       <  1:71694/140‑1 (MQ=255)                                                  ccgaaaccgaTGATGGTTTGGCACATCAGCAGGAAAGGTTTTTCAGTCACTGCGCGCGCTTCTGCAACTGCGTGTTTGATAGATGCCGGGGCATGGCCGTCGATGTCGCGAATAACGGGCCAGCCGTAAGCTTCGAAACg                                                                       <  2:87878/140‑1 (MQ=255)                                                             gatgGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGgtgt                                                            >  2:52304/1‑140 (MQ=255)                                                               tgGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTgtcg                                                          >  2:636303/1‑140 (MQ=255)                                                                ggTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTgtcgt                                                         >  1:246857/1‑140 (MQ=255)                                                                       cacaTCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGaa                                                  <  2:797359/140‑1 (MQ=255)                                                                        acaTCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAAc                                                 >  2:75720/1‑140 (MQ=255)                                                                         caTCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAAcc                                                <  1:917704/140‑1 (MQ=255)                                                                           tCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAg                                              >  2:503314/1‑140 (MQ=255)                                                                            cagcagGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAAc                                                 >  1:737711/1‑136 (MQ=255)                                                                            cagcagGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGc                                             >  1:1022845/1‑140 (MQ=255)                                                                            cagcagGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGCCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAAc                                                 <  2:737711/136‑1 (MQ=255)                                                                                    aaGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTg                                                    <  2:808282/125‑1 (MQ=255)                                                                                    aaGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTg                                                    >  1:808282/1‑125 (MQ=255)                                                                                       gTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGGCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGa                                  <  1:636303/140‑1 (MQ=255)                                                                                        tttGTCAGTCACTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGAc                                 <  1:503314/140‑1 (MQ=255)                                                                                            tcagtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATc                             <  1:898406/140‑1 (MQ=255)                                                                                            tcagtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATc                             <  2:748099/140‑1 (MQ=255)                                                                                            tcagtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATc                             <  2:24816/140‑1 (MQ=255)                                                                                               gtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGa                                  >  1:432953/1‑132 (MQ=255)                                                                                               gtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGa                                  <  2:432953/132‑1 (MQ=255)                                                                                               gtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGAt                          <  1:52304/140‑1 (MQ=255)                                                                                               gtcaCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGAt                          <  1:379715/140‑1 (MQ=255)                                                                                                 caCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAg                        >  2:508186/1‑140 (MQ=255)                                                                                                 caCTGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAg                        >  2:286381/1‑140 (MQ=255)                                                                                                    tGCGCGCGCTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGaaa                     >  2:902991/1‑140 (MQ=255)                                                                                                        cgcgcTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATc                             >  2:267986/1‑128 (MQ=255)                                                                                                        cgcgcTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATc                             <  1:267986/128‑1 (MQ=255)                                                                                                        cgcgcTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATAcc                 <  2:246857/140‑1 (MQ=255)                                                                                                         gcgcTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCg                >  1:758300/1‑140 (MQ=255)                                                                                                          cgcTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGt               >  1:618543/1‑140 (MQ=255)                                                                                                          cgcTTCTTCTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGt               <  1:286381/140‑1 (MQ=255)                                                                                                             ttcttcTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGt            <  2:497968/140‑1 (MQ=255)                                                                                                             ttcttcTACTGCGTGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGt            <  1:902991/140‑1 (MQ=255)                                                                                                                       gcgtGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCg       <  2:233491/135‑1 (MQ=255)                                                                                                                       gcgtGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCg       >  1:233491/1‑135 (MQ=255)                                                                                                                       gcgtGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGaa  >  2:393885/1‑140 (MQ=255)                                                                                                                          |                                                                                                                                         GCCCAGCGGCGCACCGTGGGAGTCGTGGGTACCGGCTTTGTTCGGGGAACCGAAACCGATGATGGTTTTGCACATCAGCAGGGAAGGTTTGTCAGTCACTGCGCGCGCTTCTTCTACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAA  >  NC_000913/3080838‑3081095

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