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%	179.7 / NA	49	R250H (CGT→CAT)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/21);  total (28/21)

ACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGC  >  NC_000913/1938262‑1938544                                                                                                                                        |                                                                                                                                                    aCTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGgtgt                                                                                                                                      <  2:443385/151‑1 (MQ=255)      gCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGaaa                                                                                                                                 <  1:131739/151‑1 (MQ=255)        ccGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAAtt                                                                                                                               >  1:175596/1‑151 (MQ=255)        ccGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAAtt                                                                                                                               <  1:264394/151‑1 (MQ=255)            cGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCg                                                                                                                           <  2:125803/151‑1 (MQ=255)                 ggCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCg                                                                                                                      <  1:631078/151‑1 (MQ=255)                   cACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCgg                                                                                                                     >  2:385461/1‑150 (MQ=255)                     cgcgATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGTACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAAc                                                                                                                  >  2:396247/1‑151 (MQ=255)                     cgcgATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAAc                                                                                                                  >  1:557329/1‑151 (MQ=255)                         aTGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGt                                                                                                              >  2:794488/1‑151 (MQ=255)                           gCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCg                                                                                                            >  2:7475/1‑151 (MQ=255)                                 tgtgATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTc                                                                                                      >  1:242447/1‑151 (MQ=255)                                    gATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAg                                                                                                    <  2:39857/150‑1 (MQ=255)                                         gAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCg                                                                                              >  2:210981/1‑151 (MQ=255)                                           aaaTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGtt                                                                                            <  2:693546/151‑1 (MQ=255)                                             aTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTg                                                                                           >  1:264399/1‑150 (MQ=255)                                             aTTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTg                                                                                           >  2:50468/1‑150 (MQ=255)                                                gttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATcc                                                                                       >  2:165217/1‑151 (MQ=255)                                                gttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATcc                                                                                       <  2:124770/151‑1 (MQ=255)                                                 ttgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                                                                      >  1:275172/1‑151 (MQ=255)                                                  tgCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                                                                      <  1:468516/150‑1 (MQ=255)                                                     cAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcgtcgt                                                                                  <  2:504796/151‑1 (MQ=255)                                                              aCGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTc                                                                           >  1:252670/1‑149 (MQ=255)                                                                 tGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCt                                                                       >  2:422353/1‑150 (MQ=255)                                                                    ggAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAAc                                                                   <  2:833770/151‑1 (MQ=255)                                                                     gAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAAcc                                                                  >  2:390688/1‑151 (MQ=255)                                                                       aGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGa                                                                <  1:344739/151‑1 (MQ=255)                                                                         ccGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGc                                                              <  2:586302/151‑1 (MQ=255)                                                                          cGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCg                                                             >  2:746805/1‑151 (MQ=255)                                                                           gTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCgg                                                            <  2:625044/151‑1 (MQ=255)                                                                                   ccAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACg                                                    >  1:567892/1‑151 (MQ=255)                                                                                     aGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACg                                                    >  1:671571/1‑149 (MQ=255)                                                                                       gATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCg                                                 >  1:383559/1‑150 (MQ=255)                                                                                        aTGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGAc                                               >  2:551672/1‑151 (MQ=255)                                                                                          gCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGAcc                                              >  2:642056/1‑150 (MQ=255)                                                                                                      atcatcatcGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGt                                 >  1:155242/1‑151 (MQ=255)                                                                                                        catcCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTc                                >  2:496847/1‑150 (MQ=255)                                                                                                        catcCTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTc                                >  1:683906/1‑150 (MQ=255)                                                                                                           ccTCGCCTCTGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTGGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATg                            >  1:232407/1‑151 (MQ=255)                                                                                                                 ctctGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACt                      >  2:564467/1‑151 (MQ=255)                                                                                                                  tctGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACt                      >  1:776380/1‑150 (MQ=255)                                                                                                                   atGAAGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTa                     >  1:725579/2‑150 (MQ=255)                                                                                                                    tGACGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACt                      <  1:35084/148‑1 (MQ=255)                                                                                                                      aCGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATGCAGAAAGCGTTGATCCGTAGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGCCCCAGATGATGGAGTCAATGATTACTa                     <  1:228706/147‑1 (MQ=255)                                                                                                                       cGTGGTAATGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCg                <  2:656278/151‑1 (MQ=255)                                                                                                                              aTGGTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGGGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCg          <  2:843964/150‑1 (MQ=255)                                                                                                                                ggTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGAc        <  2:762597/150‑1 (MQ=255)                                                                                                                                ggTTGCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGACATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGACTACTAACCCGATGCCGAcg       <  2:771369/151‑1 (MQ=255)                                                                                                                                    gCACATGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTg   <  2:103711/151‑1 (MQ=255)                                                                                                                                     cacaTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGc  <  1:590232/151‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                    ACTATGCCCGTCGCCTGGCACGCGATGCAGGATGTGATGCGAAAATTGTTGCCAAGGTTGAACGTGCGGAAGCCGTTTGCAGCCAGGATGCAATGGATGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGCGACCCAGATGATGGAGTCAATGATTACTAACCCGATGCCGACGCGTGC  >  NC_000913/1938262‑1938544

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