Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,204,480 | 0 | T | C | 66.7% | 17.1 / 12.4 | 18 | Y290H (TAT→CAT) | yehI | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (6/0); new base C (0/12); total (6/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.39e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAC > W3110S.gb/2204441‑2204544 | cgcgAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGCATGCCCATAAAGAAGATTGCGGCAGTctgct < 1:1774919/71‑1 (MQ=255) cgcgAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGCATGCCCATAAAGAAGATTGCGGCAGTctgct < 1:3854164/71‑1 (MQ=255) tCCGCTCCTCGCGATTCATGCCCATAAAGAAGATTGCGGCAGTctgct < 1:3879617/48‑1 (MQ=255) ccGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGt > 1:797455/1‑71 (MQ=255) ccGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGt > 1:339299/1‑71 (MQ=255) ccGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGt > 1:3353544/1‑71 (MQ=255) ccGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGt > 1:2888099/1‑71 (MQ=255) ccGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGt > 1:2411351/1‑71 (MQ=255) ccGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGt > 1:1881327/1‑71 (MQ=255) cgcgCCTCGCGATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTc < 1:362499/71‑1 (MQ=255) cgcgCCTCGCGATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTc < 1:3752801/71‑1 (MQ=255) cgcgCCTCGCGATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTc < 1:3796581/71‑1 (MQ=255) ccTCGCGATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATc < 1:3440845/70‑1 (MQ=255) gcgATGCATTCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAc < 1:1199948/71‑1 (MQ=255) gcgATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAc < 1:3179258/71‑1 (MQ=255) gcgATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAc < 1:3532780/71‑1 (MQ=255) gcgATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAc < 1:3559597/71‑1 (MQ=255) gcgATGCATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAc < 1:3644741/71‑1 (MQ=255) | CGCGAACAAGGACTGGCAGCCATTCCGCGCCTCGCGATGTATGCCCATAAAGAAGATTGCGGCAGTCTGCTGGTACAAATTAACCATCCGCAAGTCATCCGCAC > W3110S.gb/2204441‑2204544 |
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: ‑ |