Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,588	G→A	L65L (CTG→TTG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,588	0	G	A	94.7%	53.5 / ‑3.5	19	L65L (CTG→TTG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGT‑CGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATA  >  NC_000913/65346‑65815                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      gTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGc                                                                                                                                                                                                                               <  1:87139/251‑1 (MQ=255)                cAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATACTTTCACCCCAGCGGCGTCCAGCCAAAGCCTTGTTCAACCTGc                                                                                                                                                                                                                 >  2:197159/1‑250 (MQ=255)                       gCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGGGTCAAGCGAAAGCCGTGTTCACCCTGCCAAATATg                                                                                                                                                                                                         >  2:227313/1‑251 (MQ=255)                        cATTAACGATAGTGCTATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATACACCGCTCAATCAGATAGCGTTCTGGCGGACGCCCATCGGCCTCGTAGACGGTAACGCAACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGACTGGCGGTGAAAATCCTTTAACGCCAGCGGGGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGc                                                                                                                                                                                                        >  1:50813/1‑251 (MQ=255)                                              cATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCCAAATATGCTGAGCGCGGGGAACCTGATcgg                                                                                                                                                                                  >  2:53314/1‑251 (MQ=255)                                               atatCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATcggc                                                                                                                                                                                 <  1:202592/251‑1 (MQ=255)                                                              cacaCCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTGCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCAc                                                                                                                                                                 <  1:106807/251‑1 (MQ=255)                                                                  ccGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGAc                                                                                                                                                             >  1:179752/1‑251 (MQ=255)                                                                       gaggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGCAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTg                                                                                                                                                        <  1:135264/251‑1 (MQ=255)                                                                                                        ggCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGtt                                                                                                                          <  2:210212/248‑1 (MQ=255)                                                                                                                           cGTAGACGGTTACGCCACCATCACGCAGGCGCTTTTCGTAATTCATCAAGTGGCGATGGGCGCGACAGTAAAGGCCATCCACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAAGATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGga                                                                                                     <  2:217024/250‑1 (MQ=255)                                                                                                                                                         gcTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAg                                                                      <  1:122923/251‑1 (MQ=255)                                                                                                                                                                gtaatcaatcaATTGGCGTTGGGCGCGACAGTAAAGGCAATACACCGGCTGGCGGTGAAAATCCTTTAGCGCCAGCGGTGTCAAGGGAAAGCCTTGTTCACCCGGCAAAATATGCTGAGCGCGGGGACCCTGATCGGCGGGAATAAACG‑CCACGGCCTCTTGCGGGGCAAGCGTAACCTGCAACGGCCCGTTGTCCGGCTCCAGCCAGAGGGAGGCTTCTGTCCCTTGCGGTGTGTCCCTCCAGTGT‑Cggg                                                              <  2:141003/244‑1 (MQ=255)                                                                                                                                                                               gCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAAAG‑CCCCGGACTCTTTCGGTGCAAGCGTAACCTGCAACCGCCCGTTGTCCGTCGCCCGCCAGAAGGAGACTTCCGTCCCTTGCGGGGGGTCCCGCCCGGGT‑CGGGGTAAAAAAAAAccct                                               >  1:16594/1‑249 (MQ=255)                                                                                                                                                                                                aaaGGCCATACCCCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGT‑CGGGTTAAGATAAAACCTGCCTGCGCCACGCTGaaa                              <  1:140901/251‑1 (MQ=255)                                                                                                                                                                                                                  ggCGGTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGACAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGACCA‑GGACTCTTGCGGGGGAAACGTAACCTGCAACGGCCCCTTGTCCGTCGCCAGCCCGAAGGAGACTTCTGTCCCCTTCGGGGGGTCCCCCCAGTGTTCGGGTAAAGA‑AAAAACTCCCTGCGCCCCGCGGCAAAACCACCaaaaaaaca              >  2:159547/1‑246 (MQ=255)                                                                                                                                                                                                                     ggTGAAAATCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGGGTCCCGCCAGTGT‑CGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCtt         <  1:240557/251‑1 (MQ=255)                                                                                                                                                                                                                          aaaTCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTTCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGT‑CGGGTTAAGATAAAACCCTCCCGCGCCACGCTGAAAATCCATCAAAAAACCACGCCTGAGt     >  1:120500/1‑250 (MQ=255)                                                                                                                                                                                                                            aTCCTTTAACGCCAGCGGTGTCAAGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGT‑CGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtata  <  1:198226/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACG‑CCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGT‑CGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATA  >  NC_000913/65346‑65815

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