| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 953,088 | 1 | . | C | 20.0% | 48.9 / 14.1 | 20 | coding (467/2283 nt) | pflB | formate C‑acetyltransferase 1, anaerobic, pyruvate formate‑lyase 1 |
| Reads supporting (aligned to +/- strand): ref base . (7/9); new base C (2/2); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAACCTTCGATCATTTTGATACCACCGAACGGGAT > NC_000913/952953‑953224 | aaCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCg > 1:239640/1‑139 (MQ=255)aaCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCg < 1:97693/139‑1 (MQ=255) tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGcc < 2:498519/139‑1 (MQ=255) aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtt < 2:22778/139‑1 (MQ=255) cGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt < 2:180562/139‑1 (MQ=255) tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa > 1:495781/1‑139 (MQ=255) tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGATCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa > 1:233853/1‑139 (MQ=255) cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGccc > 1:155073/1‑118 (MQ=255) cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGccc < 2:155073/118‑1 (MQ=255) cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt > 1:60543/1‑139 (MQ=255) aaCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGccc > 1:204496/1‑115 (MQ=255) aaCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGccc < 2:204496/115‑1 (MQ=255) tgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat < 2:239640/139‑1 (MQ=255) atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc > 1:62341/1‑109 (MQ=255) atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc < 2:62341/109‑1 (MQ=255) ggccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGAt < 2:495781/139‑1 (MQ=255) ggccacggccaTATGCATCTGGCAGACCGATCAGAACACCAGATTTACGGCAACGCAGGATGTCCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGAt < 2:233853/139‑1 (MQ=255) cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg < 2:251701/79‑1 (MQ=255) cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg > 1:251701/1‑79 (MQ=255) tGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAACCTTCGATCATTTTGATACCACCGAACGGGAt > 2:449447/1‑139 (MQ=255) | AACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTGTACGCTTTGCAGGAACCTTCGATCATTTTGATACCACCGAACGGGAT > NC_000913/952953‑953224 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |