Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,080,037	C→A	L533L (CTG→CTT)	tktA ←	transketolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,080,037	0	C	A	100.0%	176.4 / NA	50	L533L (CTG→CTT)	tktA	transketolase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/22);  total (28/22)

TGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACG  >  NC_000913/3079959‑3080111                                                                               |                                                                           tGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAg                                                                          >  2:534937/1‑81 (MQ=255)    aCCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCAACGCGCGCGATGTTTGCAAGTTg                                                                       >  1:534936/1‑81 (MQ=255)     ccGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGc                                                                      <  1:534884/81‑1 (MQ=255)      cGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGct                                                                     <  2:534886/81‑1 (MQ=255)      cGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGct                                                                     <  2:534885/81‑1 (MQ=255)       ggTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGctc                                                                    <  1:534887/81‑1 (MQ=255)        gTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGctct                                                                   >  2:534945/1‑81 (MQ=255)           gCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTc                                                                 <  1:534888/80‑1 (MQ=255)                  aaaTCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAgttcgtt                                                         >  2:534941/1‑81 (MQ=255)                    aTCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTcc                                                       <  1:534889/81‑1 (MQ=255)                      cAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCctg                                                     >  2:534935/1‑81 (MQ=255)                      cAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCctg                                                     >  1:534975/1‑81 (MQ=255)                       aGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCctgc                                                    >  1:534909/1‑81 (MQ=255)                          tCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCtgctgc                                                 >  1:534951/1‑81 (MQ=255)                           ccGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTgcg                                                >  1:534954/1‑81 (MQ=255)                            cGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTgcg                                                <  1:534890/80‑1 (MQ=255)                              gCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTgcgc                                               >  2:689673/1‑79 (MQ=255)                              gCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCa                                             <  1:534891/81‑1 (MQ=255)                               cTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAg                                            >  1:534938/1‑81 (MQ=255)                               cTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAg                                            <  2:534892/81‑1 (MQ=255)                               cTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAg                                            >  1:534969/1‑81 (MQ=255)                                 gACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGt                                          >  1:534942/1‑81 (MQ=255)                                  aCCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAgg                                           <  1:534893/79‑1 (MQ=255)                                    cGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCt                                       >  2:534960/1‑81 (MQ=255)                                      gcgcAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGa                                     <  2:534894/81‑1 (MQ=255)                                        gcAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACg                                   >  2:534964/1‑81 (MQ=255)                                        gcAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACg                                   <  1:534895/81‑1 (MQ=255)                                         cAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACgg                                  <  1:534896/81‑1 (MQ=255)                                          aGTCTTTCAGCACATAACCACCGCGCGCGATATTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACggg                                 <  2:534897/81‑1 (MQ=255)                                            tCTTTCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGgag                               <  1:534898/81‑1 (MQ=255)                                               ttCAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGagag                             >  2:534958/1‑80 (MQ=255)                                                 cAGCACATAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGAt                          >  1:534962/1‑81 (MQ=255)                                                      caTAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTg                     >  2:534952/1‑81 (MQ=255)                                                       aTAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGc                    >  1:534965/1‑81 (MQ=255)                                                        tAACCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCg                   <  2:534899/81‑1 (MQ=255)                                                         aaCCACCGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCgg                  >  2:534949/1‑81 (MQ=255)                                                         aaCCACCGCGCGCGATGTTTGCAAGTTGCTCCTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCgc                  >  1:534959/1‑80 (MQ=255)                                                          accaccGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGt                 >  2:534950/1‑81 (MQ=255)                                                          accaccGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGt                 <  2:534900/81‑1 (MQ=255)                                                           ccaccGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTc                <  1:534901/81‑1 (MQ=255)                                                           ccaccGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTc                <  1:534902/81‑1 (MQ=255)                                                             accGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCgg              <  1:534903/81‑1 (MQ=255)                                                             accGCGCGCGATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCgg              >  2:534971/1‑81 (MQ=255)                                                               cgcgcgcgATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGc            >  2:534970/1‑81 (MQ=255)                                                                  gcgcgATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGt         <  1:534904/81‑1 (MQ=255)                                                                  gcgcgATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGt         >  2:534955/1‑81 (MQ=255)                                                                     cgATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCt      >  1:534948/1‑81 (MQ=255)                                                                      gATGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCt      >  1:534967/1‑80 (MQ=255)                                                                       aTGTTTGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGa    >  2:534956/1‑81 (MQ=255)                                                                          tttGCAAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACg  <  1:534905/80‑1 (MQ=255)                                                                               |                                                                           TGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACG  >  NC_000913/3079959‑3080111

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