Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,676,640	T→G	K422Q (AAA→CAA)	pntA ←	pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,676,640	0	T	G	100.0%	209.0 / NA	56	K422Q (AAA→CAA)	pntA	pyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/31);  total (25/31)

TGCAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTG  >  NC_000913/1676492‑1676777                                                                                                                                                     |                                                                                                                                          tGCAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCg                                                                                                                                         >  2:164183/1‑151 (MQ=255)   cAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGgc                                                                                                                                       >  1:766219/1‑151 (MQ=255)    aaCAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTgcggcg                                                                                                                                      <  1:755744/151‑1 (MQ=255)     aCAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGgcgc                                                                                                                                     >  1:871917/1‑151 (MQ=255)      cAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGgcgc                                                                                                                                     >  2:571524/1‑150 (MQ=255)      cAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCa                                                                                                                                    <  1:103537/151‑1 (MQ=255)       aGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCaa                                                                                                                                   <  1:550755/151‑1 (MQ=255)        gTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAAc                                                                                                                                  <  2:723470/151‑1 (MQ=255)             ccGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCt                                                                                                                               <  1:560892/149‑1 (MQ=255)             ccGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCt                                                                                                                               <  1:490452/149‑1 (MQ=255)               gACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGcc                                                                                                                           >  1:879264/1‑151 (MQ=255)                acaacaATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGcc                                                                                                                           >  1:291082/1‑150 (MQ=255)                   acaaTAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATcc                                                                                                                       >  1:753824/1‑151 (MQ=255)                    caaTAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCa                                                                                                                      <  1:652929/151‑1 (MQ=255)                        aatCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCa                                                                                                                  <  2:761702/151‑1 (MQ=255)                                aaaTCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAAt                                                                                                          >  2:433207/1‑151 (MQ=255)                                aaaTCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAAt                                                                                                          <  1:405632/151‑1 (MQ=255)                                 aaTCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATg                                                                                                         <  2:377996/151‑1 (MQ=255)                                  aTCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAAt                                                                                                          <  2:263181/149‑1 (MQ=255)                                    cgcgTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGAtt                                                                                                      <  2:429947/151‑1 (MQ=255)                                          ggTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAgc                                                                                                <  1:730789/151‑1 (MQ=255)                                          ggTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAg                                                                                                 <  2:543041/150‑1 (MQ=255)                                           gTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAg                                                                                                 <  1:575958/149‑1 (MQ=255)                                                  aCATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCa                                                                                         <  1:224703/150‑1 (MQ=255)                                                   cATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAc                                                                                       <  2:245890/151‑1 (MQ=255)                                                    aTCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAcg                                                                                      <  1:309464/151‑1 (MQ=255)                                                    aTCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAcg                                                                                      <  1:301480/151‑1 (MQ=255)                                                     tCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAAcgc                                                                                     >  2:338111/1‑151 (MQ=255)                                                        aCGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGt                                                                                   >  2:579533/1‑150 (MQ=255)                                                          ggTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTAtt                                                                                <  1:809623/151‑1 (MQ=255)                                                           gtgtATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTAttt                                                                               <  2:892188/151‑1 (MQ=255)                                                             gtATGGAGCGCGTGCGATACATTCCACACGACGTAATAACCGACAGCGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAACAAGAATGATTGCCAGCGCCATCAACGCGTATTTa                                                                              <  2:214848/150‑1 (MQ=255)                                                                  cAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCAc                                                                        >  1:464141/1‑151 (MQ=255)                                                                        gTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAg                                                                  <  2:763108/151‑1 (MQ=255)                                                                            gATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAgg                                                              <  1:3285/151‑1 (MQ=255)                                                                              tACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGt                                                             >  1:103964/1‑150 (MQ=255)                                                                               aCATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTAc                                                           <  1:906847/151‑1 (MQ=255)                                                                                    ccACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACAttt                                                       <  1:366710/150‑1 (MQ=255)                                                                                     cacaCCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACAttttt                                                     >  2:894619/1‑151 (MQ=255)                                                                                            cGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGt                                              <  1:421038/151‑1 (MQ=255)                                                                                               aataaCCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGtttt                                           >  1:821130/1‑151 (MQ=255)                                                                                               aataaCCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGtttt                                           >  2:379407/1‑151 (MQ=255)                                                                                                 taaCCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGGATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTc                                          >  2:95760/1‑150 (MQ=255)                                                                                                   aCCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACt                                       <  2:422999/151‑1 (MQ=255)                                                                                                   aCCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACt                                       <  2:732508/151‑1 (MQ=255)                                                                                                      gACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACtt                                      <  1:433217/149‑1 (MQ=255)                                                                                                         aaCGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGt                                 >  1:82435/1‑151 (MQ=255)                                                                                                             cAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGcc                              >  2:705514/1‑150 (MQ=255)                                                                                                             cAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCg                             <  2:782163/151‑1 (MQ=255)                                                                                                             cAGGCCAGCGCGAAAACGGTGAAGTGCCAAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCg                             >  2:293878/1‑151 (MQ=255)                                                                                                                    gcgcGAAAACGGTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtg                      >  2:888198/1‑151 (MQ=255)                                                                                                                              ggTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGc            >  1:663156/1‑151 (MQ=255)                                                                                                                              ggTGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGc            >  1:328630/1‑151 (MQ=255)                                                                                                                                tGAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCgg          >  1:148037/1‑151 (MQ=255)                                                                                                                                 gAAGTGCCCAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGgc         <  1:654028/151‑1 (MQ=255)                                                                                                                                        ccAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGgctgagctg  >  2:711476/1‑151 (MQ=255)                                                                                                                                         cAAGGAATTCTTGCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGgctgagctg  >  1:223053/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                          TGCAACAGTGCTCCGACAACAATAATCCCTGAAATCGCGTTGGTGACCGACATCAACGGTGTATGCAGCGCGTGCGATACATTCCACACCACGTAATAACCGACAACGCAGGCCAGCGCGAAAACGGTGAAGTGCCCAAGGAATTCTTTCGGCGCAACGCTTGCCATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTG  >  NC_000913/1676492‑1676777

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