Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,051,614 T→C 25.9% G126G (GGT→GGC)  ugpC → glycerol‑3‑phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,051,6140TC25.9% 40.5 / 15.2 27G126G (GGT→GGC) ugpCglycerol‑3‑phosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (15/5);  new base C (4/3);  total (19/8)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TGGGCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCG  >  W3110S.gb/4051553‑4051678
                                                             |                                                                
tGGGCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCaaa                                                         >  1:3213859/1‑71 (MQ=255)
tGGGCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCa                                                           <  1:745826/69‑1 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGt                                                                  >  1:63930/1‑60 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGt                                                                  >  1:3981867/1‑60 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCt                                                                >  1:2052371/1‑62 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCaaa                                                         >  1:1825615/1‑69 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACg                                                       >  1:1355777/1‑71 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACg                                                       >  1:4160744/1‑71 (MQ=255)
  ggCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACg                                                       >  1:2168047/1‑71 (MQ=255)
   gCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGt                                                      >  1:139552/1‑71 (MQ=255)
   gCAAGCAGCAAATTCCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGt                                                      >  1:917971/1‑71 (MQ=255)
        cagcaAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTTGAGCTGGACGGCCTGCTCAAACGTCGcc                                                  >  1:906202/1‑70 (MQ=255)
         agcaAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGCCTGCTCAAACGt                                                      <  1:180916/65‑1 (MQ=255)
         agcaAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGCCTGCTCAAACGt                                                      <  1:4336398/65‑1 (MQ=255)
            aaaTTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCcgcgc                                             >  1:2756645/1‑71 (MQ=255)
                tGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGc                                         >  1:2283532/1‑71 (MQ=255)
                 gCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGCCTGCTCaaa                                                         >  1:1876984/1‑54 (MQ=255)
                  ccGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCt                                        <  1:1467966/70‑1 (MQ=255)
                            tAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCa                                                           <  1:3638281/41‑1 (MQ=255)
                             aaaGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGt                              >  1:630174/1‑69 (MQ=255)
                                            aTTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGc             <  1:3284836/71‑1 (MQ=255)
                                             ttCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCg            <  1:4249605/71‑1 (MQ=255)
                                             ttCTGGAGCTGGACGGCCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCgtgt                <  1:634984/67‑1 (MQ=255)
                                               cTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAg                           >  1:1071918/1‑54 (MQ=255)
                                               cTGGAGCTGGACGGCCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGAt          >  1:4543180/1‑71 (MQ=255)
                                               cTGGAGCTGGACGGCCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGAt          >  1:2690653/1‑71 (MQ=255)
                                                tGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATg         >  1:3484827/1‑71 (MQ=255)
                                                     cTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGAGTGGCGATGGGcc     >  1:1621518/1‑70 (MQ=255)
                                                       ggACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCcgcg  >  1:234012/1‑71 (MQ=255)
                                                             |                                                                
TGGGCAAGCAGCAAATTGCCGAGCGCGTTAAAGAAGCGGCGCGCATTCTGGAGCTGGACGGTCTGCTCAAACGTCGCCCGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCG  >  W3110S.gb/4051553‑4051678

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: ‑