Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,539	T→C	S414S (TCA→TCG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,539	0	T	C	96.7%	86.9 / ‑4.2	30	S414S (TCA→TCG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (14/15);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

CAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACC‑GAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCAT  >  NC_000913/64305‑64736                                                                                                                                                                                                                                           |                                                                                                                                                                                                       cAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCt                                                                                                                                                                                         <  1:341314/250‑1 (MQ=255)  aGCGCCTGCGCCAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTATCGCGAGAACCAGGGATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGACTTCCACCAGCCCGACGGGATCACTCAGAAAGGTGCGGTTGCTCGACGGGTACAGGCTTTTATAGTCCAGAACCAGTACCGAATCATAAAGCCCTgg                                                                                                                                                                                       >  2:433952/1‑251 (MQ=255)     gcCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCg                                                                                                                                                                                    >  1:276758/1‑251 (MQ=255)         gCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTGTTCTCGCGAGAACCATGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTGTATAGTCCAGCTCCAGTACCGAATCATAAAGCCCTGGCCGTGaa                                                                                                                                                                                <  2:160475/251‑1 (MQ=255)           gACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATAAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGAACCGAATCATAAAGCACTGGCCGTGAATc                                                                                                                                                                              >  1:426816/1‑251 (MQ=255)              aGCGGATTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCAAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCac                                                                                                                                                                           >  2:25667/1‑250 (MQ=255)                        tACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAgcc                                                                                                                                                                  <  1:72954/250‑1 (MQ=255)                                 gTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGACTCATAAAGCCCTGGCCGTGAATCCATCCCGTAGCCGCCAGGGCTg                                                                                                                                                        >  2:341097/1‑251 (MQ=255)                                 gTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTg                                                                                                                                                        <  1:398662/251‑1 (MQ=255)                                             cgcgCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCggcggc                                                                                                                                            <  1:1429/251‑1 (MQ=255)                                               cgcCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGCTCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGACAGGTGCGGATGCTCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCCCGTAGCCGCCCGGGCCGGCGTGTggcggccc                                                                                                                                          >  2:129183/1‑249 (MQ=255)                                                         aGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCg                                                                                                                                   >  1:376570/1‑248 (MQ=255)                                                               tAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGTTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAgg                                                                                                                           <  1:128054/250‑1 (MQ=255)                                                                          cccGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCAGCGAGAAACCCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAAcc                                                                                                               <  2:146533/250‑1 (MQ=255)                                                                             tgcaggcaATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAgc                                                                                                            >  1:68618/2‑251 (MQ=255)                                                                                     aTGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGGGCTCTGGATCAGGATGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCa                                                                                                    <  2:340605/251‑1 (MQ=255)                                                                                         ttttCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCAtac                                                                                                 <  1:276107/250‑3 (MQ=255)                                                                                                  gagaACCAGGCATCGAGAAAACCTTCGGTACTGGGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAAGCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATa                                                                                       <  2:218783/251‑1 (MQ=255)                                                                                                     aaCCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCgag                                                                                                                                 >  1:195651/1‑206 (MQ=255)                                                                                                     aaCCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCgag                                                                                                                                 <  2:195651/206‑1 (MQ=255)                                                                                                      aCCAGGCATCGAGAAAACCTTGGGTACTGGGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGGGCGGATTATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAgaga                                                                                   <  2:399455/251‑1 (MQ=255)                                                                                                                   aaaaCCTTCGGTACTGTGATCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACgccg                                                                      >  1:198567/1‑251 (MQ=255)                                                                                                                                cTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGc                                                                                                                                                      >  1:133671/1‑158 (MQ=255)                                                                                                                                cTGTGCTCTGGATCAGGCGGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGc                                                                                                                                                      <  2:133671/158‑1 (MQ=255)                                                                                                                                           aacAGGCTGCGCCATGCCGTCCACCAGCCCGGCGGGATCAATCAGAAAGGGGCGGATGAACGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACC‑GATCCCCCGTGGCGGTcta                                             <  2:171170/249‑3 (MQ=255)                                                                                                                                                gCTGCGCAAGGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGAGGAGCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAAGCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACC‑GATCCCCCGTGGCGTTCTACCGGc                                        <  2:429252/251‑1 (MQ=255)                                                                                                                                                     gcCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCCGCACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGACAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGCTGACCAAACGCCGCCACC‑GATCCCCCGTGGCGGTCTACCGGCAGGcc                                   >  1:59492/1‑251 (MQ=255)                                                                                                                                                       cATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGCTGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGACGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTTGGCGCGACATAACCAGCGCGGTGCATCAGAGGAAAATAGAGATGACCCAACGCC‑CAACCCGATCCCCAGTGTCGGGCTACCGGCAGGCCgg                                 >  2:195725/1‑250 (MQ=255)                                                                                                                                                               ccaccaGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCGTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACC‑GATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGtc                         <  2:167099/251‑2 (MQ=255)                                                                                                                                                                                      aGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGAGGTGAATACGAGGAAAATAGAGATGACCAAACGCCGCCACC‑GATCCCCCGGGGCGGGCTACCGGCAGGCCGGTCACCGTCGGCCGTTCGAGGAAAAAAGGCAt  >  2:41186/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                       CAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACC‑GAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCAT  >  NC_000913/64305‑64736

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