Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,725,622	T→G	29.4%	V640G (GTG→GGG)	polB →	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,725,622	0	T	G	29.4%	24.4 / 4.3	32	V640G (GTG→GGG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base G (2/8);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGAGGTGCACTACCGGCGCTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGCCCTACCGTGATTATGTCCGCGAGTATGTACGGCGCACGTTGG  >  ACB143/2725489‑2725735                                                                                                                                      |                                                                                                                  ttCGAGGTGCACTACCGGCGCTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTgc                                                                                                                 >  4:287066/1‑136 (MQ=255)   cGAGGTGCACTACCGGCGCTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTgcgg                                                                                                               >  7:4683/1‑135 (MQ=255)     aGGTGCACTACCGGCGCTTCCTGATGCCCACCAACCGCGCCACGGACGAGGGCAGCAAGAAGCGCTATGCCGGCCGGGTGCCGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCAc                                                                                                             <  7:179104/136‑1 (MQ=255)     aGGTGCACTACCGGCGCTTCCTGATCCCCACCAACCGCGGCACCGACGAGGGCAGCAAGACGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCAc                                                                                                             <  6:149897/136‑1 (MQ=255)          cACTACCGGCGCTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGAAGAAAGAAGCGCTATGCCGGCGTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACt                                                                                                        <  3:210579/136‑1 (MQ=255)                   cgcATCCTGATGCCCACCACCCGCGGCCCCGACGAGGGCAGCAAGACGCGCTAGGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGGGCGCACCGACTGGTCACCGc                                                                                               <  1:112092/136‑1 (MQ=255)                     cTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCACACCGACTGGTCACCGCTg                                                                                             >  1:224053/1‑136 (MQ=255)                     cTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGGGCGCACCCACGGGTCACCGcgg                                                                                             >  1:252595/1‑134 (MQ=255)                             tGCCCACCATCCGCGGCACCGCCGCGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGGGCTGAAGAAATGGTTTACAAGGGCCTGGAGGCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCa                                                                                     <  6:11445/136‑1 (MQ=255)                                  cccaTCCGCGTCACCGCCGCGGGCAGCACGACGCGCTATGCCGGCCTGGGGCAGCGGGCCGACGGTGCTAAAGAATTGGTTTCAACGGGCCTGGAGGGGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTcc                                                                                <  4:174141/135‑1 (MQ=255)                                     agccGCGCCCCCGCCGCGGGCAGCAAGACGCGCTATGCCGGCCGGGGGCCGCGGGCCGACGGTGCTGAAGAAATGGTTTACACGGGCCGGGAGGCGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCcagc                                                                             <  3:12909/134‑1 (MQ=255)                                     acccGCGGCCCCGCCGCGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCCGCGGGCCGACGGTGCTGAAGAAATGGTTAACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCcagc                                                                             <  7:277057/134‑1 (MQ=255)                                         gcgGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCCGCAgg                                                                          >  2:266459/1‑135 (MQ=255)                                           ggCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCGCCGCTGGCCCGGCCGTTCCCGCAGGAAc                                                                       >  1:327159/1‑136 (MQ=255)                                            gCACCGCCGCGGGAAGAAAGAAGCGCTATGCCGGCGTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGCGCCTGGAGTCGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAAct                                                                      <  1:304586/136‑1 (MQ=255)                                              aCCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAActct                                                                    <  1:222250/136‑1 (MQ=255)                                               ccGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCCCCGCCTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAActcg                                                                    >  4:12909/1‑134 (MQ=255)                                                             aagaagCGCTATGCCGGCCTGGGGCCGCGGGCCGACGGGGCTGAAGAACTGGTTTACAAGGGCCGGGAGGCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTcc                                                     <  3:287066/136‑1 (MQ=255)                                                               gaagGGTTATGCGGGCCTGGTGCTGCGGGCCGAGGGTGCTTAAGAATTGGTTTACAAGGGCCTGGAGTCGGTGGGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTATGTTCCGc                                                   <  3:104025/136‑1 (MQ=255)                                                               gaagCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCCCCGCCTGGTGAACGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCgtgggtggcgc                                                   >  1:112940/1‑131 (MQ=255)                                                               gaagCGCCTGGCCGCCCGGGGGCCGCGGGCCGCCGGGGCTGAAGAAATGGTATACCCGGGCGGGGAGGCGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGc                                                   <  3:196833/136‑1 (MQ=255)                                                                  gcgcTATGCGGGCGTGGTGCTGCGGGCCGACGGTGCTGAATAATTGGTTTACACGGGCCGGGAGTGGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGc                                                <  5:137620/136‑1 (MQ=255)                                                                  gcgcTATGCCGGCCTGGTGCCGCGGGCCGACGGTGCTGAATAAATGTTTTCCAAGGGCCTGGCGGCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGc                                                <  7:41954/136‑1 (MQ=255)                                                                        tGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGcc                                           <  7:123667/135‑1 (MQ=255)                                                                                     cAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTGCGCAGCGCGCCCTAACGTGATTATg                             >  1:322661/1‑136 (MQ=255)                                                                                     cAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGCCCTACCGTGATTATg                             >  4:72537/1‑136 (MQ=255)                                                                                      cgCGGGCCGACGGTGCTGAAGAAATGGTTTACACGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGCCCTACCGTGATTATGt                            <  6:283421/135‑1 (MQ=255)                                                                                      cgCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGGCGGGGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGCCCTACCGTGATTATGt                            <  1:108784/135‑1 (MQ=255)                                                                                      aGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGCCCTACCGTGATTATGt                            <  6:199407/136‑1 (MQ=255)                                                                                                cGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGGTCCGCCGCGAGCCCGACCGTGATTAAGTCAGCGAGTAt                  >  6:389654/1‑136 (MQ=255)                                                                                                cGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGGGCGCCCCGACAGGTCCCCCCGGGCCCCGCACTTGCCGCAGGAACTCTAAGGCCGTGGGTTGCGGAGCGCGCCCCAACCGGATTATTTCCCCGGGTAt                  >  4:367432/1‑136 (MQ=255)                                                                                                                gTTTACAAGGGCCTGGAGTCGGTGCGCAACGACTGGTCACCGCAGGCCCCGCAGTTGCAGCAGGAACTCTACCGCCGTGGGTTACGCAGCCCCCCCAACCGAGAGTAATTGAGCGAGTATATAACGCGCACGATgg  >  6:78412/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                  TTCGAGGTGCACTACCGGCGCTTCCTGATGCCCACCATCCGCGGCACCGACGAGGGCAGCAAGAAGCGCTATGCCGGCCTGGTGCAGCGGGCCGACGGTGCTGAAGAAATGGTTTACAAGGGCCTGGAGTCGGTGCGCACCGACTGGTCACCGCTGGCCCGGCAGTTCCAGCAGGAACTCTACGGCCGTGTGTTCCGCAGCGAGCCCTACCGTGATTATGTCCGCGAGTATGTACGGCGCACGTTGG  >  ACB143/2725489‑2725735

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