Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,245,658 G→A 47.6% S63S (AGC→AGT xylE ← D‑xylose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,245,6580GA47.6% ‑1.7 / 26.3 21S63S (AGC→AGTxylED‑xylose transporter
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base A (8/2);  total (17/4)
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.17e-01

AATCACGACGACCGAAGCGGTTACTGCAATAACCACCGAGGGCACCGCCGATGATGCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTTGTGGAGCAACAA  >  W3110S.gb/4245588‑4245722
                                                                      |                                                                
aaTCACGACGACCGAAGCGGTTACTGCAATAACCACCGAGGGCACCGCCGATGATGCAACCAATCAGAgcg                                                                  <  1:2036899/71‑1 (MQ=255)
    acgacgACCGAAGCGGTTACTGCAATAACCACCGAGGGCACCGCCGATGATGCAACCAATCAGAGCGCTgg                                                              >  1:2023503/1‑71 (MQ=255)
                     tACTGCAATAACCACCGAGGGCACCGCCGATGATGCAACCAATCAGAGCACTGGCCACGCAAAACCCTaa                                              >  1:1423002/1‑70 (MQ=255)
                       cTGCAATAACCACCGAGGGCACCGCCGATGATGCAACCAATCAGAGCGCTGGCCACGCaaa                                                     >  1:1458367/1‑61 (MQ=255)
                               accaccGAGGGCACCGCCGATGATGCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTgg                                   >  1:2043683/1‑71 (MQ=255)
                               accaccGAGGGCACCGCCGATGATGCAACCAATCAGAGCACTGGCCACGCAAAACCCTAACAGGGAGTTgg                                   >  1:1826697/1‑71 (MQ=255)
                                     gAGGGCACCGCCGATGATGCAACCAATCAGAGCACTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCg                              >  1:1960042/1‑70 (MQ=255)
                                        ggCACCGCCGATGATGCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGAtt                          <  1:1242956/71‑1 (MQ=255)
                                         gCACCGCCGATGATGCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGAtt                          <  1:1445016/70‑1 (MQ=255)
                                              gccgATGATGCAACCAATCAGAGCACTGGCCACGCaaaa                                                    >  1:1449405/1‑39 (MQ=255)
                                              gccgATGATGCAACCAATCAGAGCACTGGCCACGCaaaa                                                    >  1:573880/1‑39 (MQ=255)
                                              gccgATGATGCAACCAATCAGAGCACTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACtt                     >  1:1614832/1‑70 (MQ=255)
                                              gccgATGATGCAACCAATCAGAGCACTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACtt                     >  1:1768338/1‑70 (MQ=255)
                                                    gatgCAACCAATCAGAGCGCTGGCCACGCAAAACCCTaa                                              >  1:1656036/1‑39 (MQ=255)
                                                     atgCAACCAATCAGAGCACTGGCCACGCAAAACCCTAACAggg                                         <  1:1684522/43‑1 (MQ=255)
                                                     atgCAACCAATCAGAGCACTGGCCACGCAAAACCCTAACAggg                                         <  1:1494447/43‑1 (MQ=255)
                                                       gCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTtgtg           >  1:1544704/1‑71 (MQ=255)
                                                       gCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTtgtg           >  1:583280/1‑71 (MQ=255)
                                                             aaTCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTTGAGGAGca      >  1:2668785/1‑70 (MQ=255)
                                                                 agagCGCTGGCCACTCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTTGTGGAGcaaca   >  1:1161019/1‑69 (MQ=255)
                                                                 agagCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTTGTGGAGcaacaa  >  1:85890/1‑70 (MQ=255)
                                                                 agagCACTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTTGTGGAGcaaca   >  1:1730159/1‑69 (MQ=255)
                                                                      |                                                                
AATCACGACGACCGAAGCGGTTACTGCAATAACCACCGAGGGCACCGCCGATGATGCAACCAATCAGAGCGCTGGCCACGCAAAACCCTAACAGGGAGTTGGCAGCGGATTCACTTAAGTTTTGTGGAGCAACAA  >  W3110S.gb/4245588‑4245722

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: