Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,840,761 A→G 74.1% 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,7610AG74.1% 98.3 / 32.6 54E372E (GAA→GAGascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (9/5);  new base G (18/22);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCA  >  W3110S.gb/2840698‑2840825
                                                               |                                                                
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACgg                                                             >  1:1084067/1‑69 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACgg                                                             >  1:701528/1‑69 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACgg                                                             >  1:267648/1‑69 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACgg                                                             >  1:2304471/1‑69 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACg                                                              >  1:2498185/1‑68 (MQ=255)
   tCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGcc                                                           <  1:2576093/68‑1 (MQ=255)
   tCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGc                                                            <  1:970255/67‑1 (MQ=255)
         tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTggg                                                       <  1:3142693/66‑1 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGgcgc                                                    >  1:643306/1‑65 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGgcgc                                                    >  1:2939102/1‑65 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGgcgc                                                    >  1:315845/1‑65 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGgcgc                                                    >  1:1344186/1‑65 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGgcg                                                     >  1:3698601/1‑64 (MQ=255)
                   aaTATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGa                                           >  1:2428738/1‑68 (MQ=255)
                      atgatgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAttt                                       <  1:1971213/69‑1 (MQ=255)
                      atgatgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAttt                                       <  1:2769745/69‑1 (MQ=255)
                      atgatgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAttt                                       <  1:3153077/69‑1 (MQ=255)
                      atgatgTACGAGCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAttt                                       <  1:1987002/69‑1 (MQ=38)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 <  1:3426267/69‑1 (MQ=255)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 >  1:3097541/1‑69 (MQ=255)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 >  1:2857910/1‑69 (MQ=255)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 <  1:2801216/69‑1 (MQ=255)
                                  cGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGc                           <  1:80668/69‑1 (MQ=255)
                                     tATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                         <  1:3322219/68‑1 (MQ=255)
                                     tATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                         <  1:3624331/68‑1 (MQ=255)
                                     tATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                         <  1:1057611/68‑1 (MQ=255)
                                      aTCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 >  1:3384601/1‑59 (MQ=255)
                                      aTCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCgag                        >  1:2467669/1‑68 (MQ=255)
                                      aTCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCgag                        >  1:527713/1‑68 (MQ=255)
                                         aGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAtt                     <  1:2634018/68‑1 (MQ=255)
                                         aGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAtt                     <  1:2073516/68‑1 (MQ=255)
                                         aGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 <  1:2139206/56‑1 (MQ=255)
                                         aGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 <  1:352871/56‑1 (MQ=255)
                                             gCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacg                 >  1:3449232/1‑68 (MQ=255)
                                              ccGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgac               <  1:906316/69‑1 (MQ=255)
                                              ccGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgac               <  1:2130830/69‑1 (MQ=255)
                                              ccGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgac               <  1:2353751/69‑1 (MQ=255)
                                                  tGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt           <  1:492809/69‑1 (MQ=255)
                                                  tGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt           <  1:3015397/69‑1 (MQ=255)
                                                  tGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt           <  1:1407015/69‑1 (MQ=255)
                                                  tGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt           <  1:3459354/69‑1 (MQ=255)
                                                   gTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATg                                            >  1:1761517/1‑35 (MQ=255)
                                                   gTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATg                                            >  1:1477748/1‑35 (MQ=255)
                                                   gTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATg                                            >  1:3134361/1‑35 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                  >  1:1098086/1‑58 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt           >  1:2933501/1‑65 (MQ=255)
                                                      tCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg       >  1:3649382/1‑69 (MQ=255)
                                                      tCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg       >  1:2444769/1‑69 (MQ=255)
                                                        tggtggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGc      <  1:1150748/68‑1 (MQ=255)
                                                        tggtggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTACCAATGGCGAGATTAACGACTACTATCGc      <  1:3059795/68‑1 (MQ=255)
                                                         ggtggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa     <  1:1177663/68‑1 (MQ=255)
                                                            ggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCa  >  1:3141268/1‑68 (MQ=255)
                                                             gAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 >  1:603077/1‑36 (MQ=255)
                                                             gAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                 >  1:930228/1‑36 (MQ=255)
                                                               |                                                                
CGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCA  >  W3110S.gb/2840698‑2840825

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: