Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,806,446	G→A	66.7%	G26E (GGA→GAA)	pfkB →	6‑phosphofructokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,806,446	0	G	A	66.7%	40.1 / 35.1	42	G26E (GGA→GAA)	pfkB	6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base G (8/6);  new base A (13/15);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTGGTTTCACTGTTGGC  >  NC_000913/1806303‑1806590                                                                                                                                                |                                                                                                                                                 atatTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAaaaa                                                                                                                                              <  4:103674/148‑1 (MQ=255)              cttcACTTTCCGCTGATTCGGTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGGCACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCg                                                                                                                                 >  2:213486/1‑148 (MQ=255)               ttcACTTTCCGCTGATTCGGTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCa                                                                                                                               >  1:173462/1‑149 (MQ=255)                        cGCTGATTCGGTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTc                                                                                                                      <  5:46176/149‑1 (MQ=255)                        cGCTGATTCGGTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTc                                                                                                                      <  7:252350/149‑1 (MQ=255)                          cTGATTCGGTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGa                                                                                                                    <  4:91019/149‑1 (MQ=255)                                 ggTGCTAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCggg                                                                                                             <  2:173462/149‑1 (MQ=255)                                     ctaGACTGACATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGggcg                                                                                                         <  8:64972/147‑1 (MQ=255)                                             aaaTCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATCGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCa                                                                                                   <  3:165197/147‑1 (MQ=255)                                                tCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAAc                                                                                              <  8:40823/149‑1 (MQ=255)                                                   gCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTTGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTc                                                                                           <  3:35028/149‑1 (MQ=255)                                                     cTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGc                                                                                         >  1:292032/1‑149 (MQ=255)                                                           gaggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAAATACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGtt                                                                                                                       >  8:154068/1‑113 (MQ=255)                                                              gaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCAt                                                                                <  4:149822/149‑1 (MQ=255)                                                               aaaTGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCAtt                                                                               <  3:80808/149‑1 (MQ=255)                                                                  tgatgGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCAt                                                                                <  6:155858/145‑1 (MQ=255)                                                                  tgatgGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAAGTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCAt                                                                                >  5:155858/1‑145 (MQ=255)                                                                    atgGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACACCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGcc                                                                            >  5:100537/1‑147 (MQ=255)                                                                     tgGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGCACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCAt                                                                         >  8:319894/1‑149 (MQ=255)                                                                        tACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCtt                                                                      <  6:4224/149‑1 (MQ=255)                                                                           gTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGa                                                                   <  3:290364/149‑1 (MQ=255)                                                                                tataCGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCCCCGGTGTTCGAACCCGGGGGCGGCGGCATCAGCGTCGCCCGCGCCATTGCCCATCTTGGAGGCAg                                                              >  2:214536/1‑149 (MQ=255)                                                                                tataCGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAg                                                              >  8:298484/1‑149 (MQ=255)                                                                                 ataCGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCAcc                                                                                                                             >  2:237217/1‑85 (MQ=255)                                                                                 ataCGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCAcc                                                                                                                             <  1:237217/85‑1 (MQ=255)                                                                                     gTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCca                                                         >  8:151314/1‑149 (MQ=255)                                                                                     gTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCca                                                         >  5:70293/1‑149 (MQ=255)                                                                                      ttGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCca                                                         >  6:173186/1‑148 (MQ=255)                                                                                           acTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGa                                                   >  4:222031/1‑149 (MQ=255)                                                                                             ttGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCGCATCTTGGAGGCAGTGCCACAGCGATc                                                 <  4:180185/149‑1 (MQ=255)                                                                                                cgcCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTc                                              <  1:163188/149‑1 (MQ=255)                                                                                                          cGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTg                                    >  4:279473/1‑149 (MQ=255)                                                                                                                     caaTTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGc                         <  2:292032/149‑1 (MQ=255)                                                                                                                          aCCCCGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAac                     <  7:151314/148‑1 (MQ=255)                                                                                                                           ccccGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGTCATCAACGTCGACCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGTGTTGCGCGACCGGCGAacac                   >  6:113658/1‑149 (MQ=255)                                                                                                                            cccGCAAATTTATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACAcc                  >  1:158972/1‑149 (MQ=255)                                                                                                                              cGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCCGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTg                >  6:168718/1‑149 (MQ=255)                                                                                                                                    tttATCCCGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTGGTTTCa          >  1:193855/1‑149 (MQ=255)                                                                                                                                         cccGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAgg                                                                 <  1:120646/89‑1 (MQ=255)                                                                                                                                         cccGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAgg                                                                 >  2:120646/1‑89 (MQ=255)                                                                                                                                         cccGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTGGTTTCACTGtt     >  3:220164/1‑149 (MQ=255)                                                                                                                                          ccGAAGAAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTGGTTTCACTGTTg    <  3:82520/149‑1 (MQ=255)                                                                                                                                            gaagaaAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTGGTTTCACTGTTGGc  <  6:70293/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 ATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCACAGCGATCTTCCCGGCGGGTGGCGCGACCGGCGAACACCTGGTTTCACTGTTGGC  >  NC_000913/1806303‑1806590

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