Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 322,502 | 0 | C | G | 54.5% | 2.7 / 13.5 | 11 | G153R (GGC→CGC) | betA | choline dehydrogenase |
Reads supporting (aligned to +/- strand): ref base C (4/1); new base G (0/6); total (4/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGC > NZ_CP009273/322413‑322582 | ttCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc > 1:406185/1‑90 (MQ=255) ttCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc < 2:263882/90‑1 (MQ=255) ttCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc < 2:44672/90‑1 (MQ=255) ttCCTGCTAATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc > 2:166669/1‑90 (MQ=255) gctgATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGtt > 1:348885/1‑90 (MQ=255) cGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCa > 1:106254/1‑90 (MQ=255) tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc < 1:11415/90‑1 (MQ=255) tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc > 1:241379/1‑90 (MQ=255) cTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATa > 2:366320/1‑90 (MQ=255) cccGCTTCAACCCCCGTTCCAAACCGCGGTTTGCGGCGGGGTTTGGGGGTAGTGCCCCCCCCCGGCCCCTCCCCGCCGTGATAGTCGttt < 2:209622/90‑1 (MQ=255) cctccAAACAGCGGTTTGACGCGGGGTTTGGAGGTAGTGCCGCCCCCCGGGCCATCTCCGCCGTGATAGTCGTTTTCACCCATATCGCGa < 2:348885/86‑1 (MQ=255) aaaCAGCGGTTTGCCGCGGGGTTTGGGGGTAGTGACGCCCCCCGGGCCACCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtctc < 2:241379/90‑1 (MQ=255) cGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCtt < 1:533865/90‑1 (MQ=255) gATTGGCGCGGGGTTGGGGGGTGGTGCCGCCCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 2:579438/90‑1 (MQ=255) gATTGACGCGGGGTTGGGGGGTGGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 2:39596/90‑1 (MQ=255) | TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGC > NZ_CP009273/322413‑322582 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |