| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,336,492 | 0 | T | G | 42.5% | 59.3 / 14.2 | 40 | T218P (ACC→CCC) | fumB | class I fumarate hydratase |
| Reads supporting (aligned to +/- strand): ref base T (9/14); new base G (17/0); total (26/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.76e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGC > NZ_CP009273/4336404‑4336581 | aaCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTa < 2:96613/90‑1 (MQ=255) tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGGGGCAGGACCGGGGGGACGCa > 2:747380/1‑90 (MQ=255) tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTAACGGGGGGACGcc > 1:156756/1‑89 (MQ=255) aGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCAGTACCGAGGGGACGCATTTTTTcgacg > 1:1150740/1‑90 (MQ=255) tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGa > 1:324219/1‑90 (MQ=255) aCCGCCAATCACAAACGCGATATGGTACGGCGGGCGGGCTGCAGGACCCAGGGGACGCATTTTTTCGACGGGGGAGGTTTtttgtttggc > 2:186970/1‑88 (MQ=255) aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc < 1:1381726/90‑1 (MQ=255) gCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGGTTTTCAGTTTGCCggg > 1:252798/1‑90 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGGCGCCGGACCCAGGGGACGCCTTTTTTCGACGGGGGAGGTTTTTTGTTTTCCggggg > 2:1355850/1‑88 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAg < 2:324219/90‑1 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAg < 2:252101/90‑1 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGGCggggg > 1:684965/1‑88 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGCCggggg > 1:29489/1‑88 (MQ=255) acaAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagc < 2:626219/90‑1 (MQ=255) acagcgcgATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagca < 1:7441/86‑1 (MQ=255) gATATGGTACGGCGGGGAGGCGGCAGGACCGAGGGGACGCCTTTTCTCGACGAGGGAGGTTTTTAGGTTGGCGGGGGGCAGCAGGGCttt > 1:726019/1‑90 (MQ=255) gATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCGGGGGGCAGCAGGGCgtt > 1:1067891/1‑90 (MQ=255) atGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGCGTTGGt > 1:1149788/1‑90 (MQ=255) tACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTTCCGGGGGGCAGCAGGGGTTTGGtctgt > 2:7441/1‑86 (MQ=255) aCGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt < 2:394672/90‑1 (MQ=255) cggcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTTGTTTTCCGGGGGGCAGCAGGGGTTTGGTTTTCTg > 1:1521052/1‑90 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGt < 1:589566/90‑1 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGt < 2:1067891/90‑1 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGt < 2:1521052/90‑1 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTTCTggg > 1:460497/1‑89 (MQ=255) gcggGCAGGCTGCAGGACCGAGGGGACGCATTTTTTTGACGGGGGAGGTTTTTTGTTTGCCGGGGGGCAGCAGGGGTTTGGGTTTCTggg > 2:1155727/1‑89 (MQ=255) gcggGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTCCTggg > 2:1492801/1‑89 (MQ=255) cAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAgaga < 1:1492801/90‑1 (MQ=255) cAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAgaga < 2:389860/90‑1 (MQ=255) cAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc > 1:921640/1‑49 (MQ=255) cAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc < 2:921640/49‑1 (MQ=255) gTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttg < 2:1150740/90‑1 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 2:58580/1‑90 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 2:1506210/1‑90 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 1:410469/1‑90 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 2:1136422/1‑90 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 1:74203/1‑90 (MQ=255) gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag > 1:664112/1‑90 (MQ=255) gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag > 1:1285437/1‑90 (MQ=255) gggTACGCATGTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAACAGGGCTTTGGTTTCCTGGTAGAGGTACGTTTTGTTGGCag > 1:664113/1‑90 (MQ=37) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGt < 2:679030/66‑1 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGt > 1:679030/1‑66 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc > 1:1283583/1‑90 (MQ=255) | AACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGC > NZ_CP009273/4336404‑4336581 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |