| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB143 | 5,251,184 | 0 | C | A | 52.0% | 18.0 / 30.3 | 25 | intergenic (‑64/+80) | CJ019_04948/fabV | hypothetical protein/Enoyl‑[acyl‑carrier‑protein] reductase [NADH] |
| Reads supporting (aligned to +/- strand): ref base C (3/9); new base A (13/0); total (16/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGA‑GGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTG > ACB143/5251052‑5251309 | gccgccGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTcgcg < 2:206790/136‑1 (MQ=255) gccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTGGTAGCGGGt > 2:202215/1‑136 (MQ=255) gccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGGTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGt < 6:47802/136‑1 (MQ=255) cgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGa < 1:247933/136‑1 (MQ=255) atgcggTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCGTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAc < 6:356404/131‑1 (MQ=255) aGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAc < 8:193945/136‑1 (MQ=255) aGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAc < 3:307669/136‑1 (MQ=255) gCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGATCCCGCGGCGACAgggg > 5:362997/1‑136 (MQ=255) gCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGAACCCCCTCCTACAgggg > 4:284066/1‑136 (MQ=255) gCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGAACCCCCTCCTACAgggg > 3:273760/1‑136 (MQ=255) cgGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGCACCAGCTCATACAGGGGgag > 3:235171/1‑135 (MQ=255) gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCtgtg < 6:188854/136‑1 (MQ=255) gCGAAATGGGGTAGCGGGGTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCtgtg < 1:40800/136‑1 (MQ=255) aCTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTTAACCCCCGCCAACAGGGGGAGAACCGGTGGGCGCGGGGTTGCCCGCCAAGA‑GGCCggga > 4:97931/1‑134 (MQ=255) aCTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCGTAGCGGGTGGAACCCCTCCGACAGGGGGGGGACCAGTGGGGGCGGGGGTGCCCCCGCAGAGGGCCGGt > 6:139451/1‑136 (MQ=255) gcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGTACCCCCTCCTACAGGGGGGGAACCTGTGGGGGCGGGGTTGGCCGCGGAGA‑GGCCGGGACCGGCGTAc > 1:209425/1‑136 (MQ=255) gcATATCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGAACCCCCGCCGACAGGGGGAGTACCGGTGGGGGCGGGGTTGGCCGCGCAGA‑GGGCGGGACAGGcgttc > 5:155346/1‑132 (MQ=255) ccTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCGTAGCGGGTGAACCCCCGCCGACACGGGGGGAACCTGTGGGGGCGGGGGTTCCCGCGCAGA‑GGCCGGGACAGGCTTAACCCAtt > 8:39034/1‑136 (MQ=255) gttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGAACCCGCGCCTACAGGGGGCGAACCGGTGGGGGCGGGGTTGCCCGCGCAGA‑GGCCGGGTCAGGCTTTCCCGATTACaac > 2:304644/1‑133 (MQ=37) gttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGGGAAACCCCGCCTAAAGGGGGAGAACCGGTGGGGGCGGGGTTGGCCGCGAAGA‑GGGCGGGAAAGGCTTAACCCATTaaacc > 8:240877/1‑131 (MQ=255) cgcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTTGACCCCGAGCGAAAGGGGGGGAACCTGTGGGAGCGGGGGTTCCCGGGCAGA‑GGGCGGTTCAGGCTTACCCGATTTCAGCAcgccgc > 2:203543/1‑136 (MQ=255) gcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCCCTCCGACAGGGGGATAACCGGGGGGAGCGGGCTTGGCCGCGAAGA‑GGCCGGTTCACGCTTACCCCATTACACCTCGACGCa > 8:121852/1‑136 (MQ=37) cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAc > 3:188497/1‑50 (MQ=255) cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAAc < 4:188497/50‑1 (MQ=255) aaaCTAGCAAAGATTTTGTAGCCTGTACTGGCCTGGTAGCCGGTGGAACCCCACCTAAAGGGGGTGAAACGGTGGGAGCGGGCTTGCCCGCGCTAA‑GGCCGCTAAATGCTTTACCGAATAACGCACGACGCAGTcc > 8:278393/1‑135 (MQ=255) gggCCTCTTCGGGGGTCACCCCGCCCATACGGGGGGATACCCTGGGGGCGGGGGCTGGCCCGCTACGA‑GGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTg < 4:273760/134‑1 (MQ=255) | GCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGA‑GGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTG > ACB143/5251052‑5251309 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |