Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,316	T→A	I155I (ATA→ATT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,316	0	T	A	96.0%	71.8 / NA	25	I155I (ATA→ATT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (11/13);  minor base C (1/0);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGT  >  NC_000913/65081‑65560                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       gcagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATAAAT‑AGAAACCCACTTg                                                                                                                                                                                                                                        >  1:145142/1‑251 (MQ=255)         cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGGTGGCAAACCAGGCGTCGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGCCGTATTAAAGTTCAAAATCGAGCGCGCAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTACAGGCCGATGCAGTACAGCTCACCGTGGCGGATGGTTGCACTATCCATCACAAACCC‑CATGAGcggcgg                                                                                                                                                                                                                                >  2:285818/1‑251 (MQ=255)                       caccacGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCCATATCAAT‑AGAAACCCACTTGAGCGGCGGGCGaaaggtggg                                                                                                                                                                                                                    >  2:17236/1‑241 (MQ=255)                               tCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑ttt                                                                                                                                                                                                          <  1:100665/250‑1 (MQ=255)                                          atCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCCATATCAAT‑AGAAAACCACTTGAGCGGCGGACGAGAGTCGGGGTGCGG‑GTTCAGAAggggctt                                                                                                                                                                                              >  2:128281/1‑247 (MQ=255)                                            cacaTCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAATCGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAAT‑AGAAACACACTTGAGCAGCGGGCGGTAGATGGGATGCGGGTTTC‑GACGGGGATTa                                                                                                                                                                                             >  2:221791/1‑250 (MQ=255)                                                              ggCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGc                                                                                                                                                                           <  2:51423/250‑1 (MQ=255)                                                              ggCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCa                                                                                                                                                                          <  1:169402/251‑1 (MQ=255)                                                                  aaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATc                                                                                                                                                                      <  1:84992/251‑1 (MQ=255)                                                                    ccAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCAc                                                                                                                                                                    <  1:74289/251‑1 (MQ=255)                                                                     cAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑ATAAACCCACTTGAGCGGCGGACGCTTATCGGGATGCGG‑TTTCAGACGGGCATTAAAGATAGCGCCATTGTGCATATCAcc                                                                                                                                                                   >  1:6345/1‑251 (MQ=255)                                                                         cGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCCATATCAAT‑AAAAACCCACTTGAGCGGCGGGCGATAGTCGGGATGGGG‑GTTCAGACGGGCCTTAACGAAAACGCCATTGTGCataaaccccccc                                                                                                                                                               >  2:170992/1‑242 (MQ=255)                                                                                 ttttccagtagatGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGGTCAATATCAAT‑AGAAAACCACCTGAGCGGGCGACGGTGGTGGGGGTGCGG‑TGTTTGAAGGGCATTAAGCATAGAGCCATTGGTCCTAACACCCCCGACCCacac                                                                                                                                                       >  2:11949/13‑251 (MQ=255)                                                                                          agcTGCGGGCGGCTGGCGACGTATGCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtgaggt                                                                                                                                              <  2:13596/249‑1 (MQ=255)                                                                                                         gCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCAAT‑AGAAACCCCCTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCCTTTACGAGAGCGCCATTGTGCATATCACCCTCGACCCCCCCCCGTGAGGGGATAAACCGCTCCAt                                                                                                                               >  2:119944/1‑251 (MQ=255)                                                                                                            aCGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAg                                                                                                                            <  2:168592/251‑1 (MQ=255)                                                                                                                        gcaaaaTCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTg                                                                                                                 <  2:99527/247‑1 (MQ=255)                                                                                                                                                         ggCCCCAGCATATAAACGATGCGCTGCCCGCAGTCTTGCAGGGCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGGTCTTGCGGGCGCACATCGGCCTCGTAGACGGTAAc                                                                                  <  1:292930/248‑1 (MQ=255)                                                                                                                                                                gCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCCATTTCAAT‑AGAAACCCCCTTGAGCGGCGGACGCTAGTCGGGGTGCGG‑GTTTAGAAGGGCATTAACGATAGCGCCACTTTGCATATCACCCTCGACCCCCACCGGTGAGGGGATAAACCGCCCCATCAGATAGCGCTCCGGCGGGCGCCCAACGGccccctcagcgggaacacccccttcc                                                                        >  2:311485/1‑227 (MQ=255)                                                                                                                                                                        aCGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCCCGTAGACGGGAACGCCACCTTTACGCAGgcg                                                                >  1:103013/1‑251 (MQ=255)                                                                                                                                                                                    ccGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGGTGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCCCACCGGGGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACAAAGGCCCCGTAGACGGGAACGCCACCCTCACGCAGGCGCTTTTCGTAAtt                                                    >  1:10112/1‑251 (MQ=255)                                                                                                                                                                                        aGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATtcatc                                                <  1:63401/251‑1 (MQ=255)                                                                                                                                                                                                    cGATGCAGTACAGCTCACCGTGGCGGGTGGTTTAAATAGCAAT‑AGAAACCGACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATgg                                    <  2:238995/251‑1 (MQ=255)                                                                                                                                                                                                             aCAGCTCACCGTGGCGGGTGGTTTCAATATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACa                            <  1:127868/250‑1 (MQ=255)                                                                                                                                                                                                                                      aaTATCAAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGt  <  1:117190/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       GCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTAT‑AGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGG‑TTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCGGT  >  NC_000913/65081‑65560

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: ‑