Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,769,471 | 1 | . | G | 51.7% | 43.2 / 48.0 | 29 | coding (398/1113 nt) | ydiK | UPF0118 family inner membrane protein |
Reads supporting (aligned to +/- strand): ref base . (5/9); new base G (8/7); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGC > NC_000913/1769339‑1769597 | tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑gggg > 1:1221528/1‑137 (MQ=255) tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑gggg < 2:1221528/137‑1 (MQ=255) aaTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTc < 2:473023/139‑1 (MQ=255) gCCATTCCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGCATGCAGGCTGGAACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCt < 1:1006916/139‑1 (MQ=255) tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGt < 2:1605482/127‑1 (MQ=255) tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGt > 1:1605482/1‑127 (MQ=255) cAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGc < 1:479030/139‑1 (MQ=255) aGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca < 1:611545/139‑1 (MQ=255) aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca > 2:730677/1‑139 (MQ=255) tGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGc < 1:153387/139‑1 (MQ=255) cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTc > 1:337463/1‑96 (MQ=255) cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg > 1:1254551/1‑139 (MQ=255) cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACCTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTc < 2:337463/96‑1 (MQ=255) ccATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg < 1:679092/139‑1 (MQ=255) ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgc > 1:299129/1‑139 (MQ=255) tCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg < 2:1254551/139‑1 (MQ=255) ggCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTg > 1:1259763/1‑139 (MQ=255) gCTGGCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc > 1:700172/1‑75 (MQ=255) gCTGGCACAACTTGCTGGATAT‑GGGGGGGTACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacc < 2:700172/75‑1 (MQ=255) cTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtctc > 1:558363/1‑139 (MQ=255) ggCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt < 1:1452875/116‑1 (MQ=255) ggCACAACTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt > 2:1452875/1‑116 (MQ=255) ggCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCtt > 2:294227/1‑139 (MQ=255) aCTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc < 1:294227/139‑1 (MQ=255) cTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt > 2:1126011/1‑102 (MQ=255) cTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt < 1:1126011/102‑1 (MQ=255) cTTGCTGGATATGGGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGc < 2:299129/139‑1 (MQ=255) cTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTg < 1:485707/53‑1 (MQ=255) cTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTg > 2:485707/1‑53 (MQ=255) cTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg > 2:187783/1‑133 (MQ=255) cTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg < 1:187783/133‑1 (MQ=255) | TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑GGGGGGG‑ACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGC > NC_000913/1769339‑1769597 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |