Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,366,654 T→C 20.0% Y237H (TAT→CAT)  puuB → gamma‑Glu‑putrescine oxidase, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,366,6540TC20.0% 49.1 / 8.6 25Y237H (TAT→CAT) puuBgamma‑Glu‑putrescine oxidase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (14/6);  new base C (3/2);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCAC  >  W3110S.gb/1366602‑1366721
                                                    |                                                                   
aaCTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTag                                                   <  1:2499575/71‑1 (MQ=255)
       aaaGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGa                                            <  1:2898022/71‑1 (MQ=255)
          ggTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACt                                         >  1:2655050/1‑71 (MQ=255)
          ggTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACt                                         >  1:187770/1‑71 (MQ=255)
                gTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAg                                                      <  1:1852574/52‑1 (MQ=255)
                gTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGa                                    >  1:921177/1‑70 (MQ=255)
                gTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGa                                    >  1:2507828/1‑70 (MQ=255)
                      gCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAATTAGAGCCGGAACTGGCGAAACgcag                             <  1:4016362/71‑1 (MQ=255)
                      gCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcag                             <  1:478095/71‑1 (MQ=255)
                        gAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTa                                                    >  1:3783332/1‑46 (MQ=255)
                             aTGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTa                                                    >  1:1822206/1‑41 (MQ=255)
                             aTGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTa                                                    >  1:1334450/1‑41 (MQ=255)
                             aTGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTa                                                    >  1:3000424/1‑41 (MQ=255)
                             aTGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTgg                                       >  1:839926/1‑54 (MQ=255)
                             aTGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTgg                                       >  1:2682977/1‑54 (MQ=255)
                              tgtgATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGaaa                                  >  1:1763186/1‑58 (MQ=255)
                               gtgATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTg                    >  1:3823670/1‑71 (MQ=255)
                                 gATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGcc                       >  1:1923815/1‑66 (MQ=255)
                                 gATTGTCGCCGGCAATGCGTATCTGTGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGcc                       >  1:1408380/1‑66 (MQ=255)
                                        gCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc            <  1:2827982/70‑1 (MQ=255)
                                        gCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc            <  1:1149443/70‑1 (MQ=255)
                                            ggAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTg       <  1:156057/71‑1 (MQ=255)
                                                  cGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcac  >  1:3153934/1‑70 (MQ=255)
                                                  cGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcac  >  1:679181/1‑70 (MQ=255)
                                                  cGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcac  >  1:2011485/1‑70 (MQ=255)
                                                    |                                                                   
AACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCAC  >  W3110S.gb/1366602‑1366721

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: