Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,367,159 C→T 31.4% T405M (ACG→ATG)  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,367,1590CT31.4% 44.4 / 26.8 35T405M (ACG→ATG) puuBgamma‑Glu‑putrescine oxidase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base T (4/7);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATT  >  W3110S.gb/1367095‑1367228
                                                                |                                                                     
ggcgACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCg                                                                 <  1:3470456/71‑1 (MQ=255)
ggcgACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCg                                                                 <  1:1923658/71‑1 (MQ=255)
ggcgACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCg                                                                 <  1:2155683/71‑1 (MQ=255)
  cgACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCg                                                                 <  1:3183193/69‑1 (MQ=255)
         gAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCGTGTGCCGtt                                                        <  1:1629376/71‑1 (MQ=255)
         gAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGtt                                                        >  1:2895448/1‑71 (MQ=255)
                  gATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGAt                                               <  1:1805352/71‑1 (MQ=255)
                      ccTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATggg                                            <  1:312296/70‑1 (MQ=255)
                      ccTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACCCTGCGTGTGCCGTTTACCGCGATggg                                            <  1:3476612/70‑1 (MQ=255)
                          cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGtgtg                                                             >  1:1058980/1‑49 (MQ=255)
                          cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGtgtg                                                             >  1:1342549/1‑49 (MQ=255)
                          cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGtgtg                                                             >  1:328116/1‑49 (MQ=255)
                          cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGtgtg                                                             >  1:2823294/1‑49 (MQ=255)
                              aaTCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTa                                   >  1:504407/1‑71 (MQ=255)
                              aaTCTGCCGCATTACCCGTTCCCCGGCGGGCGCATGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTa                                   >  1:1309489/1‑71 (MQ=255)
                                tCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACctcg                                                 <  1:1109391/55‑1 (MQ=255)
                                  tGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTa                                >  1:492985/1‑70 (MQ=255)
                                          tACCCGTTCCCCGGCGGGCGCATGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTgcg                       >  1:2461269/1‑71 (MQ=255)
                                          tACCCGTTCCCCGGCGGGCGCATGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTgcg                       >  1:947549/1‑71 (MQ=255)
                                           aCCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTgc                        >  1:3073206/1‑69 (MQ=255)
                                            cccGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGcgcg                     <  1:615659/71‑1 (MQ=255)
                                            cccGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGAGGCGTATTACAGCCTgcgc                      >  1:909835/1‑70 (MQ=255)
                                               gTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATc                  >  1:924460/1‑71 (MQ=255)
                                               gTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATc                  >  1:861451/1‑71 (MQ=255)
                                                ttCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCg                 >  1:3226846/1‑71 (MQ=255)
                                                ttCCCCGGCGGGCGAATGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCg                 <  1:4083251/71‑1 (MQ=255)
                                                         gggcgCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGt        >  1:2638231/1‑71 (MQ=255)
                                                         gggcgCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGt        >  1:948925/1‑71 (MQ=255)
                                                           gcgcATGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTAtt                                 <  1:1612207/44‑1 (MQ=255)
                                                           gcgcACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGttt      <  1:1216162/71‑1 (MQ=255)
                                                           gcgcACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGttt      <  1:1843018/71‑1 (MQ=255)
                                                           gcgcACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGttt      <  1:1973478/71‑1 (MQ=255)
                                                           gcgcACGCCGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGttt      <  1:3327497/71‑1 (MQ=255)
                                                            cgcACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGttt      <  1:3818403/70‑1 (MQ=255)
                                                               aCGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAAtt  <  1:382679/71‑1 (MQ=255)
                                                                |                                                                     
GGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATT  >  W3110S.gb/1367095‑1367228

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: