Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,840,761 A→G 89.8% 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,7610AG89.8% 176.0 / 8.8 59E372E (GAA→GAGascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (26/27);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCA  >  W3110S.gb/2840698‑2840822
                                                               |                                                             
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGcc                                                        >  1:3175167/1‑71 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGcc                                                        >  1:1306599/1‑71 (MQ=255)
cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGc                                                         >  1:3085762/1‑70 (MQ=255)
         tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGTCCTGGGCGCaa                                               <  1:989601/71‑1 (MQ=255)
         tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                               <  1:3058619/71‑1 (MQ=255)
         tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                               <  1:664937/71‑1 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGagaa                                                             >  1:905480/1‑53 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCaa                                               >  1:476304/1‑67 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCaa                                               >  1:1199818/1‑67 (MQ=255)
             aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCaa                                               >  1:551438/1‑67 (MQ=255)
                aTGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGACGCAAAAGATg                                         >  1:3015769/1‑70 (MQ=255)
                 tGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGaa                                       <  1:250722/71‑1 (MQ=255)
                 tGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGTGCGCAAAAGATGaa                                       <  1:647456/71‑1 (MQ=255)
                 tGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                       <  1:3694141/71‑1 (MQ=255)
                 tGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                       <  1:3537263/71‑1 (MQ=255)
                 tGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                       <  1:925862/71‑1 (MQ=255)
                   aaTATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAtt                                     >  1:2412432/1‑71 (MQ=255)
                   aaTATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAtt                                     >  1:330174/1‑71 (MQ=255)
                       tgatgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAt                                      >  1:263449/1‑66 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:2229636/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:3113706/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:3312889/71‑1 (MQ=255)
                          tgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:143052/71‑1 (MQ=255)
                           gTACGTTCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:1424907/70‑1 (MQ=255)
                           gTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              >  1:2748207/1‑70 (MQ=255)
                           gTACGATAGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              >  1:1241810/1‑70 (MQ=255)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCaa                            <  1:3060360/71‑1 (MQ=255)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCaa                            <  1:2890892/71‑1 (MQ=255)
                            tACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCaa                            <  1:2081117/71‑1 (MQ=255)
                                  cGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGc                        >  1:655412/1‑69 (MQ=255)
                                  cGTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGc                        >  1:2568758/1‑69 (MQ=255)
                                   gTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCaa                            >  1:2875847/1‑64 (MQ=255)
                                   gTTATCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                      <  1:2292421/70‑1 (MQ=255)
                                   gTTATAAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCaa                            >  1:2715679/1‑64 (MQ=255)
                                      aTCAGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAt                   >  1:3046835/1‑70 (MQ=255)
                                         aGAAGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac               >  1:170270/1‑71 (MQ=255)
                                           aaGCCGCTGTTTCTGGTGGAGAACGGCCTTGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacga             >  1:3588274/1‑71 (MQ=255)
                                           aaGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                      <  1:2071026/62‑1 (MQ=255)
                                           aaGCCGCTGTTTCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacga             >  1:1100832/1‑71 (MQ=255)
                                                    tttCTGGTGGAGAACGGCCTGGGCGCAAAAGATGa                                        <  1:941928/35‑1 (MQ=255)
                                                    tttCTGGTGGAGAACGGCCTGGGCGCAAAAGATGa                                        <  1:2441073/35‑1 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAAt                                      >  1:21929/1‑35 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:2801436/1‑71 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:1098744/1‑71 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:2938143/1‑71 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:1255124/1‑71 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:1268576/1‑71 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:2785949/1‑71 (MQ=255)
                                                      tCTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  >  1:1706385/1‑71 (MQ=255)
                                                       cTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  <  1:2603691/70‑1 (MQ=255)
                                                       cTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  <  1:647429/70‑1 (MQ=255)
                                                       cTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  <  1:3200956/70‑1 (MQ=255)
                                                       cTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  <  1:2779614/70‑1 (MQ=255)
                                                       cTGGTGGAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCa  <  1:2903610/70‑1 (MQ=255)
                                                           tggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:509851/38‑1 (MQ=255)
                                                           tggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:2394010/38‑1 (MQ=255)
                                                           tggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:1777060/38‑1 (MQ=255)
                                                           tggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:1439847/38‑1 (MQ=255)
                                                           tggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:1002378/38‑1 (MQ=255)
                                                           tggAGAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                              <  1:3198972/38‑1 (MQ=255)
                                                               |                                                             
CGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCA  >  W3110S.gb/2840698‑2840822

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: