Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,066,786 C→T 100% 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,7860CT80.8% 52.7 / 9.0 26intergenic (+144/‑129)asd/glgBaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base T (8/13);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 3.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCAAAAAAA  >  W3110S.gb/4066731‑4066851
                                                       |                                                                 
ttAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAg                                                    <  1:334068/71‑1 (MQ=255)
ttAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCCTACCTGAAGTGAg                                                    <  1:197364/71‑1 (MQ=255)
                   ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATg                                  <  1:444359/70‑1 (MQ=255)
                   ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATg                                  <  1:1883584/70‑1 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAgg                                     >  1:124287/1‑65 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:1446890/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:995055/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:1260340/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:207940/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:1604313/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:1548559/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTCAGGACGACGCAGAGAGGATGca                                >  1:1625237/1‑70 (MQ=255)
                     cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGca                                >  1:880401/1‑70 (MQ=255)
                        tAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACag                             <  1:1589285/70‑1 (MQ=255)
                        tAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACag                             <  1:571123/70‑1 (MQ=255)
                        tAAGGGTGGCTGAATACATGAGTATTCACAGTCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACag                             <  1:960015/70‑1 (MQ=255)
                             gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgct                       >  1:228296/1‑71 (MQ=255)
                               ggCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgctgc                     >  1:1457975/1‑71 (MQ=255)
                               ggCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgctgc                     >  1:1072169/1‑71 (MQ=255)
                                              tATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                      <  1:1330519/39‑1 (MQ=255)
                                              tATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                      <  1:1273372/39‑1 (MQ=255)
                                              tATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                      <  1:709757/39‑1 (MQ=255)
                                              tATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                      <  1:733206/39‑1 (MQ=255)
                                              tATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                      <  1:823814/39‑1 (MQ=255)
                                              tATTCACAGTCTTACCTGAAGTGAGGACGACGCagagag                                      <  1:140575/39‑1 (MQ=255)
                                                  cacaGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCaaaaaaa  >  1:782860/1‑71 (MQ=255)
                                                       |                                                                 
TTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGCTGCGCCGTTCAGGTCAAAAAAA  >  W3110S.gb/4066731‑4066851

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: