Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,373,368 G→A R480H (CGT→CAT)  ycjM → predicted glucosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,373,3680GA88.0% 64.8 / 2.4 25R480H (CGT→CAT) ycjMpredicted glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (15/7);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 5.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTTCGTCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGATACGATTAATTCATCCGT  >  W3110S.gb/1373303‑1373431
                                                                 |                                                               
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:916633/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:513406/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:1447929/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:3656848/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:3532193/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:2931960/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:1251720/1‑69 (MQ=255)
aCTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTtcatcg                                                              >  1:2332462/1‑69 (MQ=255)
          gAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAAc                                                    >  1:284392/1‑69 (MQ=255)
          gAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAAc                                                    >  1:1597260/1‑69 (MQ=255)
            aGCTACATTAAGTCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTTCGTCGCAGCCATAACGAg                                                 <  1:3398267/70‑1 (MQ=255)
            aGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTTCGTCGCAGCCATAACGAg                                                 <  1:282636/70‑1 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGAtaataa                                >  1:3900540/1‑70 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:1371557/1‑71 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:2183242/1‑71 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:2891330/1‑71 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:3637745/1‑71 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:3535270/1‑71 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:1610367/1‑71 (MQ=255)
                             cGGTATATCATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:3530386/1‑71 (MQ=255)
                             cGGTATATAATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:4013785/1‑71 (MQ=255)
                             cGGTATAACATGAGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATaataat                               >  1:169738/1‑71 (MQ=255)
                                         aGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGAt                   <  1:3499346/71‑1 (MQ=255)
                                         aGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGAt                   <  1:3358952/71‑1 (MQ=255)
                                         aGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGAt                   <  1:1645221/71‑1 (MQ=255)
                                         aGTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGAt                   <  1:2531887/71‑1 (MQ=255)
                                          gTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGAt                   <  1:2081213/70‑1 (MQ=255)
                                          gTTGTCGCGTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGAt                   <  1:1422320/70‑1 (MQ=255)
                                                  gTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTATCATTGATACGATTAAt          >  1:3908580/1‑71 (MQ=255)
                                                  gTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGATACGATTAAt          >  1:1909216/1‑71 (MQ=255)
                                                  gTTTAATTACACTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGATACGATTAAt          >  1:1614672/1‑71 (MQ=255)
                                                           cacTTCGTCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGATACGATTAATTCATCCGt  <  1:897042/70‑1 (MQ=255)
                                                           cacTTCATCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGATACGATTAATTCATCCGt  <  1:1926976/70‑1 (MQ=255)
                                                                 |                                                               
ACTGAACGATGAAGCTACATTAAGGCATGCGGTATATCATGAGTTGTCGCGTTTAATTACACTTCGTCGCAGCCATAACGAGTTTCATCCGGATAATAATTTTACCATTGATACGATTAATTCATCCGT  >  W3110S.gb/1373303‑1373431

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: