Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	714,866	G→A	16.7%	A437T (GCA→ACA)	pgm →	phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	714,866	0	G	A	16.7%	56.6 / 4.9	24	A437T (GCA→ACA)	pgm	phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base A (2/2);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCACCCTGGCAGGTGACCCGATCACCGCGCGCCTGACTG  >  NC_000913/714728‑715007                                                                                                                                           |                                                                                                                                              gagagTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCttt                                                                                                                                    >  2:142978/1‑150 (MQ=255)         aGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGcc                                                                                                                            <  1:175971/150‑1 (MQ=255)            ggCTTCCTTCCTGCGTTTCGACGGCACGGCGTGGTCCAGCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACATGTAAGAACCCGCAGGAACACTACTACGAACTGGCAAAACGCTTTGGCGCGCCGAg                                                                                                                         <  1:155695/150‑1 (MQ=255)                   ttcctGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGACAAAACGCTTTGGTGCGCCGAGCTACAAc                                                                                                                  >  2:66863/1‑150 (MQ=255)                   ttcctGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGACAAAACGCTTTGGTGCGCCGAGCTACAAc                                                                                                                  >  2:66857/1‑150 (MQ=255)                             tCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACCAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTgcagg                                                                                                        <  2:169125/150‑1 (MQ=255)                              cGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTgcaggc                                                                                                       >  2:19833/1‑150 (MQ=255)                                aCGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTgcaggcag                                                                                                     >  2:15130/1‑150 (MQ=255)                                 cGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGc                                                                                                    >  1:200994/1‑150 (MQ=255)                                    cACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGc                                                                                                 <  1:168989/150‑1 (MQ=255)                                         cGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACtt                                                                                            <  1:136028/150‑1 (MQ=255)                                                     aCAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAg                                                                                <  1:19833/150‑1 (MQ=255)                                                        aaGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAgcgg                                                                             <  2:200994/150‑1 (MQ=255)                                                                  catcatGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTaa                                                                   >  1:122691/1‑150 (MQ=255)                                                                     catGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGTAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCt                                                                <  2:150240/150‑1 (MQ=255)                                                                             tgctgGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTGTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCgg                                                        >  1:11388/1‑150 (MQ=255)                                                                                    ggcggAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGt                                                 <  2:122691/150‑1 (MQ=255)                                                                                     gcggAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTg                                                >  1:122614/1‑150 (MQ=255)                                                                                      cggAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGa                                               <  2:119067/150‑1 (MQ=255)                                                                                               ccGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCa                                      >  1:181584/1‑150 (MQ=255)                                                                                                                   ccGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCACCCTGGCAGGTGACCCGATc                  >  2:174977/1‑150 (MQ=255)                                                                                                                                 cAACGAACTGACAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCTGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCACCCTGGCAGGTGACCCGATCACCGCGCGCCTGAc    <  1:66863/150‑1 (MQ=255)                                                                                                                                 cAACGAACTGACAAAACGCTGTGGTGCGCCGAGCTACAACCGTTTGCCGGCAGCTGCGACTTCCGCACAAAACGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCACCCTGGCAGGTGACCCGATCACCGCGCGCCTGAc    <  1:66857/150‑1 (MQ=255)                                                                                                                                   aCGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCACCCTGGCAGGTGACCCGATCACCGCGCGCCTGActg  <  2:240538/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                              GAGAGTGCAGGGGCTTCCTTCCTGCGTTTCGACGGCACGCCGTGGTCCACCGACAAAGACGGCATCATCATGTGTCTGCTGGCGGCGGAAATCACCGCTGTCACCGGTAAGAACCCGCAGGAACACTACAACGAACTGGCAAAACGCTTTGGTGCGCCGAGCTACAACCGTTTGCAGGCAGCTGCGACTTCCGCACAAAAAGCGGCGCTGTCTAAGCTGTCTCCGGAAATGGTGAGCGCCAGCACCCTGGCAGGTGACCCGATCACCGCGCGCCTGACTG  >  NC_000913/714728‑715007

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