Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,152,727 | 0 | T | G | 54.5% | 21.1 / 28.7 | 22 | G165G (GGT→GGG) | ycfH | metal‑dependent hydrolase |
Reads supporting (aligned to +/- strand): ref base T (0/10); new base G (12/0); total (12/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152640‑1152816 | cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTaa < 1:843385/90‑1 (MQ=255) tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTgg > 1:128579/1‑90 (MQ=255) tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATgg > 1:1286098/1‑90 (MQ=255) tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGAt < 2:902611/90‑1 (MQ=255) cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAAAGGACAGAGAAACGGCGGGGGAAATAATGGATTTCGGGttt > 2:798188/1‑90 (MQ=255) aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGTTTCGGATTTTaaa > 1:190366/1‑88 (MQ=255) aaaaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACAtct < 1:1274739/90‑1 (MQ=255) ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCgg < 1:450489/90‑1 (MQ=255) gTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGa < 2:417278/90‑1 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCt < 1:1215700/90‑1 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCt < 1:667889/90‑1 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATATTGGGTCTCGGGTTTTTAATTTCCTTTTTCGGCATTTTGGCCt > 1:417278/1‑90 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATCTCGGATTTTAAATCTCCTTTTTCTGCATTTTGAACt > 1:621283/1‑90 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTCGGGTTTTAAATCTCCTTTTTCGGCATTTTGACCTTCCGTAATGc > 1:519427/1‑90 (MQ=255) aggTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGc < 1:732020/88‑1 (MQ=255) ttACAGAGGACAGAGAAACGGCGGGGAAAATATTGGGTCTCGGGTTTTTTATCTCCTTTTTCCGCATTGTGGCCTTTCGTAATggggggg > 2:343005/1‑87 (MQ=255) ggACAGAGAAACGGCGGGGAAAATATTTGGTTTCGGGTTTTTAATCTCCTTTTTCGGCATTGTGACCTTCCGTTATGCGGGGGAACTgcg > 2:25298/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATGGGTTTCGGGTTTTTCATCTTCTTTTTCTGCCTTGTGGCCTTTTGTTATGCGGGGGAACTgcg > 2:704673/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATGGGTCTCGGGTTTTTCATCTCCTTTTTCCCCGTTGTGGGGTTCCGTAATGCGGGGGAAATgct > 2:732020/1‑89 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATGGATCTCGGGGTTTAAATCTCCTTTTTCCGCAATTTGACCTTCCCGAAAGCGGGGCAACTgcg > 1:865895/1‑90 (MQ=255) cggcggGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTTACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcgc < 1:47258/90‑1 (MQ=255) gcggGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGt < 2:1213374/90‑1 (MQ=255) ggggAAAATACTGGGTCTCGGGGTTTTCATCTTCTTTTTCGGGATTGTGACCGTCCGTAATGCGGGGCAACTGCGCGATTCTTcgcggtt > 2:349233/1‑87 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCTGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:704673/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 2:621283/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:991919/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:919043/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:1285981/90‑1 (MQ=255) | CGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152640‑1152816 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |