Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,336,492 | 0 | T | G | 51.3% | 37.0 / 59.3 | 37 | T218P (ACC→CCC) | fumB | class I fumarate hydratase |
Reads supporting (aligned to +/- strand): ref base T (4/14); new base G (19/0); total (23/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.01e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGC > NZ_CP009273/4336405‑4336581 | aCTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTAc < 1:965497/90‑1 (MQ=255) ttGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGc < 1:1375938/90‑1 (MQ=255) tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAGGACCCGGGGGACGCa > 2:1022977/1‑90 (MQ=255) tGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGCCGCAGGACCCAGGGGACGCa > 2:1171481/1‑90 (MQ=255) tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGGACCGGGGGGACGCATTTTTTCGACGAggg > 2:168010/1‑89 (MQ=255) tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGGGCAGGACCCAGGGGACGCCATTTTTCGACGAggg > 1:62329/1‑89 (MQ=255) tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGa > 2:194944/1‑90 (MQ=255) tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGACGAggg > 2:238632/1‑89 (MQ=255) tCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGGGCAGGACCGAGGGGGCGCATTTTTTTGACGGGGGAGgtt > 1:1402847/1‑90 (MQ=255) cGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGGTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGttttt > 2:1321777/1‑90 (MQ=255) gACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGt < 2:418224/90‑1 (MQ=255) aCGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtt < 1:1265371/90‑1 (MQ=255) aCCGCCAATCACAAACGCGATATGGTACGGCGGGGAGGGCGCCGGACCCGGGGGGCGCATTTTTTCGACGAGGGAGGTTTTTAGTTTTcc > 1:48094/1‑90 (MQ=255) aCCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGcc < 2:177171/90‑1 (MQ=255) cAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCACTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTTCCggggg > 2:903685/1‑88 (MQ=255) aaaCGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagcag < 2:288891/90‑1 (MQ=255) atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCa < 2:962781/40‑1 (MQ=255) atatGGTACGGCGGGCAGGCTGCAGTACCGAGGGGACGCa > 1:962781/1‑40 (MQ=255) atGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt < 2:1262898/90‑1 (MQ=255) aCGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt < 1:549195/90‑1 (MQ=255) aCGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCt < 1:986215/90‑1 (MQ=255) gcggGCAGGCTGCAGTACCGGGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTGGTTTTCTggg > 2:733947/1‑89 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTCCTggg > 2:559711/1‑89 (MQ=38) gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGGAGGTTTTCAGTTTGCCGGGGGGCAGCAGGGGTTTGGTTTTCTGGt > 2:383831/1‑90 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTggg > 2:194467/1‑89 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTCCTGGt > 2:64961/1‑90 (MQ=255) gcggGCAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTTAGTTTTCCGGGGGGCAGCAGGGCTTTTGGTTTCTGGt > 2:464603/1‑90 (MQ=255) gcggGCAGGCTGCAGTACCCAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGGTTTTGGTTTCTGGt > 1:545347/1‑90 (MQ=255) gcggGCAGGCTGCAGTACCCAGGGGACGCATTTTTTCGACGAGGAAGGTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGGTTTCTggg > 2:669495/1‑89 (MQ=255) cggGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTa < 1:781194/90‑1 (MQ=255) gggCAGGCTGCAGGACCGAGGGGACGCATTTTTTCGGCGGGGGAGGTTTTTTGTTTGTCGGGGGGCAGCAGGGGTTTTGTTTTCTGGGag > 1:382308/1‑90 (MQ=255) cTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACg < 1:595966/90‑1 (MQ=255) gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGtt < 1:464603/90‑1 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag > 1:143032/1‑73 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTag < 2:143032/73‑1 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 1:876768/1‑90 (MQ=255) gAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGc > 2:569087/1‑90 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc > 2:400835/1‑90 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc > 2:844567/1‑90 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc > 2:1284378/1‑90 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc > 2:302228/1‑90 (MQ=255) tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgc > 1:484272/1‑90 (MQ=255) | ACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGC > NZ_CP009273/4336405‑4336581 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |