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	92.9%	73.3 / NA	28	I155I (ATA→ATT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (11/15);  minor base C (2/0);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTG‑CATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTC  >  NC_000913/65070‑65504                                                                                                                                                                                                                                                       |                                                                                                                                                                                               ggCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAg                                                                                                                                                                                            <  1:35623/251‑1 (MQ=255)    aTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCCATAGaa                                                                                                                                                                                          >  1:79202/1‑250 (MQ=255)     tGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCACTATCAATAGaaa                                                                                                                                                                                         >  2:78878/1‑250 (MQ=255)      gTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAAcc                                                                                                                                                                                       <  2:233250/251‑1 (MQ=255)       tttttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCCATATCAATAGAAAccc                                                                                                                                                                                      >  1:132190/1‑251 (MQ=255)         tttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCa                                                                                                                                                                                     <  2:207725/250‑1 (MQ=255)                caTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAgcg                                                                                                                                                                             <  2:79202/251‑1 (MQ=255)                   tCGCACATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTTGCAACCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCTTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAgcggc                                                                                                                                                                           <  1:259108/250‑1 (MQ=255)                       aGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAAACCACTTGAGcggcgggcg                                                                                                                                                                      >  1:114530/1‑248 (MQ=255)                              aCTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGGTTCCATATCAATATAAACCCACTTGAGCGGCGGACGATA‑GTCg                                                                                                                                                              >  1:88722/1‑251 (MQ=255)                                 gCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGa                                                                                                                                                           <  1:90344/251‑1 (MQ=255)                                     cacGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCTGGATGCg                                                                                                                                                       <  2:151641/251‑1 (MQ=255)                                               ccGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTTAATATCAATAGAAACCCACCTGAGCGGCGGACGGGA‑GTCGGGGGGCGGGTTCAGACg                                                                                                                                             >  2:207397/1‑251 (MQ=255)                                                       cacctCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTAGTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCTGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACTGGCATTa                                                                                                                                      <  2:204400/246‑1 (MQ=255)                                                                      gTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGGCGGTA‑GGCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATtgt                                                                                                                      >  1:187211/1‑251 (MQ=255)                                                                      gTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATtgt                                                                                                                      <  1:181734/251‑1 (MQ=255)                                                                        tGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGCCGGGCATTAAAGATAGCGCCATTGTG‑c                                                                                                                   >  1:270358/1‑251 (MQ=255)                                                                                 aGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCACTATCCATAGAAAACCACTTGAGCGGCGGGCGGTA‑GTCGGGATGCGGTTTCATCAGGGCATTAACGAGAGCGCCATTTTG‑CCTATCAccc                                                                                                          >  2:34942/1‑251 (MQ=255)                                                                                               tccagtagcTGCGGGCGGCTGGCGACGTATTCCAGTGCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTG‑CATATCACCCTCGACCCACACCg                                                                                             <  1:239723/242‑1 (MQ=255)                                                                                                 cagtagcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTG‑CATATCACCCTCGACCCACACCggtg                                                                                          <  1:88148/245‑1 (MQ=255)                                                                                                     agcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTTAATATCAATAGAAAACCACTTTAGCGGCGGACGATA‑GTCGGGGTGGGGGTTCCGACGGGGCTTAACGATAGCGCCCTT‑TGTCATATAACCCTCGACCCCCACCggggaggg                                                                                      >  2:52929/3‑250 (MQ=255)                                                                                                      gcTGCGGGCGGCTGGCGACGTATTCAAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTG‑CATATCACCCTCGACCCACACCggtgagggg                                                                                     >  1:58604/2‑249 (MQ=255)                                                                                                                                 gTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTG‑CATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCt                                                          <  2:41316/251‑1 (MQ=255)                                                                                                                                         tcgaGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGAAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GCCGGGATGCGGTTTCAGACGGGCATTAACGGTAGCGACATTGTG‑CATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCgggcg                                                  >  1:10328/4‑248 (MQ=255)                                                                                                                                                 ggaggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCAAAAGAAACCCACTTGAGCGGCGGACGAAAAGT‑GGGATGCGGTTTTAGACGGGCATTAAAGAAAGCGGCCTTGGG‑CATATCACCCTCGACCCCCACCGGTGGGGTGATAAACCGCTCCCTCAGATAGCGTTCTGGCGGGCGCACAttgg                                          >  2:124863/1‑248 (MQ=255)                                                                                                                                                    ggCGTCGCCATTCTCTGGCCCCAGCATGTAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTGCAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCGTTGTG‑CGTATCACCCTCGACCCACACCGGTGAGGTGATAAGCCGCTCCGTCAGATAGCGTTCTGGCGGGCGCACGTCGGCCt                                       <  2:106304/251‑1 (MQ=255)                                                                                                                                                                                        gcgcTGCCCGCAGCCGTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTG‑CATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCttt    <  1:208998/250‑1 (MQ=255)                                                                                                                                                                                         cgcTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTG‑CATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTc  <  1:53938/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                               GGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATA‑GTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTG‑CATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTC  >  NC_000913/65070‑65504

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