Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	61,962	A→G	12.3%	F499F (TTT→TTC)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	61,962	0	A	G	12.3%	149.2 / 11.5	57	F499F (TTT→TTC)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (31/19);  new base G (2/5);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

TGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCA  >  Exported/61898‑62025                                                                 |                                                                tGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCa                                                           <  1:4895270/71‑1 (MQ=255) tGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCa                                                           <  1:4050852/71‑1 (MQ=255)      aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:529948/71‑1 (MQ=255)      aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:4960578/71‑1 (MQ=255)      aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:1848242/71‑1 (MQ=255)      aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:254312/71‑1 (MQ=255)        cgacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga                                                    >  1:848268/1‑71 (MQ=255)        cgacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga                                                    >  1:2847933/1‑71 (MQ=255)         gacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga                                                    <  1:2108599/70‑1 (MQ=255)               gCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt                                             >  1:358330/1‑71 (MQ=255)               gCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt                                             >  1:4742047/1‑71 (MQ=255)               gCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt                                             >  1:2344276/1‑71 (MQ=255)                ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:3155193/71‑1 (MQ=255)                ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:3149026/71‑1 (MQ=255)                ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:4210221/71‑1 (MQ=255)                ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:4216314/71‑1 (MQ=255)                ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:2149633/71‑1 (MQ=255)                   gCCGCGGATCGAAGAAGCGTCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc                                         >  1:2545034/1‑71 (MQ=255)                   gCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc                                         >  1:3147449/1‑71 (MQ=255)                   gCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc                                         >  1:2699352/1‑71 (MQ=255)                       cgGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCt                                     <  1:4821529/71‑1 (MQ=255)                       cgGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCt                                     <  1:2774116/71‑1 (MQ=255)                       cgGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCt                                     <  1:2998960/71‑1 (MQ=255)                         gATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGc                                   <  1:2332710/71‑1 (MQ=255)                          aTCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:2889953/51‑1 (MQ=255)                            cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:1743065/1‑70 (MQ=255)                            cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:1225392/1‑70 (MQ=255)                            cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:4422257/1‑70 (MQ=255)                            cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:3991775/1‑70 (MQ=255)                            cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTAATGATGATTTTCAGCGCCTGCGa                                 >  1:3824803/1‑70 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:3820330/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:4069394/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:4431695/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:1185120/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:2297980/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:2522947/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:1553632/1‑52 (MQ=255)                              aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:1290761/1‑52 (MQ=255)                                gaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatgat                                                   <  1:252617/48‑1 (MQ=255)                                  agCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGt                          >  1:381994/1‑71 (MQ=255)                                   gCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATTATTTTCATCGCCt                                     <  1:4741629/59‑1 (MQ=255)                                      gCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGt                      <  1:2124368/71‑1 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:101698/1‑46 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:2733028/1‑46 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3792336/1‑46 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3657449/1‑46 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3137110/1‑46 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3415326/1‑46 (MQ=255)                                           cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:814777/1‑46 (MQ=255)                                                 tgCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTTCGACAGCGGTTTGTTACCCTGGCGt           <  1:615969/71‑1 (MQ=255)                                                  gCCGAGCACGCCGTAGAAGGCGTTCATGATTATCTTCAGCGCCTGCGa                                 >  1:531144/1‑48 (MQ=255)                                                  gCCGAGCACGCCGTAGAAGGCGTTCATGATTATCTTCAGCGCCTGCGa                                 >  1:922626/1‑48 (MQ=255)                                                  gCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCa                                           >  1:2499409/1‑38 (MQ=255)                                                    cGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTg        >  1:1223149/1‑71 (MQ=255)                                                      aGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc      <  1:4997668/71‑1 (MQ=255)                                                          cGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCa  <  1:3058615/71‑1 (MQ=255)                                                               tAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGTCGTTTg        <  1:670793/60‑1 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTa                    >  1:3349075/1‑46 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTAc                   >  1:3105215/1‑47 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGtt          >  1:4124873/1‑56 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGtt          >  1:319388/1‑56 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGtt          >  1:3171794/1‑56 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc      >  1:2624123/1‑60 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc      >  1:1547859/1‑60 (MQ=255)                                                                 aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc      >  1:2630115/1‑60 (MQ=255)                                                                 |                                                                TGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCA  >  Exported/61898‑62025

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