Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_017635 | 4,454,267 | 0 | A | T | 30.0% | 61.8 / 14.0 | 30 | intergenic (+48/+193) | rpoC/thiH | DNA‑directed RNA polymerase subunit beta'/2‑iminoacetate synthase ThiH |
Reads supporting (aligned to +/- strand): ref base A (10/11); new base T (9/0); total (19/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.06e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGAAAGGC > NC_017635/4454220‑4454315 | tCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATg < 1:1603441/50‑1 (MQ=255) tAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGgt‑g > 2:1687971/1‑50 (MQ=21) aaTCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTTTTGgt‑gt > 2:1306504/1‑50 (MQ=21) aTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGGTTTTTTTTTGGgg‑gt > 2:1507264/1‑50 (MQ=18) gCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTAtt > 1:1036778/1‑49 (MQ=255) cAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGGG‑GTttttta > 2:1086255/1‑48 (MQ=17) aaaTAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTTTTGGG‑GTTAtt > 2:1723074/1‑47 (MQ=18) aaTAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTTTTGGTGGTTAtta > 2:618131/1‑47 (MQ=17) aaTAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCaga < 2:1168698/50‑1 (MQ=255) aTAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCagag < 1:44486/50‑1 (MQ=255) gTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCa > 1:1842426/1‑49 (MQ=255) tAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCaaa > 1:1765038/1‑50 (MQ=255) tAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGGG‑GTTTTTTATAGGGCaaa > 2:301777/1‑50 (MQ=9) aaaaaCCCGCTTCGGC‑GGGTTTTTTTTTGGGG‑GTTTTTTAGAGGGCaaaa > 2:2004644/1‑49 (MQ=17) aaaaCCCGCTTCGGC‑GGGTTTTTTTTTTGGG‑GTTTTTTAGAGGGCaaaat > 2:1433399/1‑48 (MQ=17) aaaaCCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAAt < 1:125981/49‑1 (MQ=255) aaaaCCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAAt < 1:549977/49‑1 (MQ=255) cTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACgg < 2:924180/50‑1 (MQ=255) cggc‑ggGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGa < 2:1452896/48‑1 (MQ=255) cggc‑ggGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGa < 2:1569639/48‑1 (MQ=255) ggc‑ggGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAAtt > 1:1546826/1‑50 (MQ=255) gc‑ggGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAAttt < 2:763599/50‑1 (MQ=255) c‑ggGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAAtttt < 2:1671866/50‑1 (MQ=255) gggTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAAtttt > 2:1140060/1‑49 (MQ=255) ggTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAAttttt > 2:537503/1‑49 (MQ=255) ggTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTa < 1:431114/50‑1 (MQ=255) tttttATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGa > 1:1665900/1‑49 (MQ=255) ttttATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGa > 1:71951/1‑48 (MQ=255) ttATGGGT‑GTTATTCAGAGGGCGAATTCTTCACGGAATTTTTAAGAAAgg > 2:646731/1‑50 (MQ=255) tATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGAAAgg > 1:1777044/1‑49 (MQ=255) aTGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGAAAGGc < 1:991809/49‑1 (MQ=255) aTGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGAAAGGc < 2:1929188/49‑1 (MQ=255) | TCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGT‑GTTATTCAGAGGGCAAATTCTTCACGGAATTTTTAAGAAAGGC > NC_017635/4454220‑4454315 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |