Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 322,525 | 0 | A | C | 58.8% | 13.8 / 15.6 | 17 | V145G (GTG→GGG) | betA | choline dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (5/2); new base C (0/10); total (5/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.39e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGA > NZ_CP009273/322436‑322603 | cGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCa < 2:293845/90‑1 (MQ=255) tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc < 2:229794/90‑1 (MQ=255) gcgcGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGcc > 1:232420/1‑90 (MQ=255) cccGCCCGCTTAAATCACCGCTTCAAACAGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTc < 1:234948/88‑1 (MQ=255) cccGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGttt > 1:328983/1‑90 (MQ=255) cccGCTTCAATCATCGCTTAAAACAGCGGATTGACGCGGGGTTGGGAGGTAGTGACCCTCACCGGGCCATCACCGCCGTGATAGTCGttt < 1:216006/90‑1 (MQ=255) cccGCTTAAATCACCGTTTCAAACGGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGttt < 2:82827/90‑1 (MQ=255) ttaaaTCATCGCTTCAAACAGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCGGGGCCATCACCGCCGTGATAGTCGGTTTCAcc < 1:320954/87‑1 (MQ=255) taaaTCATCCCTTAAAACAGCGGTTTGCCGCGGGGTTTGGGGGTAGTGCCGCCCCCGGGGCCATCCCCGCCGTGATAGTCGTTTTCAccc < 2:328983/88‑1 (MQ=255) taaaaCAGCGGATTGCCGCGGGGTTGGGGGGTAGTGACCCCCCCGGGGCCATCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtc < 2:196331/88‑1 (MQ=255) tCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtc > 1:234540/1‑90 (MQ=255) aaaCAGCGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtctc < 1:33483/90‑1 (MQ=255) gATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc < 1:111203/90‑1 (MQ=255) gggAGGTAGTGCCGCCCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggcaggc < 1:291834/89‑1 (MQ=255) ggTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTc < 1:177260/90‑1 (MQ=255) ggTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTc < 1:16757/90‑1 (MQ=255) tAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGa > 2:162484/1‑90 (MQ=255) tAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGa > 1:150754/1‑90 (MQ=255) | CGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGA > NZ_CP009273/322436‑322603 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |