Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,594	A→C	F63V (TTT→GTT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,594	0	A	C	93.5%	90.6 / NA	31	F63V (TTT→GTT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (11/18);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGC  >  NC_000913/65352‑65836                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   aTGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGtt                                                                                                                                                                                                                                            <  2:270358/251‑1 (MQ=255)              cGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAat                                                                                                                                                                                                                               <  1:151712/251‑1 (MQ=255)                     tAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGACTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGGGAAAATCCTTTAAAGCCAGCGGTGTCAGGCGAAAGCCTTGGTCACCCTGCCAAATATGCTg                                                                                                                                                                                                                         >  2:36312/1‑250 (MQ=255)                            agcgCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGAGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCTCATCGGCCTCGTAGACGGTAACGCCACCGTCACGCGGGCGCTTTTCGTAATTCATTAATTGTCGATGGGCGCGGCAGTAAAGGCCATACACCGGCTGGCGGGGAAAATCCTTTAACCCCAGAGGGTTCCGGCGAACGCCCTGTTAATCCTGTCAAATAATCTGAGcgcgtg                                                                                                                                                                                                                 >  1:136158/4‑249 (MQ=255)                                     gtgCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGa                                                                                                                                                                                                         <  1:78878/250‑1 (MQ=255)                                             cACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCgggg                                                                                                                                                                                                >  1:183205/1‑250 (MQ=255)                                                  tCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGGATaa                                                                                                                                                                                            >  1:179638/1‑250 (MQ=255)                                                            ccGGTGCGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGCAGGGAACGCCACATTCACGAAGGCGCTTTTCGTAATTCATCCATTGGCGATGGGCGCGACAGGAAAAGCCATACAACGGATGGCGGTGAAAATCCTTTTACGCAAGCGGGGTCAGGCGGAAGGCTTGTTTACCCCGGAAAAAATGCTGAGCGCGGGGAAACCGAGCGGCgggggtaaaacgccccgg                                                                                                                                                                                  >  2:53194/1‑235 (MQ=255)                                                                 gaggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTg                                                                                                                                                                            <  2:218773/251‑1 (MQ=255)                                                                   ggtgaTAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCg                                                                                                                                                                          <  1:4000/251‑1 (MQ=255)                                                                         aaaCCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAAAGCCAccggatctttgggtgca                                                                                                                                                                     >  1:112806/1‑234 (MQ=255)                                                                                   aTCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGGCTAAAACCCACGGACTCTTGCGGGGCAAGCGTAAccat                                                                                                                                                          >  1:261408/1‑249 (MQ=255)                                                                                     cAGATAGCGTTCTGGCGGGCGCGCATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGc                                                                                                                                                         <  2:84616/250‑1 (MQ=255)                                                                                      aGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCaa                                                                                                                                                       <  2:101386/251‑1 (MQ=255)                                                                                                                        aGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAACGCCATACACCGGCTGGCGGCGAAAATCCTTTAACGCCAGCGGCGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACCCCACGGGCTCTTTCGGCGGTAGCGTAACCCGCAACGGCCCGATTTCCCTCgccagcccacacgggaat                                                                                                                     >  1:254997/1‑240 (MQ=255)                                                                                                                           cGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCt                                                                                                                  <  2:179638/251‑1 (MQ=255)                                                                                                                           cGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCt                                                                                                                  <  2:92243/251‑1 (MQ=255)                                                                                                                               aaCGCCCCCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTc                                                                                                               <  1:208897/250‑1 (MQ=255)                                                                                                                                   ccACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAAAGCCACGGACTCTTGCGGGGCAAGGGGAACCTGCAAAGGGCCGTTGGTCGTCGCCCGCCCGAAGGAGACTTTTGTACCTTg                                                                                                          >  1:219703/1‑251 (MQ=255)                                                                                                                                            cGCAGGCGCGGTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAACGGCCAGACACCGGCTGGCGGTGAAAAGCCTTTAACGCCAGCGGTGTCAGGCGAACGACTTTTTCACCCTGCAAAATATGCTGAGCGCGGGGAAACTGATCGGCGGGAATAAACGCCTCGTACTCTTGCTGTGCAAGAGTAACCTGCATCGGCCCGTTGTCCGTCGCCGTCCAGAAGTAGACTTTTGTTCCTGGCGGGTTGTc                                                                                                 <  2:249418/251‑1 (MQ=255)                                                                                                                                              cAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTccc                                                                                               <  2:58604/251‑1 (MQ=255)                                                                                                                                                ggCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAAAGCCACGGACTCTTGCGGGGCAAGCGTAACCTGCAACGGCCCGGTGTCCGTCGCCCACCCGAAGGGAAATTCTGTCCCTTGCGGGGGGGCCCGc                                                                                             >  1:190167/1‑251 (MQ=255)                                                                                                                                                                            aTGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAAcctgcct                                                                 <  1:207397/251‑1 (MQ=255)                                                                                                                                                                                     aCAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCCCGGGCTCTTTCGGTGCAAGCGTAACCTGCCACGGCCCGGTGTCCGCCGCCAGCCAGAAGGAGCCTTTTGTCCCTTGCGGGGGGTCCCGCCAGTGGCGGGttgaaaaaaaaaccgccagcggcccgc                                                        >  1:259259/1‑226 (MQ=255)                                                                                                                                                                                      aagaaaagaccatccgccGGCTGGCGGTGAAAATCCTTTACCGCCAGCGGTGTCAGGCCAACGCCTTGTGCCCCCTGCAAAATATGGGGAGCGCGGGGACCCTGATCGGCGGGAATAAGCGCCACGGCCTCTTGCGGTGCAAGCTTAACCGGCACCGGCCCGTTGGCCGTCGCCACCCAGAAGGAGACTTCTGTCCCTTGCGGGGTTTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCt                                                       <  2:136158/235‑1 (MQ=255)                                                                                                                                                                                              gCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCGGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTgtaattcc                                               <  1:129934/251‑8 (MQ=255)                                                                                                                                                                                                                  gAAAATCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGGGTCCCGCCAGTGTCGGGGTAAGATAAAACCCGCCTGCGCCACCCTGTAATTCCCTCCAAAAAACAGGCTTGat                           >  1:38851/1‑250 (MQ=255)                                                                                                                                                                                                                        ccTTTAACGCCAGGGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTTCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGGCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCAAAAAACCAGGCTTGATTATAg                      <  1:166372/250‑1 (MQ=255)                                                                                                                                                                                                                                gCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTc             <  1:204193/251‑1 (MQ=255)                                                                                                                                                                                                                                       gtgtCAGGCGAACGCCTTGTTCCCCCTGCAAAATATGCGGGGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTACCCTGCAACGGCCCGTTGTCCGTCCCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGTAATTCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGAt      <  2:88722/251‑1 (MQ=255)                                                                                                                                                                                                                                           cAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACATGATCGGCGGGAATAAACACCACGGACTCTTGCGGTGCAAGCGTAAACTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAAAATTCTGTCCCTTGCGGGGGGTCCCGCCAGTGTCGGGGTAAGATAAAACCTGCCTGCGCCCCGCTGTAATTTCCTCAAAAAACCAAGGTTGATTTTATCCTGGTTTCGTTTGAtttgg  >  1:159591/1‑248 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   ATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGTGAAAATCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGC  >  NC_000913/65352‑65836

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