Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	501,402	C→T	S426S (AGC→AGT)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,402	0	C	T	89.5%	44.7 / ‑2.9	19	S426S (AGC→AGT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (11/6);  total (13/6)
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.26e-01

GGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGG  >  NC_000913/501207‑501619                                                                                                                                                                                                    |                                                                                                                                                                                                                          ggCAGCGTTACTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTCCCAAACTCCAGCCAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGCGTGTAACTATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTATTACTGGGATCGTTAAGTTATCGTCGGt                                                                                                                                                                                    >  1:379275/1‑235 (MQ=255)                                      cNTCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGAt                                                                                                                                                                                                      >  1:393310/1‑180 (MQ=255)                                      cATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTgg                                                                                                                                                 >  1:477934/1‑233 (MQ=255)                                            ngCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGAGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGa                                                                                                                                                                                                       >  1:308378/2‑173 (MQ=255)                                            aGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCGAACCGTGTGAGCTAGCAGGTACGGAGCCAGCATTCTCCTCGTGTAACGCGCTGCTTGCCGGGTCGTGAATATAGTCTGGAAGTGTCTTACGGGCAGCTTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATg                                                                                                                          <  1:467518/250‑1 (MQ=255)                                              cAACGTACCAAACTCCAGCAAACATAAATTCCCCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCg                                                                                                                        >  1:418467/1‑250 (MQ=255)                                                      cAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTACATATGCTAACCGTGTGAGCTATCAGGAACTGCACCAGCATTCACTTCGTTTAACGCGCGGCTTGCCGGGGCGTGACGACAGTCTGGAAGCGTCTTACTGGGATCGTTAAGt                                                                                                                                                                                              >  1:458907/1‑172 (MQ=255)                                                                   nnnnnnaTNNNNNTGGNTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGa                                                                                                                                                                                                       <  2:308378/144‑1 (MQ=255)                                                                              ccTGGTTAACTTATTCATCGTTAGCGCAAGTGTGTAAATCTGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCg                                                                                                                                                                 >  1:419569/1‑177 (MQ=255)                                                                                             cATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTACCGAGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATACGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATACCTCCATCTTTTTCATATTTTACATcc                                                                        >  1:194951/1‑251 (MQ=255)                                                                                             cATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATcc                                                                        <  1:233935/251‑1 (MQ=255)                                                                                                             gtgtAAATATGCTAACCGTGTGAGCTATCAGGTACTGCACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTg                                                                                                                                                  >  1:219647/1‑161 (MQ=255)                                                                                                                                                     cAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGCAGAGTCTTACTGGGCTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCAAGCAGtt                >  1:283822/1‑251 (MQ=255)                                                                                                                                                        naTTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTAt                                                                                                                                                                                           >  1:305909/2‑77 (MQ=255)                                                                                                                                                        cATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTAt                                                                                                                                                                                           <  2:305909/77‑1 (MQ=255)                                                                                                                                                                  gTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTGTCATATTTTACATCCGGCAACCACCGTTTAcc                                                       <  1:177521/199‑1 (MQ=255)                                                                                                                                                                  gTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGAATCTGGCATTTGTCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTCTCATGCCTACATAtttt                                                                                       >  1:64417/1‑167 (MQ=255)                                                                                                                                                                    ttAACGCGCGGCTTGCGGGGGCTTGAAGACTTTCTGAAAGAGTCTTACTGGGATCGTTAAGTTACCGTCGGTTCGTAGGCCAGATAAGGCTTTCACGCCGCATCTGCCATTTGGCTCTCGAGGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTCCATATTTGACACCCGCCACCCCCCGTTGCCCCCGTCACCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg  <  1:325757/250‑1 (MQ=255)                                                                                                                                                                         gngcgGATTGCCGGAGCGTGAAGACAGTCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCg                                                                       >  1:395268/1‑176 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                                          GGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGG  >  NC_000913/501207‑501619

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