Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,071,836	A→G	27.4%	F263L (TTC→CTC)	pgk ←	phosphoglycerate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,071,836	0	A	G	27.4%	14.5 / 4.8	16	F263L (TTC→CTC)	pgk	phosphoglycerate kinase
Reads supporting (aligned to +/- strand):  ref base A (5/6);  new base G (2/2);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-02

CGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAGCGATACCACCACCAACAATCAGCTGGTCAGCGATTTTAGACAGGGAGTCCAGAACGGTCAGTTTGGTAGATACTTTAGAACCACCAACGATAGCCACCATCGGGCGAGCAGGTTCTTTCAGTGCTTTACCCAGCG  >  NC_000913/3071626‑3072105                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                              cGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGCACTCGGTTGCTACGCGACCATCAGACGGAACCGGGATGTTGAAGGTGGTCAGCAGACGTTTAGCt                                                                                                                                                                                                               >  1:132707/1‑275 (MQ=255)                           aTTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTTCTCTTCAGCTTTCACATCGTTAACAGAGTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAg                                                                                                                                                                                                                 <  1:121977/247‑1 (MQ=255)                                                      aaTTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGt                                                                                                                                                                                                  <  1:316735/235‑1 (MQ=255)                                                                   cccccgccccGTTCCACACAATGGTTTACGCATGCTTCAGGATGTCGGCCAGTTACTGAGCTGTAGAAGCACCGAGATCCAGGGTCTGCTCCCCAGATTTCCCCTCGTTAACAGATTTGAGGGTAGACGGTGCCAGTGCAGAGATCTCGGTTGCTCCGCGACCAGCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAAGCGTTCAGCTGCGTCAACCAGGTCAGCATCGTACAGGGATTAACCCACAGCGTGGCCTTGTGCCGCGATAAAGgtg                                                                                                                                             <  2:132706/267‑1 (MQ=255)                                                                           ccGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCCGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGACCATCAGACGGAACCGGGAGGTTGCAGGTGGTc                                                                                                                                                                                                                              >  1:349524/1‑186 (MQ=255)                                                                                                                    cAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACGCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTAAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCcg                                                                                                                                                        <  1:310145/215‑1 (MQ=255)                                                                                                                                 gAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAGCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTa                                                                                                                                                                              >  2:24354/1‑180 (MQ=255)                                                                                                                                 gAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAGCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTa                                                                                                                                                                              >  2:24368/1‑180 (MQ=255)                                                                                                                                 gAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATATCAGTGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTCCGCGAACATCAGACGGAACCGGGGTGTTGCAGGTGGTCAGCAGACGTTTTGCTTCGTCAACCCGGTc                                                                                                                                                                                                 >  1:46927/1‑161 (MQ=255)                                                                                                                                     cATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAGCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCca                                                                                                                                                                          <  1:24354/180‑1 (MQ=255)                                                                                                                                     cATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAGCTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCca                                                                                                                                                                          <  1:24368/180‑1 (MQ=255)                                                                                                                                                                   cgttCACATCGTTAACAGATTTCATGGTAGCCCGTGCAGTTTCAGAGAACTTGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATaa                                                                                                                                                  <  1:138167/172‑1 (MQ=255)                                                                                                                                                                                        tCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGa                                                                                                                                                                                                                        >  1:407399/1‑83 (MQ=255)                                                                                                                                                                                            ggTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGANATTGCAGGTGGTCAGCAGACNtt                                                                                                                                                                                                                    <  2:407398/83‑1 (MQ=255)                                                                                                                                                                                               aGCCGGTGCAGTTTCAGAGTACTCGGGTGCTACGCGAACATCGGGCGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGGCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAGCGATACCACCACCAACAATCAGCTGGTCAGCGATTATAGACAGGGAGTGCAGAACGGTCAGTTTGGGAGATACTATAGAACCACCAACGATAGCCACCATCGGGCGAGCAGGTTCTTTCa                 <  1:99398/275‑1 (MQ=255)                                                                                                                                                                                                               gagaACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTGGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGACTTGTGACGCGATAAAGGTGTTAGCGATACCACCACCAACAATCAGCTGGTCAGCGATTTTAGACATGGAGTCCAGAACGGTCAGTGAGGTAGATACTTTAGAACCACCAACGATACCCACCATCTGGCGAGCAGGTTCTTTCAGTTCTTTACCCAgcg  >  1:167519/1‑274 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                              CGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACCCACATCGTGGCCTTGTGCCGCGATAAAGGTGTTAGCGATACCACCACCAACAATCAGCTGGTCAGCGATTTTAGACAGGGAGTCCAGAACGGTCAGTTTGGTAGATACTTTAGAACCACCAACGATAGCCACCATCGGGCGAGCAGGTTCTTTCAGTGCTTTACCCAGCG  >  NC_000913/3071626‑3072105

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