Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 539,323 C→T 18.7% P318L (CCA→CTA)  allB → allantoinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb539,3230CT18.7% 66.8 / 10.1 32P318L (CCA→CTA) allBallantoinase
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base T (4/2);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

AACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAG  >  W3110S.gb/539278‑539389
                                             |                                                                  
aaCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGt                                           <  1:2972473/71‑1 (MQ=255)
aaCTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGt                                           <  1:2474722/71‑1 (MQ=255)
             tGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAg                              <  1:3173446/71‑1 (MQ=255)
                       tGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAAt                                                     <  1:2662568/38‑1 (MQ=255)
                          cTGGTTTCCGACCACTCTCTATGCCCGCCGGAAATGAAAGCCGGt                                           <  1:1247815/45‑1 (MQ=255)
                          cTGGTTTCCGACCACTCTCTATGCCCGCCGGAAATGAAAGCCGGt                                           <  1:902767/45‑1 (MQ=255)
                            ggTTTCCGACCACTCTCTATGCCCGCCGGAAATGAAAGCCGGTaa                                         >  1:280361/1‑45 (MQ=255)
                            ggTTTCCGACCACTCTCTATGCCCGCCGGAAATGAAAGCCGGTaa                                         >  1:2731677/1‑45 (MQ=255)
                            ggTTTCCGACCACTCTCTATGCCCGCCGGAAATGAAAGCCGGTaa                                         >  1:2481204/1‑45 (MQ=255)
                            ggTTTCCGACCACTCTCTATGCCCGCCGGAAATGAAAGCCGGTaa                                         >  1:2141841/1‑45 (MQ=255)
                            ggTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:1974748/1‑58 (MQ=255)
                            ggTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATg                          >  1:2146067/1‑60 (MQ=255)
                            ggTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATg                          >  1:2053105/1‑60 (MQ=255)
                            ggTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATg                          >  1:1072601/1‑60 (MQ=255)
                            ggTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATg                          >  1:302969/1‑60 (MQ=255)
                            ggTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGACGGTAACATCATGAAAGCATg                          >  1:868460/1‑60 (MQ=255)
                                 ccGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg           <  1:650028/70‑1 (MQ=255)
                                 ccGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg           <  1:756534/70‑1 (MQ=255)
                                 ccGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg           <  1:409806/70‑1 (MQ=255)
                                 ccGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg           <  1:3286901/70‑1 (MQ=255)
                                 ccGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg           <  1:2467215/70‑1 (MQ=255)
                                  cGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCgg           <  1:1657558/69‑1 (MQ=255)
                                          ctcCATGCCCGCCGGATATGAAAGCCGGTAACATCATGAAAGCa                            >  1:3174283/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATg                                  >  1:2277493/1‑38 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:2428184/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:683035/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:943690/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:577385/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:1662782/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:3171418/1‑44 (MQ=255)
                                          ctcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCa                            >  1:2633393/1‑44 (MQ=255)
                                           tcCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAg  <  1:498220/69‑1 (MQ=255)
                                             |                                                                  
AACTGTTTAACGGTGAAATCGACTGCCTGGTTTCCGACCACTCTCCATGCCCGCCGGAAATGAAAGCCGGTAACATCATGAAAGCATGGGGCGGTATCGCCGGTCTGCAAAG  >  W3110S.gb/539278‑539389

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: