Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,152,727 | 0 | T | G | 53.3% | 14.8 / 20.0 | 15 | G165G (GGT→GGG) | ycfH | metal‑dependent hydrolase |
Reads supporting (aligned to +/- strand): ref base T (0/7); new base G (8/0); total (8/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152638‑1152816 | gCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGt < 2:667259/90‑1 (MQ=255) ttCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCg < 2:701483/90‑1 (MQ=255) cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATATTGGATCTCTGGttt > 1:869345/1‑90 (MQ=255) aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTT‑ACa < 2:568617/90‑1 (MQ=255) gTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTT‑ACATCTCCttt < 2:253663/90‑1 (MQ=255) aCGGATCGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATTATGGGAT‑TTGGGTTTT‑TCATCTCCtttttt > 1:611615/1‑88 (MQ=255) gTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATAATGGATTTCGGGTTTT‑TTATTTCCTTTTTCCGCCTTTTGGGCTTTCCTAATGTgg > 1:929767/1‑90 (MQ=255) cAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAAc < 2:929767/90‑1 (MQ=255) gagGACAGAGAAACGGC‑GGGGAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTTCGGCATTGTGACCTTCCGTAATGCGGAGGAACTg > 2:708154/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGGAAATTA‑TGGGTTTCTGGTTTT‑ACATCTCCCTTTTCGGGATTTTGACCTTCCCTAATGCGGGGGAAcatcg > 2:552915/1‑86 (MQ=255) ggACAGAGAAACGGC‑GGGGAAATTAATTGATCTCGGATTTT‑ACATCTCCTTTTTCGGCATTGTGGCCTTCCCGAATGTGGAGGAACTtcg > 2:1065607/1‑90 (MQ=255) ggACAGAGAAACGGC‑GGGGAAAATAATGGATCTCTGGTTTT‑TCATCTCCTTTTCCCGCATTTTGACCTTCCGTAATGCGGAGCAAATgcg > 2:398160/1‑90 (MQ=255) ggACAGAGAAACGGC‑GGGGAAAATAAATGGTCTCGGGTTTTAAAAATTCCTTTTCCGG‑ATTGGGGCCTTCCCGAATGCGGAGAAACTtcg > 1:148475/1‑90 (MQ=255) ggc‑ggGTAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGcgcg < 2:93280/90‑1 (MQ=255) gggTAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTa < 1:525885/90‑1 (MQ=255) ggTAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTAt < 2:620921/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 2:844681/90‑1 (MQ=255) | GCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTT‑ACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152638‑1152816 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |