Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,250,252 | 0 | T | C | 66.7% | 28.5 / 25.0 | 30 | F179S (TTC→TCC) | yeiH | conserved inner membrane protein YeiH |
Reads supporting (aligned to +/- strand): ref base T (7/3); new base C (6/14); total (13/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.62e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTACTGTG > NC_000913/2250174‑2250331 | aGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGAtttt < 1:368470/79‑1 (MQ=255) aGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTccc > 1:368506/1‑81 (MQ=255) tAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTAcc < 2:368471/81‑1 (MQ=255) gAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGc > 1:368508/1‑81 (MQ=255) aaaGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGc < 2:368472/80‑1 (MQ=255) aaaGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCg < 1:368473/81‑1 (MQ=255) gCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGata < 2:368474/81‑1 (MQ=255) ccAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGc < 1:368475/81‑1 (MQ=255) agtaaagtaaCCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCt > 1:368509/1‑80 (MQ=255) taaagtaaCCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGAt > 2:368520/1‑81 (MQ=255) gtaaCCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGtct > 2:368524/1‑81 (MQ=255) gtaaCCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGtct < 2:368476/81‑1 (MQ=255) taaCCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGtctc < 2:368477/81‑1 (MQ=255) aaacGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGtctc > 2:368516/4‑80 (MQ=255) tGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGt > 1:368513/1‑81 (MQ=255) ggCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGtt > 1:368512/1‑81 (MQ=255) ggCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGtt > 2:368517/1‑81 (MQ=255) gCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGttt < 1:368478/81‑1 (MQ=255) ttGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAg < 1:368479/79‑1 (MQ=255) gCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCg > 2:368532/1‑81 (MQ=255) aCCGTTGTTATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGaa < 2:368480/81‑1 (MQ=255) ttgttATCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACtt < 1:368481/81‑1 (MQ=255) tCTTCGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGt < 1:368482/80‑1 (MQ=255) tCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTAATTCCGGAACCTTTCGGTacc > 2:368515/1‑78 (MQ=255) cGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTAt < 1:368484/78‑1 (MQ=255) cGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTAt < 1:368485/78‑1 (MQ=255) cGGGACCGTCGCGATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTAt < 2:368483/78‑1 (MQ=255) gTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGt > 2:368521/1‑81 (MQ=255) gcgATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCt < 2:368486/81‑1 (MQ=255) gATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTAc > 2:368531/1‑81 (MQ=255) gATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTAc < 2:368487/81‑1 (MQ=255) gATTTCCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTAc > 1:368514/1‑81 (MQ=255) ttCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTACtgtg > 2:368540/1‑81 (MQ=255) ttCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTACtgt < 2:368488/80‑1 (MQ=255) | AGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTTCGGTATCTACATTGGTTCTACTGTG > NC_000913/2250174‑2250331 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |