Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 322,525 | 0 | A | C | 53.8% | 12.4 / 13.7 | 13 | V145G (GTG→GGG) | betA | choline dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (4/2); new base C (0/7); total (4/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.10e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTA > NZ_CP009273/322436‑322607 | cGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCa > 1:264755/1‑90 (MQ=255) gTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCggg > 1:1764/1‑90 (MQ=255) aCGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCg > 1:112495/1‑90 (MQ=255) tCAATCACCGCTTCAAACGGGGGTTGGGCGCGGGGTTGGGGGGTAGTGACCCTCCCGGGGCCATCCCCGCCGTGATAGTCGTTTTCAccc < 2:502498/90‑1 (MQ=255) ccGCTTCAAACAGCGGATTGCCGCGGGGTTGGGAGGTAGTGACGCCCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATcgc < 2:206353/89‑1 (MQ=255) cAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGc < 1:375139/90‑1 (MQ=255) cGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc > 1:313605/1‑37 (MQ=255) cGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc < 2:313605/37‑1 (MQ=255) gtttGACGCGGGGTTGGGGGGTAGTGCCGCCCCCGGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 2:478376/88‑1 (MQ=255) gATTGCCGCGGGGTTTGGAGGTGGTGCCGCCCCCCGGCCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 2:502327/90‑1 (MQ=255) gATTGAGGCGGGGTTGGGAGGTAGTGCCGCTCCCCGGGCCATCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 2:1764/90‑1 (MQ=255) gATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 1:348012/90‑1 (MQ=255) aGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAg > 1:260111/1‑90 (MQ=255) gACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTa < 1:93783/90‑1 (MQ=255) | CGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTA > NZ_CP009273/322436‑322607 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |