Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,150,316	G→A	E219K (GAG→AAG)	yghA →	NADP(+)‑dependent aldehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,316	0	G	A	100.0%	107.8 / NA	32	E219K (GAG→AAG)	yghA	NADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/15);  total (17/15)

AGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTG  >  NC_000913/3150077‑3150561                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      aGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGCCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTAtt                                                                                                                                                                                                                                            >  1:179866/1‑251 (MQ=255)          ttCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCACAACAGGTCGCGAAGAAAGGTATTCGGGTGAat                                                                                                                                                                                                                                   >  2:10051/1‑251 (MQ=255)              ggATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGCCCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCGACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACGGGTCGCGAAGAAAGGTATTCGGGGGAATATTg                                                                                                                                                                                                                               <  2:145323/251‑1 (MQ=255)                      cAGACCTCACCAGCGAACAGTTTCACAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGcc                                                                                                                                                                                                                        >  2:202799/1‑250 (MQ=255)                               ccAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATc                                                                                                                                                                                                              >  2:135942/1‑251 (MQ=255)                                    gAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTGCCGAAAGGTGCAAGTATCATCACCACTTCGGCAACCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGTCTACGAAGGCGGCGATTCTGAACTACAGGCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCTCGTCAGGCCTGTTCTGGAc                                                                                                                                                                                                         <  2:179866/251‑1 (MQ=255)                                         gTTTCAAAAGACATTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCAc                                                                                                                                                                                                    >  1:210318/1‑251 (MQ=255)                                                       ttGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcgg                                                                                                                                                                                       <  2:177103/250‑1 (MQ=255)                                                                aCGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCATCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACg                                                                                                                                                                             >  2:153905/1‑251 (MQ=255)                                                                  gTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGc                                                                                                                                                                            >  1:244589/1‑250 (MQ=255)                                                                          gcTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATa                                                                                                                                                                      >  1:213889/1‑248 (MQ=255)                                                                            tGTTCTGGCTAACCCAGGAAGCGATCCACCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTAATGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGATCGCGAAGACAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATa                                                                                                                                                                      >  2:144526/1‑246 (MQ=255)                                                                                   gCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTCCTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGt                                                                                                                                                          >  1:184424/1‑251 (MQ=255)                                                                                      aaCCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGttt                                                                                                                                                        <  2:244589/250‑1 (MQ=255)                                                                                      aaCCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTg                                                                                                                                                       <  1:268212/251‑1 (MQ=255)                                                                                       aCCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTgg                                                                                                                                                      <  2:57376/251‑1 (MQ=255)                                                                                          cAGGAAGCGATCCCCCTGCTAACGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGTTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCa                                                                                                                                                   <  1:117981/251‑1 (MQ=255)                                                                                                      ccccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGAt                                                                                                                                       <  2:51781/251‑1 (MQ=255)                                                                                                           gCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAAc                                                                                                                                  <  1:256709/251‑1 (MQ=255)                                                                                                              aCCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACGGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCGGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTGTGGTCAGCAAACGCCGATGAAACGt                                                                                                                                <  2:141463/250‑1 (MQ=255)                                                                                                               ccGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGc                                                                                                                              >  1:160474/1‑251 (MQ=255)                                                                                                                cGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCg                                                                                                                             >  2:303566/1‑251 (MQ=255)                                                                                                                cGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCg                                                                                                                             <  1:156670/251‑1 (MQ=255)                                                                                                                                            gTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGACTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGtat                                                                                                 <  1:248084/251‑1 (MQ=255)                                                                                                                                                 tCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGtt                                                                                            >  1:301098/1‑251 (MQ=255)                                                                                                                                                                           ttACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCt                                                                   <  1:153405/250‑1 (MQ=255)                                                                                                                                                                                   cTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCg                                                          >  2:1258/1‑251 (MQ=255)                                                                                                                                                                                            tACGAAGGCGTCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCCAATTTCCGGCGGACAAACGCAGGCTAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAgagc                                                 >  2:255627/1‑248 (MQ=255)                                                                                                                                                                                                                       ccGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAg                      <  1:202799/251‑1 (MQ=255)                                                                                                                                                                                                                            gCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCCGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTa                  >  2:260060/1‑250 (MQ=255)                                                                                                                                                                                                                                        aGGTCGCGAAGAAAGGTATTCGGGAGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGAGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGCGCATTTAGGTTaaaacaaatcccgg     >  2:61729/1‑241 (MQ=255)                                                                                                                                                                                                                                            cgcgAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTtg  <  1:168280/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      AGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTG  >  NC_000913/3150077‑3150561

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