Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,366,654 T→C 25.2% 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,6540TC25.2% 68.0 / 21.5 41Y237H (TAT→CAT) puuBgamma‑Glu‑putrescine oxidase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (12/18);  new base C (4/6);  total (17/24)
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

GTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCA  >  W3110S.gb/1366597‑1366723
                                                         |                                                                     
gtgCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATa                                                           >  1:4326460/1‑70 (MQ=255)
gtgCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATa                                                           >  1:2620182/1‑70 (MQ=255)
          ccAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAg                                                        <  1:2004704/63‑1 (MQ=255)
          ccAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCg                                                >  1:3928794/1‑71 (MQ=255)
           cAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCgg                                               <  1:4044004/71‑1 (MQ=255)
            aaaGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATaaa                                                         <  1:1055092/60‑1 (MQ=255)
            aaaGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTa                                                      >  1:214589/1‑63 (MQ=255)
            aaaGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGt                                                       >  1:1078920/1‑62 (MQ=255)
            aaaGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGa                                              <  1:2950326/71‑1 (MQ=255)
             aaGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGaa                                             <  1:3661253/71‑1 (MQ=255)
               ggTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAAc                                            >  1:4527151/1‑70 (MQ=255)
               ggTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACt                                           >  1:2150436/1‑71 (MQ=255)
                     gTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATaaa                                                         >  1:2686534/1‑51 (MQ=255)
                     gTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTg                                          <  1:3974823/66‑1 (MQ=255)
                     gTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGa                                                             <  1:1950540/47‑1 (MQ=255)
                      tGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGa                                                             >  1:1234433/1‑46 (MQ=255)
                      tGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGaaa                                    >  1:1331120/1‑71 (MQ=255)
                        aCGGCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTagag                                                   >  1:991145/1‑54 (MQ=255)
                          ggCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATaa                                                          >  1:258730/1‑45 (MQ=255)
                          ggCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACg                                  >  1:1571179/1‑69 (MQ=255)
                           gCGAAGTATGTGATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcag                               <  1:4049304/71‑1 (MQ=255)
                             gAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcagca                             <  1:2370562/71‑1 (MQ=255)
                                gTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCg                                                <  1:149884/49‑1 (MQ=255)
                                   tgtgATTGTCGCCGGGAATGCGTATCTGGGCGATAAAg                                                        <  1:139410/38‑1 (MQ=255)
                                     tgATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTgg                                         >  1:1540965/1‑51 (MQ=255)
                                     tgATTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGc                     >  1:2352839/1‑71 (MQ=255)
                                      gATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcag                               <  1:3333691/60‑1 (MQ=255)
                                      gATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcag                               <  1:3780097/60‑1 (MQ=255)
                                      gATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcag                               <  1:543344/60‑1 (MQ=255)
                                       aTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGt                                                       <  1:4150580/35‑1 (MQ=255)
                                       aTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCgg                                               <  1:4570969/43‑1 (MQ=255)
                                       aTTGTCGCCGGGAATGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCgg                   <  1:3556586/71‑1 (MQ=255)
                                       aTTGTCGCCGGGAATGCGAATCTGGGCGATAAAGTa                                                      >  1:4449933/1‑36 (MQ=255)
                                             gCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc              <  1:2385332/70‑1 (MQ=255)
                                             gCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCa             <  1:355924/71‑1 (MQ=255)
                                                  gAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGa        <  1:4101401/71‑1 (MQ=255)
                                                    aTGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCgg                   <  1:1911991/58‑1 (MQ=255)
                                                    aTGCGCATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGAt       <  1:1889229/70‑1 (MQ=255)
                                                      gcgcATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcag                               <  1:3408324/44‑1 (MQ=255)
                                                      gCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCActac             <  1:3821369/62‑4 (MQ=255)
                                                       cGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcac    >  1:3495054/1‑70 (MQ=255)
                                                        gTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcacca  <  1:1435921/71‑1 (MQ=255)
                                                         |                                                                     
GTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCA  >  W3110S.gb/1366597‑1366723

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: