| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_004302 | 2,843 | 0 | C | G | 46.7% | 2.7 / 29.0 | 40 | intergenic (‑86/‑144) | rro/tec | Repressor protein/tec |
| Reads supporting (aligned to +/- strand): ref base C (16/0); new base G (14/0); total (40/0) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCACAAACAACAAAACTTGTCAATTTGGAAAG‑TTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGACGGTA > NC_004302/2769‑2912 | aGCACAAACAACAAAACTTGTCAATTTGGAAAG‑TTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑c > 1:5591819/1‑75 (MQ=255) tatacggtGTCAATTTGGAAAG‑TTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑aaa > 1:1999334/8‑75 (MQ=255) gcaattttttGGAAAG‑TTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCg > 1:6508547/8‑75 (MQ=255) caaggtGGAACG‑TTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACa > 1:17406969/6‑75 (MQ=255) attGGAAAGTTTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAAt > 1:16044243/2‑75 (MQ=255) tgaaatgaaCTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAATCATCGAg > 1:5540552/8‑75 (MQ=255) gtaaacgtAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGcggacgga > 1:16654241/8‑75 (MQ=255) tgattttAATGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGaa > 1:6979524/6‑75 (MQ=255) atgtaatAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGaa > 1:2345348/7‑75 (MQ=255) tcattaaaGGTTATAATTAATA‑CATAAAGTC‑AAACAACCGAACAATCATGGAGCATTACAGTAAGGCGGACGGaa > 1:15429414/7‑75 (MQ=255) gcattaaaGGTTATAATTAATT‑CATAAAGTC‑TAACAAGCGAACAATCATGGACAATTACAGTCCGGcggacggac > 1:14622910/7‑74 (MQ=255) gcggaattaact‑TATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAAt > 1:16319749/5‑75 (MQ=255) tgtggacaatt‑TATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATg > 1:14789891/8‑75 (MQ=255) tcttataaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATg > 1:16373733/8‑75 (MQ=255) tccaattaatt‑‑ATAAAGTC‑AAACAAGCGAACAATCATGGAGCACT‑CAGTACGGCGGACGGAACAGGCTCAAATGa > 1:7386401/4‑75 (MQ=255) tcaccacaatt‑TATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATg > 1:1649463/8‑75 (MQ=255) tataGTTAATT‑‑ATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGa > 1:9357105/1‑75 (MQ=255) gtatatgaatt‑TATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATg > 1:11932342/8‑75 (MQ=255) cggttttaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATg > 1:3725138/6‑75 (MQ=255) attggaaaatt‑TATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATg > 1:13326218/8‑75 (MQ=255) atgacaaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGa > 1:7474474/7‑75 (MQ=255) agtgttaatt‑‑ATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:9091376/5‑75 (MQ=255) aggcggaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGa > 1:16257530/7‑75 (MQ=255) aggcggaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGa > 1:6363290/7‑75 (MQ=255) tgggtaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:14048434/5‑75 (MQ=255) tctcaaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:7065910/6‑75 (MQ=255) taattaacT‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:4419612/1‑75 (MQ=255) gaattagtt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:13495421/2‑75 (MQ=255) gaagtaatt‑aatAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:11091432/5‑75 (MQ=255) ctgtaaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:12100135/6‑75 (MQ=255) cccagaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:11314658/6‑75 (MQ=255) cacttaaat‑CATAAAGTC‑AAACAAGCGAACAATCATCGATCAAT‑CAGAACGGCTGACCTAACAGGCTCAAATGAc > 1:12472795/4‑75 (MQ=255) agtcgaatt‑AATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:4461086/6‑75 (MQ=255) agaggtttt‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:4236709/8‑75 (MQ=255) agacgaatt‑TATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:11406523/6‑75 (MQ=255) aatttaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:10648665/4‑75 (MQ=255) aatttaatt‑GATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:11231399/4‑75 (MQ=255) aatgtagtt‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:15452635/8‑75 (MQ=255) aatgatatt‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:7688429/7‑75 (MQ=255) aaagaaatt‑GATAAAGTCAAAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGa > 1:3760300/6‑75 (MQ=255) tgaaatgt‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGACg > 1:15558393/8‑75 (MQ=255) gtttaattgCATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGAc > 1:1259052/3‑75 (MQ=255) ggctcaacATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGACgg > 1:10238359/8‑75 (MQ=255) ggaagacATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGACGGt > 1:1724182/7‑75 (MQ=255) tacggcATATAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGACGGTa > 1:17701517/6‑75 (MQ=255) | AGCACAAACAACAAAACTTGTCAATTTGGAAAG‑TTTTAACTTGACAAATTGGATAGTTAAGGTTATAATTAATT‑CATAAAGTC‑AAACAAGCGAACAATCATGGAGCATT‑CAGTACGGCGGACGGAACAGGCTCAAATGACGGTA > NC_004302/2769‑2912 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |