Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,755	G→A	I342I (ATC→ATT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,755	0	G	A	95.7%	55.6 / ‑6.2	23	I342I (ATC→ATT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/11);  total (12/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 = 1.00e+00

TTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGT  >  NC_000913/64518‑64996                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   tttaaagtccggacccagtaccgAATCATAAAGCCCGGGCCGTGAATCCATCACGTGGCCGCCAGGGCTGGCGTGGGGCGGCAATTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGGTCGAGTAAAAATGGCATGATTTCAGTGTTATGGAAAATCTGCGTCAcc                                                                                                                                                                                                                                        <  2:220103/228‑1 (MQ=255)               ccagtaccgAATCATAAAGCCCTGGCCGTGAATCCATCAAGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCCGCGCGGTGCATACGAGGAAAATAGCGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCAACCGGCAGGCCGTTCACCGTCGATCGTTAGAGTAAAAATGGCATGATTTCAGTTTTATGGGAAATCTGAGTCACCAGCT‑CGCAATCTTTc                                                                                                                                                                                                                        >  1:324921/9‑251 (MQ=255)                              aaaGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGGGTAAAAATAGCCTGATTTCCGTTTTATGGAAAATCTTCGTCACCAGCT‑CGCAATCTTTTCGGTTTTAAGTCGcc                                                                                                                                                                                                         >  2:45512/1‑251 (MQ=255)                                aGCCCGGGCCGTGAATCCATCCCGTAGCCGCCAGGGCTGTCGTGTGGCGGCACTTCGCCGAGAGCAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAg                                                                                                                                                                                                       <  2:174818/251‑1 (MQ=255)                                         cGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGttt                                                                                                                                                                                              <  2:451883/251‑1 (MQ=255)                                           tGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCAACCGGCAGGCCGGTCACCGTCGGTCCTTCGGGTGAAAATGGCAGTATTTCAGTTTTATGGAAAATTTGCGTCACCATCC‑CGCTCACTTTCAGGCTATAAGtaccccacgccggtttat                                                                                                                                                                                            >  2:371396/1‑234 (MQ=255)                                                         aGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGGGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCCATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAAcgg                                                                                                                                                                              >  1:242779/1‑251 (MQ=255)                                                           ccgccAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGGCTACCGGCAGGGCGTTCACCGTCGCACGTTCGAGTAAAAATTGCATGATTTTAGTTTTATGGGAAATATTCGTCCACAGCTACG‑AAACTTTCAGGTTATCAGTCGCCCGCGCAGGTTTTTATTTGGCGAAAcggcg                                                                                                                                                                            >  2:248861/1‑251 (MQ=255)                                                                                                               cAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGAATGATTTCAGTTTTATGGCAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCACCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTAAATTTTGTCCAATCGGACCCCCGGGTTATCGAGAGATTTTTCTTTGccc                                                                                                                        >  1:421477/1‑250 (MQ=255)                                                                                                                              gaggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCg                                                                                                         <  2:276758/249‑1 (MQ=255)                                                                                                                                                                  cgccgTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTAGCGTAAAAATGGCATGATTTCAGTTTTGTGGAGGATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGGCCATTCGGTCCCACGGGGTATCGATTGATTTTCCTTCGCCTAATAGCTCCTGCGGGACAGGTTCCAGCGAGAaagaaggataactccagaa                                                                      >  2:43893/1‑237 (MQ=255)                                                                                                                                                                                  cGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGCCAATTTCGTCCCTTCGATCCCCCGGGTGTTCGATAAATTTTCCTTTGCCTAATAACTCCTGGGCGACGGGTTCCCGAGAGAAAGAAGAGAAAATCCAAAAcgagggttttagacccc                                                     >  2:246427/1‑236 (MQ=255)                                                                                                                                                                                                 cgtcgctcGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATaa                                      <  1:41186/244‑1 (MQ=255)                                                                                                                                                                                                            cGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCG‑CGCAATCTTTCAGGGTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCt                           <  2:68618/251‑1 (MQ=255)                                                                                                                                                                                                              aGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTTAGGTTATAAGGCGCCAGCGCATGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCttt                         <  1:268797/251‑1 (MQ=255)                                                                                                                                                                                                                  aaaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAAATCCTGAGCGACGGTTTCCAGCGCGAATGAAGAGAAATTCCAGAACGCGGATTTAAGCGCCTCGCAACCGTCACTAATTAGCCGCCCTTTaacc                     >  1:330391/1‑248 (MQ=255)                                                                                                                                                                                                                   aaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTc                                                            <  1:187497/211‑1 (MQ=255)                                                                                                                                                                                                                     aTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTgg                  <  1:341097/251‑1 (MQ=255)                                                                                                                                                                                                                      tGGCATGGTTTCAGATTTATGGCAAATCTGCGTCACCAGCT‑CGCAATATTTCAGTTTATAAGTCGCCAGCGAAGGTTTATCTTCGGCGAAACGGCGGTCAATTTAGTCCCTTCGATCCCACGGGTTCTCGATCGATTTTCCTTTGCCTAATTACACCTGTGCGCCGGTTTCCAGGGAGACTGTAGATAAAATCCAGGAGGCGGATTTCAGCGCCACGATACCGTTCATGATTACACGACCTTTATCCtagg                 >  2:373119/1‑248 (MQ=255)                                                                                                                                                                                                                       ggCATGATTTAAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAACCGGCGGTCACTTTCGTCCATTCGCTCCCACGGGGTATCGGAAGATTTTCCTTCGCCTAATAACTCCTGCGCGACGGTTTCCCGCGAGAATGAAGAGAAATTCCTGAACGCGGAGTTCATCGCATCGAGACCGTCGATAATTAACCGTCCtttcgcatgcgc                >  2:81881/1‑242 (MQ=255)                                                                                                                                                                                                                            gATTTCAGTTTTATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTAGATAATTAGCCGACCTTTAGCCTGGGCaaaa            <  1:45512/250‑1 (MQ=255)                                                                                                                                                                                                                                    ttttATGGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGATAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGTGCAAAAAAGACGcc    <  1:137962/250‑1 (MQ=255)                                                                                                                                                                                                                                     tttatgGAAAATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGCCGGTTTCCAGCGCGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAACAAGACGCCGt  >  1:124193/5‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   TTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCT‑CGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGT  >  NC_000913/64518‑64996

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