Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,171,889 C→A 42.8% A21A (GCG→GCT)  ycfS ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,171,8890CA42.8% 9.9 / 9.3 14A21A (GCG→GCT) ycfSconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base A (2/4);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATC  >  W3110S.gb/1171842‑1171959
                                               |                                                                      
cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaa                                                  >  1:4146432/1‑70 (MQ=255)
   gCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaaa                                              <  1:912742/71‑1 (MQ=255)
   gCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGAGCCACGGCAGCGGCGAACGTaaaaaa                                              <  1:3963763/71‑1 (MQ=255)
    cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGAGCCACGGCAGCGGCGAACGTAAAAAAc                                             <  1:2840252/71‑1 (MQ=255)
    cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGAGCCACGGCAGCGGCGAACGTAAAAAAc                                             <  1:445161/71‑1 (MQ=255)
      ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGAGCCACGGCAGCGGCGAACGTAAAAAACGt                                           >  1:3304544/1‑71 (MQ=255)
       gcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACg                                            <  1:399710/69‑1 (MQ=255)
                       gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGt                                           >  1:4452516/1‑54 (MQ=255)
                       gtTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAAc                          >  1:2321709/1‑71 (MQ=255)
                                 gCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGa                 <  1:1411007/70‑1 (MQ=255)
                                 gCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGa                 <  1:798528/70‑1 (MQ=255)
                                     gTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcatc            <  1:3619545/71‑1 (MQ=255)
                                     gTAGCGCCAGAGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcat             >  1:3346289/1‑70 (MQ=255)
                                     gTAGCGCCAGAGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                <  1:1508982/67‑1 (MQ=255)
                                       agcgccagcgccACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAc          >  1:3849665/1‑71 (MQ=255)
                                       agcgccagcgccACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAc          >  1:1161337/1‑71 (MQ=255)
                                            cagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatca              <  1:416683/62‑1 (MQ=255)
                                               cgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATc  <  1:3702816/71‑1 (MQ=255)
                                               |                                                                      
CCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATC  >  W3110S.gb/1171842‑1171959

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑