Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,428,262:1 +G 100% coding (735/855 nt) ulaE → L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,428,2621.G81.8% 57.7 / 6.6 22coding (735/855 nt)ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (10/8);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GCCGTTTGGCGAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGC  >  W3110S.gb/4428196‑4428321
                                                                   |                                                           
gCCGTTTGGCGAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCgg                                                           >  1:600047/1‑69 (MQ=255)
     ttGGCGAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGt                                                      <  1:896895/69‑1 (MQ=255)
          gAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACTGAGTGGCTATTGCGGGGCCGTACCt                                                  <  1:3422130/69‑1 (MQ=255)
          gAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  <  1:1770919/69‑1 (MQ=255)
          gAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  <  1:2525841/69‑1 (MQ=255)
                   gTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGAt                                         >  1:1891476/1‑69 (MQ=255)
                       aTTTCGGACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGc                                                              >  1:472795/1‑44 (MQ=255)
                            gAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGa                                <  1:1272070/68‑1 (MQ=255)
                             aaCGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGaa                               >  1:656110/1‑69 (MQ=255)
                               cGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  >  1:3092848/1‑48 (MQ=255)
                               cGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  >  1:867537/1‑48 (MQ=255)
                               cGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  >  1:3146197/1‑48 (MQ=255)
                               cGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  >  1:1131119/1‑48 (MQ=255)
                               cGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGATGTgg                                     <  1:405002/61‑1 (MQ=255)
                                   gttTCGAAACGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCt                                                  <  1:1231977/44‑1 (MQ=38)
                                      tCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAg                     <  1:3639599/69‑1 (MQ=255)
                                            cGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGAccc                 <  1:2637992/68‑1 (MQ=255)
                                            cGCTCAAACAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGATCTGGAGCGAAACGGCGGAAGAccc                 <  1:1550100/68‑1 (MQ=255)
                                             gCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCcggc              >  1:256109/1‑69 (MQ=255)
                                              cTCAAACAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGAAAc                             >  1:2885949/1‑54 (MQ=38)
                                                   aCAGAGTGGCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAg         <  1:663039/69‑1 (MQ=255)
                                                        gTGGCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAGTg    >  1:3709576/1‑69 (MQ=255)
                                                        gTGGCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAGTg    >  1:1595906/1‑69 (MQ=255)
                                                           gCTATTGCGGGGCCGTACCTGATTGAGATGTGGAGCGCAACGGCGGAAGACCCGGCGGCAGAAGTGGc  >  1:1322617/1‑68 (MQ=255)
                                                                   |                                                           
GCCGTTTGGCGAAGGTGTAGTGGATTTCGAACGTTGTTTCGAAACGCTCAAACAGAGTGGCTATTGCGGGCCGTACCTGATTGAGATGTGGAGCGAAACGGCGGAAGACCCGGCGGCAGAAGTGGC  >  W3110S.gb/4428196‑4428321

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: