Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,129,117 C→T 45.2% 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,1170CT45.2% 7.1 / 38.3 31A71T (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 (14/3);  new base T (9/5);  total (23/8)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCAT  >  W3110S.gb/2129050‑2129185
                                                                   |                                                                    
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:2379707/71‑1 (MQ=255)
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:1561922/71‑1 (MQ=255)
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:612337/71‑1 (MQ=255)
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:2219334/71‑1 (MQ=255)
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:2133653/71‑1 (MQ=255)
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:3606058/71‑1 (MQ=255)
gcgTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCCGATAGGTGTTGTTGGtaac                                                                   <  1:1113896/71‑1 (MQ=255)
  gTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGtaacaa                                                                 <  1:3389578/71‑1 (MQ=255)
   tGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGtaacaa                                                                 <  1:2369488/70‑1 (MQ=255)
         atgTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGGCGGATAGGTGTTGTTGGTAACAATgccgcc                                                          >  1:850161/1‑71 (MQ=255)
          tgTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATgccgcc                                                          >  1:2718173/1‑70 (MQ=255)
          tgTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATgccgcc                                                          >  1:2928515/1‑70 (MQ=255)
             tGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACt                                                      >  1:523457/1‑71 (MQ=255)
                     atcatcATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:513377/50‑1 (MQ=255)
                     atcatcATGTTCTCGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaac                                                                   <  1:1014871/50‑1 (MQ=255)
                        atcatGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGtaa                                                                    >  1:1116498/1‑46 (MQ=255)
                          catGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTcgccgccgc                                           >  1:93287/1‑69 (MQ=255)
                              ttCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACt                                                      <  1:3184887/54‑1 (MQ=255)
                              ttCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACt                                                      <  1:949304/54‑1 (MQ=255)
                                cTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGcaacaa                                                                 >  1:3533701/1‑41 (MQ=255)
                                cTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAg                                    <  1:2542620/70‑1 (MQ=255)
                                cTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAg                                    <  1:1986891/70‑1 (MQ=255)
                                    tAGATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTg                               >  1:742156/1‑71 (MQ=255)
                                    tAGATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTg                               >  1:1986473/1‑71 (MQ=255)
                                     aGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATgccg                                                            >  1:141176/1‑41 (MQ=255)
                                      gATGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGa                                       <  1:1521060/61‑1 (MQ=255)
                                      gATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGt                             >  1:849849/1‑71 (MQ=255)
                                      gATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGt                             >  1:3253956/1‑71 (MQ=255)
                                      gATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGt                             >  1:3344754/1‑71 (MQ=255)
                                      gATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCACACTTTCGCCGCCGCCAGATAGACCTGGt                             >  1:2569303/1‑71 (MQ=255)
                                        tGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGAt                                      >  1:268726/1‑60 (MQ=255)
                                        tGAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATa                                     >  1:1986408/1‑61 (MQ=255)
                                         gAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCaa                          >  1:2055433/1‑71 (MQ=255)
                                         gAAATCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCaa                          >  1:3163973/1‑71 (MQ=255)
                                         gAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCaa                          >  1:306758/1‑71 (MQ=255)
                                          aaaTCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGa                                       >  1:3717540/1‑57 (MQ=255)
                                          aaaTCCGCCGGATAGGTGTTGTTGGTAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCaa                          >  1:1574068/1‑70 (MQ=255)
                                             tCCGCCGGATAGGTGTTGTTGGCAACAATgccgcc                                                          >  1:1871653/1‑35 (MQ=255)
                                             tCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGa                                       >  1:3345072/1‑54 (MQ=255)
                                                      tAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTgg              >  1:2660029/1‑70 (MQ=255)
                                                              gttgGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAAt                         >  1:335653/1‑51 (MQ=255)
                                                               ttgGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAAt     <  1:3239781/70‑1 (MQ=255)
                                                               ttgGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAAt     <  1:1430936/70‑1 (MQ=255)
                                                                  gCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCAt  <  1:2805923/70‑1 (MQ=255)
                                                                   |                                                                    
GCGTGAATGATGTTGCTCTCAATCATCATGTTCTGGTAGATGAAATCCGCCGGATAGGTGTTGTTGGCAACAATGCCGCCCACTTTCGCCGCCGCCAGATAGACCTGGTCAATACGTTCGCTGGCAAAGAAATCAT  >  W3110S.gb/2129050‑2129185

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: