Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,840,761 A→G 71.9% E372E (GAA→GAG ascB → cryptic 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,840,7610AG71.9% 50.2 / 20.2 32E372E (GAA→GAGascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (14/9);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCG  >  W3110S.gb/2840698‑2840820
                                                               |                                                           
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                      >  1:1213045/1‑71 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                      >  1:938000/1‑71 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                      >  1:648271/1‑71 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                      >  1:669999/1‑71 (MQ=255)
     tGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGgc                                                 >  1:1498345/1‑71 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                             >  1:427212/1‑69 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                             >  1:587331/1‑69 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                             >  1:1547515/1‑69 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                             >  1:1451773/1‑69 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                             >  1:1296628/1‑69 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCa                                              >  1:1341958/1‑68 (MQ=255)
           tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCATGGGCGCaa                                             >  1:1854468/1‑69 (MQ=255)
                aTGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATg                                       >  1:713523/1‑70 (MQ=255)
                aTGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATg                                       >  1:1822134/1‑70 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:1674967/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:922986/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:977780/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:498328/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:1907849/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGATGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:1275478/71‑1 (MQ=255)
                          tgTACGAGCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:1912188/71‑1 (MQ=255)
                                             gCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                            <  1:1553476/52‑1 (MQ=255)
                                                gCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt      <  1:369712/71‑1 (MQ=255)
                                                gCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt      <  1:271488/71‑1 (MQ=255)
                                                gCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt      <  1:507741/71‑1 (MQ=255)
                                                   gTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCaaat                        >  1:311496/1‑48 (MQ=39)
                                                   gTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATc   >  1:1526470/1‑71 (MQ=255)
                                                   gTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATc   >  1:1277087/1‑71 (MQ=255)
                                                   gTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATc   >  1:1264681/1‑71 (MQ=255)
                                                    tttGTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg  <  1:1409383/71‑1 (MQ=255)
                                                    tttCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg  <  1:33804/71‑1 (MQ=255)
                                                    tttCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg  <  1:1817503/71‑1 (MQ=255)
                                                               |                                                           
CGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCG  >  W3110S.gb/2840698‑2840820

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: