Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,286,501 C→A 48.2% A261S (GCG→TCG)  yejK ← nucleotide associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,286,5010CA48.2% 20.7 / 29.6 29A261S (GCG→TCG) yejKnucleotide associated protein
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base A (6/8);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

CGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCG  >  W3110S.gb/2286433‑2286562
                                                                    |                                                             
cGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCa                                                              >  1:1048443/1‑70 (MQ=255)
    aaCTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTAAGCCGCAAAc                                                           >  1:1618365/1‑69 (MQ=255)
            cAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGc                                                 <  1:1956215/71‑1 (MQ=255)
              gtgtGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTcc                                                        >  1:776328/1‑62 (MQ=255)
                 tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACt                                            <  1:439685/71‑1 (MQ=255)
                      cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCg                                         <  1:1430597/69‑1 (MQ=255)
                      cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCt                                       >  1:819273/1‑71 (MQ=255)
                      cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCt                                       <  1:2338627/71‑1 (MQ=255)
                      cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCt                                       >  1:816179/1‑71 (MQ=255)
                       ggTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                       <  1:563804/70‑1 (MQ=255)
                        gTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTa                                     >  1:2573534/1‑71 (MQ=255)
                         tCTGCTGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGaaaa                                                            >  1:1269437/1‑47 (MQ=39)
                         tCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTAc                                    <  1:989100/71‑1 (MQ=255)
                         tCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTAc                                    <  1:193598/71‑1 (MQ=255)
                           tGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACtt                                           >  1:1958043/1‑62 (MQ=255)
                             ccGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCt                                <  1:738519/71‑1 (MQ=255)
                                    gCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGc                          <  1:1723963/70‑1 (MQ=255)
                                     ctctctTCCAGTTCGTAGCCTTTTTCAGCCGaaaa                                                            >  1:504904/1‑35 (MQ=255)
                                        tctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCttt                     >  1:1469714/1‑71 (MQ=255)
                                         ctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAg                                                  >  1:68338/1‑41 (MQ=255)
                                                  cGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                       <  1:1147680/43‑1 (MQ=255)
                                                   gTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAtt          >  1:799720/1‑71 (MQ=255)
                                                      gCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt         >  1:912457/1‑69 (MQ=255)
                                                      gCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt         >  1:26066/1‑69 (MQ=255)
                                                       ccTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAg      <  1:984280/71‑1 (MQ=255)
                                                        cTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGct                         <  1:477845/51‑1 (MQ=255)
                                                        cTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGct                         <  1:928581/51‑1 (MQ=255)
                                                           tttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTc   >  1:185587/1‑70 (MQ=255)
                                                           tttCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTa                                     >  1:624439/1‑36 (MQ=255)
                                                           tttCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTc   >  1:2249842/1‑70 (MQ=255)
                                                            ttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCg  <  1:2359356/70‑1 (MQ=255)
                                                                    |                                                             
CGTCAACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCG  >  W3110S.gb/2286433‑2286562

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: