Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,066,786 C→T 66.7% intergenic (+144/‑129) asd → / → glgB aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,066,7860CT66.7% 56.0 / 45.1 51intergenic (+144/‑129)asd/glgBaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (11/6);  new base T (21/13);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCAAAAA  >  W3110S.gb/4066726‑4066849
                                                            |                                                               
gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTg                                                         <  1:1197606/69‑1 (MQ=255)
gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTg                                                         <  1:815955/69‑1 (MQ=255)
gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTg                                                         <  1:1827686/69‑1 (MQ=255)
gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTg                                                         <  1:164239/69‑1 (MQ=255)
    gTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGt                                                     <  1:2624546/69‑1 (MQ=255)
       aaaTAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAg                                                  <  1:2938417/69‑1 (MQ=255)
          tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGgac                                               <  1:739917/69‑1 (MQ=255)
           agagCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                                              >  1:3688771/1‑53 (MQ=255)
           agagCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                                              >  1:95481/1‑53 (MQ=255)
           agagCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                                               >  1:3253520/1‑52 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAgg                                                 >  1:2119584/1‑58 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCagag                                      >  1:161936/1‑69 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCagag                                      >  1:1485086/1‑69 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCagag                                      >  1:2083543/1‑69 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTAACTGAAGTGAgg                                                 >  1:3548340/1‑58 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCAAAGTCTTACCTg                                                         >  1:892856/1‑50 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGATTCACAGCCTTACCTGAAGTGAGGACGACGCagag                                      >  1:1545108/1‑68 (MQ=255)
                   ggCTATTCTTTAAGGGTGGCTGAATACATGAGATTCACAGCCTTACCTGAAGTGAGGACGACGCagag                                      >  1:2744748/1‑68 (MQ=255)
                      tATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                    <  1:474670/68‑1 (MQ=255)
                      tATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                    <  1:1150779/68‑1 (MQ=255)
                      tATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                    <  1:1929793/68‑1 (MQ=255)
                       aTTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGa                                  >  1:2336785/1‑69 (MQ=255)
                             tAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                              >  1:2619873/1‑67 (MQ=255)
                                  gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtg                       <  1:2839147/69‑1 (MQ=255)
                                  gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtg                       <  1:3317941/69‑1 (MQ=255)
                                     gCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgctg                    >  1:2504813/1‑69 (MQ=255)
                                     gCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgctg                    >  1:2762502/1‑69 (MQ=255)
                                          aTACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCg               >  1:967182/1‑69 (MQ=255)
                                          aTACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCg               >  1:1439812/1‑69 (MQ=255)
                                          aTACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCg               >  1:1585697/1‑69 (MQ=255)
                                           tACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCg               >  1:115318/1‑68 (MQ=255)
                                           tACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCg               >  1:844601/1‑68 (MQ=255)
                                             cATGAGTATTCACAGTCTTACCTGAAGTGAGGacgacg                                           >  1:3037346/1‑38 (MQ=255)
                                             cATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTc            >  1:1845399/1‑69 (MQ=255)
                                             cATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTc            >  1:325959/1‑69 (MQ=255)
                                             cATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTc            >  1:3373570/1‑69 (MQ=255)
                                             cATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTc            >  1:71085/1‑69 (MQ=255)
                                                 aGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAgg         >  1:1730094/1‑68 (MQ=255)
                                                 aGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAgg         >  1:948786/1‑68 (MQ=255)
                                                   tATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCa      <  1:3042673/69‑1 (MQ=255)
                                                     ttCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaa     >  1:1639430/1‑68 (MQ=255)
                                                       cacaGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  <  1:2288347/69‑1 (MQ=255)
                                                       cacaGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  >  1:1427551/1‑69 (MQ=255)
                                                       cacaGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  >  1:478300/1‑69 (MQ=255)
                                                       cacaGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  <  1:1923747/69‑1 (MQ=255)
                                                       cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATg                                <  1:568502/39‑1 (MQ=255)
                                                       cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATg                                <  1:3672454/39‑1 (MQ=255)
                                                       cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCTTTCAGGTCaaaaa  <  1:417861/69‑1 (MQ=255)
                                                       cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  >  1:2226293/1‑69 (MQ=255)
                                                       cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  >  1:1012271/1‑69 (MQ=255)
                                                       cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaa  <  1:349472/69‑1 (MQ=255)
                                                            |                                                               
GGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCAAAAA  >  W3110S.gb/4066726‑4066849

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: