Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,730,552	C→T	A489V (GCC→GTC)	lhr →	ATP‑dependent helicase Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,730,552	0	C	T	90.9%	18.0 / NA	11	A489V (GCC→GTC)	lhr	ATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (3/7);  minor base G (0/1);  total (3/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCTGGCGCTCAATTGTTGGCGGTTACCAGCGGCGGCACCATTCCGGATC  >  NC_000913/1730311‑1730740                                                                                                                                                                                                                                                  |                                                                                                                                                                                             gggTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGAGGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTg                                                                                                                                                                                     <  1:104578/251‑1 (MQ=255)                       ccGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGAAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGAAAGATc                                                                                                                                                                             >  1:60941/1‑237 (MQ=255)                                                 aaaGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCACATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGTATGCATTACAGGGAGACGGATGGTACTCCCGCGTACGCCGTGGCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGc                                                                                                                                         <  1:221781/247‑1 (MQ=255)                                                     ggCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCCCCGTCATTGTAGAGTGTATGTTCCCAGGCACGCTGGAAAAACTGACAACACCGCCTAATCCTCTCGAAGTCCTTGCGCGGCCAACCGTTGCCGCCGCGGAGATGGATGCATTACAGGTAGACGGATGGGACTCCCGCGTAAGCCGTGTCGCCCCGTGGAAAGATCTGCCACGCCGTGTTTTTGACGCCCCGCTGGATATGTTTTCCggg                                                                                                                                 >  1:98415/1‑251 (MQ=255)                                                            tttCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGc                                                                                                                                         <  1:108677/236‑1 (MQ=255)                                                                                     gTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCTACGTCCTGGCGCAGCAAACCGGTGCCGCCGCGGCGAGGGGTGCAGGACAGGTAGACTAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGaa                                                                                                                                                                                  <  1:202812/170‑1 (MQ=255)                                                                                                                            aggcTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGAc                                                          <  1:62394/251‑1 (MQ=255)                                                                                                                                 ggAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGccc                                                     >  1:115433/1‑251 (MQ=255)                                                                                                                                 ggAAAACCTGACACCACCGCATAAGCCTCGCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGAAGCATTACAGGTAGACGAATGGTACTCCCGCGTACGGCGTGTCGCACCGGGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCga                                                       <  1:118731/249‑2 (MQ=255)                                                                                                                                          gACACCACCGCATAATCCT‑TCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTGCCGGGCGCTATCCCACTGGCGATTTTTCGGCTTTTCGCCCCAAACTGGTCTTGAACAGGGAGACCGGGATATTGACCGCCCGACCTGgcn                                            <  1:40651/250‑2 (MQ=255)                                                                                                                                                                                    cGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGTCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCTGGCGCTCAATTGTTGGCGGTTACCAGCGGCGGCACCATTCCGGAtn  <  1:135010/251‑2 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                             GGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCTCGACGTCCTTGCGCAGCAAACCGTTGCCGCCGCGGCGATGGATGCATTACAGGTAGACGAATGGTACTCCCGCGTACGCCGTGCCGCACCGTGGAAAGATCTGCCAAGACGTGTTTTTGACGCCACGCTGGATATGCTTTCCGGGCGCTATCCCTCTGGCGATTTTTCTGCTTTTCGCCCCAAACTGGTCTGGAACAGGGAGACCGGGATATTGACCGCCCGACCTGGCGCTCAATTGTTGGCGGTTACCAGCGGCGGCACCATTCCGGATC  >  NC_000913/1730311‑1730740

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