Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,150,497	A→G	Y279C (TAC→TGC)	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,497	0	A	G	100.0%	101.0 / NA	29	Y279C (TAC→TGC)	yghA	NADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/13);  total (16/13)

TGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGAAAACGCCTTAT  >  NC_000913/3150368‑3150622                                                                                                                                  |                                                                                                                              tGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGt                                                                                                                            <  1:112445/133‑1 (MQ=255)        ttCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCa                                                                                                                     >  1:424870/1‑133 (MQ=255)          ccGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGa                                                                                                                   <  2:150062/133‑1 (MQ=255)              cggACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTg                                                                                                               <  2:134094/133‑1 (MQ=255)                  cAAACGAAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACg                                                                                                           <  2:141915/133‑1 (MQ=255)                               aGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTgcggcgg                                                                                              >  1:93062/1‑133 (MQ=255)                               aGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTgcggcgg                                                                                              >  1:523980/1‑133 (MQ=255)                                 aTCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTgcggcggcg                                                                                            <  2:93062/133‑1 (MQ=255)                                 aTCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTgcggcggcg                                                                                            <  2:523980/133‑1 (MQ=255)                                     cGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCa                                                                                        >  1:95696/1‑133 (MQ=255)                                     cGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGCGCACGGCGAGTGCGGCGGCGAGCa                                                                                        >  2:250583/1‑133 (MQ=255)                                                   cAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTaaaaaa                                                                          <  2:424870/133‑1 (MQ=255)                                                        gCCGATGATACGTGCCGGTCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGcc                                                                     <  1:171926/133‑1 (MQ=255)                                                        gCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGGGGCGGCGGGCAGTTTGGTTAAGAAAATGcc                                                                     >  1:169579/1‑133 (MQ=255)                                                          cGATGAAACGTGCGGTGCAACCGGTGGACCTGTCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCg                                                                   <  2:169579/133‑1 (MQ=255)                                                                            aaCCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCgg                                                 >  1:432738/1‑133 (MQ=255)                                                                                      aCTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATtgt                                       >  1:622357/1‑133 (MQ=255)                                                                                               tGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAAtt                                  >  1:403383/1‑129 (MQ=255)                                                                                               tGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAAtt                                  <  2:403383/129‑1 (MQ=255)                                                                                               tGTATATGGTTATCTGGCTAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCAct                              <  1:260346/133‑1 (MQ=255)                                                                                                         tATCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCAGTCTTACCTGAt                    >  1:323770/1‑133 (MQ=255)                                                                                                           tCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGAtgt                  >  1:462897/1‑133 (MQ=255)                                                                                                           tCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGAtgt                  >  1:431420/1‑133 (MQ=255)                                                                                                           tCTGGCAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGAtgt                  >  2:523660/1‑133 (MQ=255)                                                                                                                cAAGTCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGa             >  1:502222/1‑133 (MQ=255)                                                                                                                 aaGTCAGGCGTCGAGCTGTGTCCCCGCAGATGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGaa            <  1:250583/133‑1 (MQ=255)                                                                                                                    tCAGGAGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGAAAAc         >  2:555609/1‑133 (MQ=255)                                                                                                                         aGTCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGAAAACGCCtt    >  1:408970/1‑133 (MQ=255)                                                                                                                           tCGAGCTGCGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGAAAACGCCTTAt  <  2:502222/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                              TGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTTGTGAAAAGCAACCGGGCATCATTGTGAATTCACTCTTACCTGATGTGGCGAAAACGCCTTAT  >  NC_000913/3150368‑3150622

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