Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,247,148 | 0 | C | T | 31.6% | 17.4 / 12.3 | 19 | A77T (GCT→ACT) | treA | periplasmic trehalase |
Reads supporting (aligned to +/- strand): ref base C (6/7); new base T (4/2); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTT > NC_000913/1247069‑1247214 | tcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc > 1:221577/1‑80 (MQ=255) cGGCAGGGCGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCCGGTTTTGCTGCATCCGATAATCAGTAAGGat > 2:221587/1‑81 (MQ=255) aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatcatc < 2:221562/81‑1 (MQ=255) aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatcat > 1:221582/1‑80 (MQ=255) gggTGAAATTGACGTTAACGGAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGTAAGGatcat > 2:221578/1‑79 (MQ=255) gTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGc > 2:221579/1‑81 (MQ=255) aaTTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGAt > 2:221590/1‑81 (MQ=255) tGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGTAAGGATCATCAGCGGATCGc > 2:221586/1‑81 (MQ=255) gACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGTAAGGATCATCAGCGGATCGc < 2:221563/80‑1 (MQ=255) aCGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggc < 1:221564/81‑1 (MQ=255) aTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGTAAGGATCATCAGCGGATCGCTGTTcggcacgg > 1:221593/1‑80 (MQ=255) cgcAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCg < 2:221565/81‑1 (MQ=255) cGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGtttttt < 2:221566/81‑1 (MQ=255) ctGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTc > 1:221591/1‑81 (MQ=255) ggTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg > 1:221584/1‑81 (MQ=255) gTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg < 2:221567/80‑1 (MQ=255) tgctgcATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaa > 2:221581/1‑80 (MQ=255) gctgcATCCGATAATCAGTAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa < 1:221568/80‑1 (MQ=255) gctgcATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa < 1:221569/80‑1 (MQ=255) cATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGttt < 1:221570/81‑1 (MQ=255) | TCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTT > NC_000913/1247069‑1247214 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |