Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,347,163 A→G 10.0% P191P (CCA→CCG melA → alpha‑galactosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,347,1630AG10.0% 184.2 / 10.1 70P191P (CCA→CCGmelAalpha‑galactosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (12/51);  new base G (2/5);  total (14/56)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCAAAACC  >  W3110S.gb/4347095‑4347232
                                                                    |                                                                     
cAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGc                                                                     <  1:1105250/71‑1 (MQ=255)
cAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGc                                                                     <  1:2789693/71‑1 (MQ=255)
    tCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACg                                                                 >  1:3857099/1‑71 (MQ=255)
    tCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACg                                                                 >  1:1201008/1‑71 (MQ=255)
          tgtgCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGt                                                           >  1:1475682/1‑71 (MQ=255)
              ccATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGc                                                                     <  1:464728/57‑1 (MQ=255)
                aTTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGt                                                     >  1:1845971/1‑71 (MQ=255)
                      tGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCa                                               <  1:1104536/71‑1 (MQ=255)
                          gggAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTa                                           <  1:2012970/71‑1 (MQ=255)
                           ggAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTgcg                                                 >  1:2598449/1‑64 (MQ=255)
                            gAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATc                                         <  1:2664296/71‑1 (MQ=255)
                            gAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATc                                         <  1:1031829/71‑1 (MQ=255)
                             aaCGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTa                                                                   >  1:868738/1‑44 (MQ=255)
                             aaCGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTAcc                                                       <  1:1918294/56‑1 (MQ=255)
                             aaCGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTAcc                                                       <  1:3857608/56‑1 (MQ=255)
                             aaCGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTAcc                                                       <  1:3853438/56‑1 (MQ=255)
                                ggcggAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTAcc                                                       <  1:3726180/53‑1 (MQ=255)
                                ggcggAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAAcc                                     <  1:589114/71‑1 (MQ=255)
                                ggcggAAGAGCTGGCGCGTGACCTCAATATCGACCCGGCTACGCTGc                                                             <  1:631911/47‑1 (MQ=38)
                                ggcggAAGAGCTGGCGCGTGACCTCAATATCGACCCGGCTACGCTGc                                                             <  1:692979/47‑1 (MQ=38)
                                ggcggAAGAGCTGGCGCGTGACCTCAATATCGACCCGGCTACGCTGc                                                             <  1:2016712/47‑1 (MQ=38)
                                ggcggAAGAGCTGGCGCGTGACCTCAATATCGACCCGGCTACGCTGc                                                             <  1:2435843/47‑1 (MQ=38)
                                ggcggAAGAGCTGGCGCCTGACCTCAATATCGACCCGGCTACGCTGc                                                             <  1:1171026/47‑1 (MQ=38)
                                 gcggAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:574258/71‑1 (MQ=255)
                                 gcggAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:934856/71‑1 (MQ=255)
                                  cggAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCat                                   <  1:585930/71‑1 (MQ=255)
                                    gAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCg                                                      <  1:1571980/50‑1 (MQ=255)
                                     aaGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATg                                <  1:2438611/71‑1 (MQ=255)
                                     aaGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATg                                <  1:1699292/71‑1 (MQ=255)
                                        agTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCg                             <  1:1942630/71‑1 (MQ=255)
                                        agTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCg                             <  1:1128791/71‑1 (MQ=255)
                                        agTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCg                             <  1:3339676/71‑1 (MQ=255)
                                         gTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCg                             <  1:699610/70‑1 (MQ=255)
                                          ttGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAAcaa                                    <  1:199360/62‑3 (MQ=255)
                                          ttGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:3637504/62‑1 (MQ=255)
                                          ttGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:3633071/62‑1 (MQ=255)
                                          ttGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTAACGTTGCGCAGGTATCAACCa                                    <  1:3852198/62‑1 (MQ=255)
                                           tGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCat                                   <  1:2648517/62‑1 (MQ=255)
                                            ggCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGttt                          <  1:937651/70‑1 (MQ=255)
                                            ggCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGttt                          <  1:870079/70‑1 (MQ=255)
                                            ggCGCGTGACATCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGttt                          <  1:3061859/70‑1 (MQ=255)
                                             gcgcGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGtttt                         <  1:56220/70‑1 (MQ=255)
                                              cgcgTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTAc                       <  1:2050226/71‑1 (MQ=255)
                                              cgcgTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTAc                       <  1:580934/71‑1 (MQ=255)
                                               gcgTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTAcc                      <  1:3503731/71‑1 (MQ=255)
                                               gcgTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTAc                       <  1:3296822/70‑1 (MQ=255)
                                                 gTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCata                                  >  1:2743535/1‑57 (MQ=255)
                                                 gTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCata                                  >  1:1296350/1‑57 (MQ=255)
                                                 gTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGAAGGTATCAACCata                                  >  1:1971731/1‑57 (MQ=255)
                                                 gTGACCTCAATATCGACACAGCTACGCTGCGTTACCGTTGCGCAg                                              >  1:181450/1‑45 (MQ=255)
                                                      cTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:316205/50‑1 (MQ=255)
                                                      cTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGtttt                         <  1:232645/61‑1 (MQ=255)
                                                      cTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACctgg                   <  1:327283/67‑1 (MQ=255)
                                                      cTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACctgg                   <  1:1750778/67‑1 (MQ=255)
                                                      cTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACctgg                   <  1:1533806/67‑1 (MQ=255)
                                                        cAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGtttt                         <  1:1134357/59‑1 (MQ=255)
                                                          atatCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACctgg                   <  1:651952/63‑1 (MQ=255)
                                                            atCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAgcg         >  1:2439156/1‑71 (MQ=255)
                                                             tCGACCCGGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAACTGGAgcg         >  1:3099301/1‑70 (MQ=255)
                                                             tCGACCCGGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAACTGGAgcg         >  1:455572/1‑70 (MQ=255)
                                                               gACCCAGCTACGCTGCGTTACCGTTGCTCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCa       >  1:3770825/1‑70 (MQ=255)
                                                               gACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCaa      <  1:583122/71‑1 (MQ=255)
                                                               gACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCaa      <  1:2851244/71‑1 (MQ=255)
                                                               gACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCaa      <  1:2611066/71‑1 (MQ=255)
                                                               gACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCaa      <  1:1097378/71‑1 (MQ=255)
                                                                   cAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:540272/37‑1 (MQ=255)
                                                                   cAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:516196/37‑1 (MQ=255)
                                                                   cAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCa                                    <  1:3858588/37‑1 (MQ=255)
                                                                   cAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCAAAAcc  <  1:815255/71‑1 (MQ=255)
                                                                   cAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCAAAAcc  <  1:3232344/71‑1 (MQ=255)
                                                                    |                                                                     
CAGGTCGGGCTGTGCCATTCGGTGCAGGGAACGGCGGAAGAGTTGGCGCGTGACCTCAATATCGACCCAGCTACGCTGCGTTACCGTTGCGCAGGTATCAACCATATGGCGTTTTACCTGGAGCTGGAGCGCAAAACC  >  W3110S.gb/4347095‑4347232

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: