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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb65,1350CT89.3% 71.6 / 1.9 28D216N (GAT→AAT) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (11/14);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAA  >  W3110S.gb/65067‑65205
                                                                    |                                                                      
cTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGt                                                                      <  1:3561202/71‑1 (MQ=255)
cTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGt                                                                      <  1:1869335/71‑1 (MQ=255)
cTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGt                                                                      <  1:3305414/71‑1 (MQ=255)
cTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGt                                                                      <  1:187736/71‑1 (MQ=255)
           tttGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAAc                                                           >  1:971116/1‑71 (MQ=255)
                agcaTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAgg                                                       <  1:1078216/70‑1 (MQ=255)
                      ggCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGt                                                                   >  1:867354/2‑52 (MQ=255)
                      cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAg                                                                    >  1:2248312/1‑51 (MQ=255)
                      cGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAg                                                                    >  1:3239380/1‑51 (MQ=255)
                             cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                                     >  1:2716719/1‑43 (MQ=255)
                             cGAACTGCACCACGTTCCAACCGATGATCACATCAGGATTGTa                                                                     >  1:1191607/1‑43 (MQ=255)
                               aaCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGtt                                                                  <  1:162612/44‑1 (MQ=255)
                               aaCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGtt                                                                  <  1:1789447/44‑1 (MQ=255)
                               aaCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGtt                                                                  <  1:1080067/44‑1 (MQ=255)
                               aaCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGtt                                                                  <  1:3228882/44‑1 (MQ=255)
                               aaCTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGtt                                                                  <  1:3495826/44‑1 (MQ=255)
                                aCTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                                      <  1:797395/71‑1 (MQ=255)
                                 cTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                                      >  1:1664850/1‑70 (MQ=255)
                                 cTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                                      >  1:1455672/1‑70 (MQ=255)
                                 cTGCACCACGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGc                                      >  1:2343542/1‑70 (MQ=255)
                                    caccacGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACt                                  <  1:2939665/71‑1 (MQ=255)
                                     accacGTTCCAACCGATGATCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACt                                  <  1:1069823/70‑1 (MQ=255)
                                                     tgatCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACg                 <  1:1641291/71‑1 (MQ=255)
                                                        tCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTAt              <  1:2785541/71‑1 (MQ=255)
                                                        tCACATCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTAt              <  1:956842/71‑1 (MQ=255)
                                                           caTCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACt                                  <  1:942677/48‑1 (MQ=255)
                                                           caTCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGggcggc                         >  1:2196047/1‑57 (MQ=255)
                                                           caTCAGGATTGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGggcggc                         >  1:162058/1‑57 (MQ=255)
                                                                    cGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGAGGGCGGCTGGCGACGTATTCCAGTTCGaaa  >  1:360707/1‑71 (MQ=255)
                                                                    |                                                                      
CTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAA  >  W3110S.gb/65067‑65205

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