Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,495,194	A→T	F339Y (TTT→TAT)	ahr ←	NADPH‑dependent aldehyde reductase Ahr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,495,194	0	A	T	100.0%	58.8 / NA	20	F339Y (TTT→TAT)	ahr	NADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (10/10);  total (10/10)

GGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCG  >  NC_000913/4494997‑4495422                                                                                                                                                                                                      |                                                                                                                                                                                                                                     ggTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTgg                                                                                                                                                                                 >  2:98006/1‑251 (MQ=255)                          gAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCCGGCTGTGCTGGCGATACGACAAAACAGACTATGTGCGCAAGAGGGCAGCGCCTCAGATCAGCGCTGCGACTGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCg                                                                                                                                                        >  1:223841/1‑251 (MQ=255)                                cTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGCCATCGGGAAc                                                                                                                                                   >  2:74002/1‑250 (MQ=255)                                         ccACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTGCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCggtgg                                                                                                                                         <  1:44825/251‑1 (MQ=255)                                                       acgggATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTtg                                                                                                                           <  1:6721/248‑1 (MQ=255)                                                                tgtgCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCATATCAGCGCTGCGAATGATTTACAATAATCGGCTTTTAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTgcgggc                                                                                                                   <  1:138653/250‑1 (MQ=255)                                                                   gcgcATTGATATCAACCATGATATTGCCAATGTCACCGAACAAGGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTTGATGGCGTCGTTAATTTTCGACATCGGGGACAGTTCGGTGGTCGGCGCAACCTTGCTGCGggcggc                                                                                                                <  2:45704/250‑1 (MQ=255)                                                                                 aCCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCtt                                                                                                 >  1:102316/1‑251 (MQ=255)                                                                                         atTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGCTGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAACACGCATCAGCTTACGCAGCt                                                                                         >  2:206798/1‑251 (MQ=255)                                                                                                            aaTGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAACGAGGGCAGCGCCTCAGATCAGCGCTGCGACTGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGc                                                                      >  1:244552/1‑251 (MQ=255)                                                                                                                 cGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGaa                                                                  <  1:74002/250‑1 (MQ=255)                                                                                                                                                gAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCTATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGATCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGt                                  >  1:253246/1‑251 (MQ=255)                                                                                                                                                         cGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTGTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACAGCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGCACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCgg                         <  1:32825/251‑1 (MQ=255)                                                                                                                                                         cGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATAGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGCGACGCTGCGATCGCCCGCAATTAGAGTAAAGGCCgg                         >  2:237730/1‑251 (MQ=255)                                                                                                                                                          gAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCgg                         >  2:208701/1‑250 (MQ=255)                                                                                                                                                                   cAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgc               <  1:34649/251‑1 (MQ=255)                                                                                                                                                                    aGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgcg              >  1:79243/1‑251 (MQ=255)                                                                                                                                                                           cAGATCGGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTg         <  2:8998/249‑1 (MQ=255)                                                                                                                                                                            aGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTgaga      <  1:194752/251‑1 (MQ=255)                                                                                                                                                                                cAGCGCTGCGAATGATTTTCAATAAGCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGTAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCg  <  2:215164/251‑1 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                                                     GGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCG  >  NC_000913/4494997‑4495422

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