Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,978,077 | 0 | G | C | 61.1% | 7.1 / 19.3 | 18 | P292R (CCG→CGG) | yqeF | acetyl‑CoA C‑acetyltransferase |
Reads supporting (aligned to +/- strand): ref base G (2/5); new base C (11/0); total (13/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.45e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTG > NZ_CP009273/2977995‑2978166 | aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCgcgccaa > 2:369739/1‑87 (MQ=255) cGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAAATcc > 1:56651/1‑90 (MQ=255) aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt < 1:369991/90‑1 (MQ=255) aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt < 1:645632/90‑1 (MQ=255) aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAACCCGCGCCAATCCCCTCCATGCCGGATCt > 2:597534/1‑90 (MQ=255) gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGGGCATAAACCCGCGCCAATCCCCCCCAAACCCGAGCTACACCAAc > 1:177087/1‑90 (MQ=255) ccAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCAATCCCCTCAATGCCCGGTCTACACCCACCCCGGCAAAAGcgcg > 1:326835/1‑90 (MQ=255) acacGCTCCAGGCAACGGCGGGTCGCCTACACCGGGGCAAATCCCATCAAAGCCGGGGCTACACCAAAGCCGGCCAAAGCGCGAATGCgg > 2:49220/1‑90 (MQ=255) cTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCAATGCCGGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAg > 2:69580/1‑90 (MQ=255) cTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTCCCCCCAATGCCGGGTCTACACCAACGCCGGCAAAAGCGCGAAAGCGGGCCAg > 1:99103/1‑90 (MQ=255) cAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCAc < 2:177087/90‑1 (MQ=255) cAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCAc < 1:597534/90‑1 (MQ=255) aaCGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTa < 2:161912/90‑1 (MQ=255) cggcggGTCGCATACACCCGCGCAATTCCCCTCAATGCCGGATCTACACCCACGCCGGCAAAAGCGCGAATGCGGGCCAACACgggggaa > 1:4866/1‑86 (MQ=255) cggcggGTCGCATACACCCGCGCAATTCCCATCAATGCCGGGTCTACACCAACGCTGGCAAAAGCGCGGATGCGGGGCCGCCCGGGGaaa > 2:404605/1‑90 (MQ=255) cggcggGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCCGCACGGGGaaa > 2:45018/1‑90 (MQ=255) cggcggGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCCGCACGGGGaaa > 2:314280/1‑90 (MQ=255) ccGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt > 2:357054/1‑90 (MQ=255) ggCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTg > 2:124039/1‑90 (MQ=255) ggCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTg > 2:60802/1‑90 (MQ=255) | AAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTG > NZ_CP009273/2977995‑2978166 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |