Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 408,831 | 0 | C | T | 54.5% | 0.9 / 25.6 | 22 | pseudogene (223/342 nt) | ykiA | DUF2773 domain‑containing protein YkiA |
Reads supporting (aligned to +/- strand): ref base C (8/2); new base T (6/6); total (14/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.04e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAACGCG > NC_000913/408753‑408904 | gAGCGTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCca > 1:79114/1‑81 (MQ=255) gCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATTcatc < 1:79094/81‑1 (MQ=255) ggATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTAct > 2:79135/1‑81 (MQ=255) tCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATTCATCATTActct > 1:79126/1‑80 (MQ=255) tCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTActct > 2:79122/1‑80 (MQ=255) gCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGaa < 2:79095/81‑1 (MQ=255) tggctgCCGATGTGAAAACGGCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGt > 1:79130/1‑81 (MQ=255) gctgCCGATGTGAAAACGGCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGtt < 2:79096/80‑1 (MQ=255) gctgCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGtt < 1:79097/80‑1 (MQ=255) gctgCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGt < 1:79098/79‑1 (MQ=255) ccGATGTGAAAACGGCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGt > 1:79136/1‑80 (MQ=255) aaaCGGCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTg > 2:79137/1‑80 (MQ=255) aCGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGaa > 1:79142/1‑80 (MQ=255) ggCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAAc > 2:79138/1‑81 (MQ=255) gCGTGGTTAAAAGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAAcc < 2:79099/81‑1 (MQ=255) ttAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGc > 2:79140/1‑80 (MQ=255) aaaaGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCTCTGTTACGTGATTTAGTGAAAACCGGGGCAAc > 1:79127/1‑81 (MQ=255) aaaGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAAcg < 2:79101/81‑1 (MQ=255) aaaGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAAc < 1:79100/80‑1 (MQ=255) aaGAGGATTCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAAcgc > 1:79132/1‑81 (MQ=255) aaGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAAcgc > 2:79131/1‑81 (MQ=255) gagGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAAcgcg > 1:79133/1‑80 (MQ=255) | GAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAGTGAAAACCGGGGCAACGCG > NC_000913/408753‑408904 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |