Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	123,050	A→T	I12F (ATC→TTC)	aceE →	pyruvate dehydrogenase E1 component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	123,050	0	A	T	100.0%	124.0 / NA	47	I12F (ATC→TTC)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (25/22);  total (25/22)

AAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAG  >  NC_000913/122910‑123176                                                                                                                                              |                                                                                                                                aatcgcgGCAACTAAACGTAGAACCTGTCTTAGGGAGCTTTCCAGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACCCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTc                                                                                                                               <  1:389201/140‑1 (MQ=255) aaGCGCGGCACCGAACCGTAGAGCCTGTCTTTTTCAGCGTTACAG‑AAAAGTTTATGGTTGCCAGGTTGCCGGAGTCTCAACGTGATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTc                                                                                                                               <  2:630244/143‑1 (MQ=255) aaGCGCGGCAACTAAGCGTAGAACCCGTCTTAGTGGGCTCTGCGGCTGGAGTTCAATG‑GGACAGGTTCCAGAAATCTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTc                                                                                                                               <  2:767993/143‑1 (MQ=255)      cgGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTCCGAAAc                                                                                                                          >  2:269120/1‑143 (MQ=255)            cgagaCGTAGAACCCGTCGTAGTGAGCTTTGCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGa                                                                                                                    <  2:378110/139‑1 (MQ=255)            cTAAACTTATAACCCCTCTTATTTAGCTTTCCGGCTATAGTTCAATG‑GGACAAGTTCCAGAAACCTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGCTCCGTTCGAAACTCGCGa                                                                                                                    <  2:614422/143‑1 (MQ=255)            cTAAACGTAGAACCTTTCTTATTGATCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGATCGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGa                                                                                                                    <  1:893681/143‑1 (MQ=255)            cTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTCCGAAACTCGCGa                                                                                                                    >  2:168183/1‑143 (MQ=255)               aaCGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTAGAAACTCGCGACTg                                                                                                                 >  1:70241/1‑143 (MQ=255)                  gTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGCATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGAATGGCt                                                                                                              >  1:314940/1‑143 (MQ=255)                    aGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAAGCCATGTCAGAACGTTTCCCAAATGAAGTGGATTCGTTCGAAACGCACGACTGGCTcc                                                                                                            >  2:387165/1‑143 (MQ=255)                       aCCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAgg                                                                                                         <  2:113534/143‑1 (MQ=255)                        ccTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGc                                                                                                        >  1:266597/1‑143 (MQ=255)                         cTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACCGGTTCCAGAAAACTCAACGCTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGCGGATCCGTTCGAAACCCGCGACTGGCTCCAGGCg                                                                                                       >  2:260843/1‑143 (MQ=255)                                agtGAGCTTTGCGGCGAGCGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAGCCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAAt                                                                                                <  2:695254/141‑1 (MQ=255)                                aTTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAAt                                                                                                <  1:75011/143‑1 (MQ=255)                                 ttGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGGATc                                                                                               >  1:322853/1‑143 (MQ=255)                                 ttGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATc                                                                                               >  1:266407/1‑143 (MQ=255)                                 ttGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTACCGTCTGGCACCAGCAGATAGaaac                                                                                               >  2:656000/1‑141 (MQ=255)                                    aGCTTTCCGGCGAGGCTTCAATG‑GGACGGGTTCCAGAAAACTCAACGTTATTAGCTAGATAAGGAATAACCCATGTCAGACCGTTTCCCAGATGACGTTGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGt                                                                                            <  1:107546/143‑1 (MQ=255)                                          ccGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCGGGCGATCGAATCGGTCATCCg                                                                                      >  1:472047/1‑143 (MQ=255)                                            ggCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑g                                                                                   >  1:613567/1‑143 (MQ=255)                                             gCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGGCTGGCTCCAGGCGATCGAATAGGTCAGCCggac                                                                                   >  2:857999/1‑140 (MQ=255)                                                agagTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑gaaga                                                                               >  1:842167/1‑143 (MQ=255)                                                  agTTCAAGG‑GGACAGGTCCGGGAAAACTCTACGTTATTAGATAGATCGGGGATGGCCCATGCCACCACGTTCCCAAAATTACGTGGCTCCGTTCGAATCTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑Gaagaag                                                                             <  2:266189/143‑1 (MQ=255)                                                       aaTG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAGCGGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGCAGGTGtt                                                                        >  1:124862/1‑143 (MQ=255)                                                       aaTG‑GGACAGGGTCCAGAAACCTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGtt                                                                        <  1:244884/143‑1 (MQ=255)                                                                aGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGc                                                                >  1:419237/1‑143 (MQ=255)                                                                   ttCGAGAGAACTCAGCGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCa                                                             <  2:19777/143‑1 (MQ=255)                                                                            aCTCAACGTTAGTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGAt                                                    <  1:195678/143‑1 (MQ=255)                                                                             cTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGCTc                                                   >  2:164432/1‑143 (MQ=255)                                                                              tCAACGTTATTAGAGAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCg                                                  <  1:429479/143‑1 (MQ=255)                                                                                aaCGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGAc                                                >  1:280038/1‑143 (MQ=255)                                                                                aaCGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGAAGTGGATCCGTTCGAAACTCACGACTGACTCCAGGCGATCAAATCGGTCATCCGA‑GAAGAAGGAGATGAGCATAATCAGTATCTGATCGAc                                                >  2:74803/1‑143 (MQ=255)                                                                                 aCGTTATTAGATAGCTAAGGAATAACCCATGTCAGACCGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATAAGT‑CTAGATGTTTTTGAGCGTGCTCAGTATCTGATCGAcc                                               <  2:272917/143‑1 (MQ=255)                                                                                 aCGTTATTAGACAGATACGGAATGACCCATGTCCGAACCGTTCCCACAGGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGTATCGGTCATCAGTAGGGTATTATCCTG‑GCGTGCTCAGTATCTGATCGAcc                                               <  2:730527/143‑1 (MQ=255)                                                                                    ttattaGATAGATAAGGAATAACCCAGGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAAc                                            <  1:470398/143‑1 (MQ=255)                                                                                             agataAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATTCGTTCGAAACTCGCGAATGGATAGGGAGTAGCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCCACTGCTTGCTg                                   >  1:652202/1‑143 (MQ=255)                                                                                               ataAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGaa                                 >  1:145010/1‑143 (MQ=255)                                                                                                          cgcATGTCAGAACGTTGCCCAAATGACGTAGATCCGTTCGAAACGAGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAgg                      <  2:597614/141‑1 (MQ=255)                                                                                                              tGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGgt                  <  1:269120/143‑1 (MQ=255)                                                                                                                   gAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTaaa             >  1:432914/1‑143 (MQ=255)                                                                                                                           cccAAATTACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCg     <  1:257967/143‑1 (MQ=255)                                                                                                                           cccAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGACCGTGCTCAGACTCTGATCCGCCAACGGCTTGCTGAAGCCCGCAACGGCGGAGTAAACGtgcccg     >  2:718676/1‑138 (MQ=255)                                                                                                                           ccaaaaTGACGTGGATACGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCAT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCg     <  1:779691/140‑1 (MQ=255)                                                                                                                              aaaTGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGGATCGGTCATCCGT‑GCAGGAGGTGTTGAGCGTGGTCAGTATCCGATCGGCCAACTGCTGGCTGAAGCCCGCAGAGGCGGTGGCAACGTAGcaggag  >  2:695517/1‑138 (MQ=255)                                                                                                                              aaaTGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGccg  >  2:695870/1‑141 (MQ=255)                                                                                                                                              |                                                                                                                                AAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATG‑GGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGT‑GAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAG  >  NC_000913/122910‑123176

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