Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,938,397	G→A	R250H (CGT→CAT)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,938,397	0	G	A	100.0%	182.9 / NA	47	R250H (CGT→CAT)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/22);  total (25/22)

GGCGAAGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTG  >  NC_000913/1938249‑1938543                                                                                                                                                     |                                                                                                                                                   ggCGAAGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATg                                                                                                                                                  <  1:833806/151‑1 (MQ=255)                 aTGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTgg                                                                                                                                   <  1:139840/150‑1 (MQ=255)                      cGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTg                                                                                                                             <  1:613057/151‑1 (MQ=255)                                     gATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTgg                                                                                                              >  2:665164/1‑151 (MQ=255)                                        gCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGAGGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCg                                                                                                           <  2:556253/151‑1 (MQ=255)                                        gCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCg                                                                                                           >  2:331371/1‑151 (MQ=255)                                          aGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGc                                                                                                         >  2:85505/1‑151 (MQ=255)                                          aGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGc                                                                                                         >  1:22318/1‑151 (MQ=255)                                            gATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAt                                                                                                       >  2:373292/1‑151 (MQ=255)                                            gATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCAt                                                                                                       >  2:805361/1‑151 (MQ=255)                                              tgtgATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTc                                                                                                     >  2:26149/1‑151 (MQ=255)                                                     cGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGc                                                                                              >  1:397858/1‑151 (MQ=255)                                                         aaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTg                                                                                          <  1:837405/151‑1 (MQ=255)                                                                   aaGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcg                                                                                  <  1:510792/149‑1 (MQ=255)                                                                        tGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGcgcg                                                                            <  1:794744/150‑1 (MQ=255)                                                                         gAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTc                                                                          >  1:375340/1‑151 (MQ=255)                                                                          aaCGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCa                                                                         <  1:601934/151‑1 (MQ=255)                                                                                   aaGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCg                                                                <  1:331397/151‑1 (MQ=255)                                                                                     gCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAg                                                              <  2:741336/151‑1 (MQ=255)                                                                                      ccGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGc                                                             >  2:167530/1‑151 (MQ=255)                                                                                      ccGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGc                                                             <  1:116549/151‑1 (MQ=255)                                                                                       cGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCg                                                            <  1:854953/151‑1 (MQ=255)                                                                                        gTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCgg                                                           >  2:144941/1‑151 (MQ=255)                                                                                         tttGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCgg                                                           <  2:386151/150‑1 (MQ=255)                                                                                               gCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCAc                                                    >  1:779253/1‑151 (MQ=255)                                                                                                ccAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACg                                                   <  1:204935/151‑1 (MQ=255)                                                                                                 cAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACgg                                                  <  1:521617/151‑1 (MQ=255)                                                                                                  aGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGc                                                 >  2:512429/1‑151 (MQ=255)                                                                                                   ggATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCg                                                <  1:204907/151‑1 (MQ=255)                                                                                                   ggATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCg                                                >  1:393118/1‑151 (MQ=255)                                                                                                       gCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGAccc                                            >  2:163139/1‑151 (MQ=255)                                                                                                         aaTGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGAccc                                            >  2:452751/1‑149 (MQ=255)                                                                                                              aTGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAgatgat                                     >  2:856396/1‑151 (MQ=255)                                                                                                                 aCATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATgg                                   >  2:843770/1‑150 (MQ=255)                                                                                                                  catcatCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAg                                 <  2:719707/151‑1 (MQ=255)                                                                                                                      atcCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCaa                             <  2:10327/151‑1 (MQ=255)                                                                                                                         cTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATg                           >  1:467665/1‑150 (MQ=255)                                                                                                                               tctGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACt                     >  2:206207/1‑150 (MQ=255)                                                                                                                                 tGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTaa                   >  2:153439/1‑150 (MQ=255)                                                                                                                                  gACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTa                    <  1:14162/148‑1 (MQ=255)                                                                                                                                       ggTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATg            <  2:835832/151‑1 (MQ=255)                                                                                                                                         tAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGcc          >  1:909960/1‑151 (MQ=255)                                                                                                                                           aTGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGa        <  1:82154/151‑1 (MQ=255)                                                                                                                                            tGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCg         >  2:221805/1‑149 (MQ=255)                                                                                                                                             ggTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcg      <  1:777353/151‑1 (MQ=255)                                                                                                                                              gTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAcgc     >  2:455819/1‑151 (MQ=255)                                                                                                                                                 gCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTg  >  2:287198/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGCGAAGATCTGAACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTG  >  NC_000913/1938249‑1938543

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