Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,129,117 C→T 40.0% A71T (GCC→ACC)  fcl ← bifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase and GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,129,1170CT40.0% 13.9 / 25.3 25A71T (GCC→ACC) fclbifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase and GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase
Reads supporting (aligned to +/- strand):  ref base C (4/11);  new base T (3/7);  total (7/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAA  >  W3110S.gb/2129052‑2129180
                                                                    |                                                               
gTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:759125/69‑1 (MQ=255)
gTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:1836921/69‑1 (MQ=255)
gTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:1956605/69‑1 (MQ=255)
gTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:3216975/69‑1 (MQ=255)
gTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:2863084/69‑1 (MQ=255)
          ttGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:2948104/59‑1 (MQ=255)
          ttGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCGccc                                                    >  1:1249908/1‑69 (MQ=255)
           tGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaaca                                                             <  1:3244812/59‑1 (MQ=255)
            gCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGtaac                                                              <  1:7592/57‑1 (MQ=255)
                 cAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCGCCCACTTTcg                                             >  1:3433213/1‑69 (MQ=255)
                  aaTCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgc                                            >  1:1704390/1‑69 (MQ=255)
                  aaTCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑AGTTGTTGGCAACAATGCCGCCCACTTTcgc                                            >  1:594403/1‑69 (MQ=255)
                    tcatcaTGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgc                                            <  1:824526/67‑1 (MQ=255)
                    tcatcaTGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgc                                            <  1:3318131/67‑1 (MQ=255)
                     catcatGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGcaac                                                              <  1:1280952/48‑1 (MQ=255)
                     catcatGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGcaac                                                              <  1:3488772/48‑1 (MQ=255)
                           gTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCa                                    <  1:1697923/68‑1 (MQ=255)
                           gTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCa                                    <  1:3230505/68‑1 (MQ=255)
                                    gATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCa                                    <  1:1517241/59‑1 (MQ=255)
                                    gATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCa                                    <  1:237027/59‑1 (MQ=255)
                                    gATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCa                                    >  1:2455243/1‑59 (MQ=255)
                                       gAAATCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTc                       <  1:3232522/69‑1 (MQ=255)
                                          aTCCGCCGGATAGGTGTTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTc                       <  1:2929853/69‑1 (MQ=255)
                                          aTCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCaa                     <  1:3182488/68‑1 (MQ=255)
                                           tCCGCCGGATAGG‑‑‑TGTTGTTGGTAACAATGCCTCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATa                   <  1:1592089/69‑1 (MQ=255)
                                              gccgGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACttt                                               <  1:3021374/38‑1 (MQ=255)
                                               ccgGATAGG‑‑‑TGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGtt               <  1:3224052/69‑1 (MQ=255)
                                                    tAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgccgccgc                                      <  1:309146/41‑1 (MQ=255)
                                                    tAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgccgccgc                                      <  1:2127516/41‑1 (MQ=255)
                                                    tAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgccgccgc                                      <  1:2677658/41‑1 (MQ=255)
                                                    tAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTcgccgccgc                                      <  1:3538098/41‑1 (MQ=255)
                                                    tAGG‑‑‑TGTTGTTGGCAACAAAGCCGCCCACTTTcgccgccgc                                      <  1:2685267/41‑1 (MQ=255)
                                                               gttgGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAg                               <  1:2152779/40‑1 (MQ=255)
                                                               gttgGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGaa  >  1:1428098/1‑69 (MQ=255)
                                                               gttgGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATAc                  >  1:12822/1‑53 (MQ=255)
                                                                    |                                                               
GTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGG‑‑‑TGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAA  >  W3110S.gb/2129052‑2129180

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: