Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	65,135	C→T	51.2%	D216N (GAT→AAT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	65,135	0	C	T	51.2%	‑5.6 / 63.1	43	D216N (GAT→AAT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (5/17);  total (13/30)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACG  >  W3110S.gb/65069‑65190                                                                   |                                                        cGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGt                                                       <  1:3564281/69‑1 (MQ=255) cGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGt                                                       <  1:813693/69‑1 (MQ=255) cGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGt                                                       <  1:3452102/69‑1 (MQ=255) cGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGt                                                       <  1:1288693/69‑1 (MQ=255) cGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGt                                                       <  1:1932691/69‑1 (MQ=255)     aTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   <  1:523810/69‑1 (MQ=255)     aTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   <  1:2776799/69‑1 (MQ=255)     aTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   <  1:3511729/69‑1 (MQ=255)     aTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGtt                                                   <  1:1222008/69‑1 (MQ=255)       gTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTgg                                                 <  1:907328/69‑1 (MQ=255)             gcagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAAc                                            <  1:922083/68‑1 (MQ=255)             gcagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAAc                                            <  1:297619/68‑1 (MQ=255)             gcagcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAAc                                            <  1:957605/68‑1 (MQ=255)                        aGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGt                                <  1:1612197/69‑1 (MQ=255)                        aGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGt                                <  1:1705822/69‑1 (MQ=255)                        aGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGt                                <  1:3507028/69‑1 (MQ=255)                         gATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGt                                <  1:143717/68‑1 (MQ=255)                            cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCa                                               >  1:3527591/1‑50 (MQ=255)                            cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCa                                               >  1:2696904/1‑50 (MQ=255)                               aCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTcc                          <  1:1141820/68‑1 (MQ=255)                               aCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTcc                          <  1:3003028/68‑1 (MQ=255)                                        cGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcg                <  1:1185397/69‑1 (MQ=255)                                        cGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgc                 >  1:3293927/1‑68 (MQ=255)                                          ttCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcggg              >  1:2094659/1‑69 (MQ=255)                                          ttCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcggg              >  1:600928/1‑69 (MQ=255)                                            ccAACCGATGATCACATCAGGATTGTAGTTGGCAAAcc                                           <  1:899624/38‑1 (MQ=255)                                            ccAACCGATGATCACATCAGGATTGTAGTTGGCAAAcc                                           <  1:651297/38‑1 (MQ=255)                                            ccAACCGATGATCACATCAGGATTGTAGTTGGCAAAcc                                           <  1:1918581/38‑1 (MQ=255)                                                  gatgatCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       >  1:1064614/1‑52 (MQ=255)                                                  gatgatCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       >  1:3311344/1‑52 (MQ=255)                                                  gatgatCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       >  1:3330385/1‑52 (MQ=255)                                                  gatgatCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       <  1:618761/52‑1 (MQ=255)                                                     gatCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGAc   <  1:2437698/69‑1 (MQ=255)                                                      atCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg  <  1:2085501/69‑1 (MQ=255)                                                              aGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcgggcg            >  1:3582300/1‑51 (MQ=255)                                                              aGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTgcgggcg            >  1:2089701/1‑51 (MQ=255)                                                              aGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGa    <  1:487444/59‑1 (MQ=255)                                                              aGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGa    <  1:2839233/59‑1 (MQ=255)                                                              aGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGa    <  1:2427701/59‑1 (MQ=255)                                                               ggATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTcc                          <  1:1357625/36‑1 (MQ=255)                                                                 aTCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       >  1:2021077/1‑37 (MQ=255)                                                                 aTCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       >  1:1396195/1‑37 (MQ=255)                                                                 aTCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                       >  1:2672518/1‑37 (MQ=255)                                                                   |                                                        CGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACG  >  W3110S.gb/65069‑65190

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