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.6% | 23.3 / 35.6 | 28 | G165G (GGT→GGG) | ycfH | metal‑dependent hydrolase |
Reads supporting (aligned to +/- strand): ref base T (3/10); new base G (15/0); total (18/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.18e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152646‑1152816 | gCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATACt > 2:271441/1‑90 (MQ=255) cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATACTg > 1:1005354/1‑90 (MQ=255) tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTGAAATaattg > 2:1107005/1‑86 (MQ=255) cGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGTGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCg > 2:1042202/1‑90 (MQ=255) cGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATTCTGGGTTTCGGGTTTTAAATTTCCtttttcc > 2:419727/1‑87 (MQ=255) ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCgg < 1:1107005/90‑1 (MQ=255) gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTTTCGGGTTTTATATTTCCTTTTTCTGGATtg > 2:486361/1‑90 (MQ=255) gTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGa < 1:1168635/90‑1 (MQ=255) ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATAATGGATTTCTGGTTTTAAAAATCCTTTTTCCGCCTtgtggac > 1:898536/1‑87 (MQ=255) gCGTACTACACTGTTTTACCGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATGTTACATCTCCTTTTCCGGCATTGTGACCt < 1:949993/90‑1 (MQ=255) gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCt < 1:1260615/90‑1 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATTGATCTCGGGTTTTTAATTACCTTTTTTGGCGGTGTGAGCTTCCGTAATgg > 1:656880/1‑89 (MQ=255) cTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTCGGGTTTTTTATCTCCTTTTTCGGCATTTTGGGCCTCCGTAATGc > 1:1376688/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTTTATCTCCTTTTTCGGCATTGTGACCTTCCGTAATGCgg > 2:589998/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGATCTCGGGTTTTTTATCTCCTTTTTCCGCATTGTGACCTTCCCTAATTCgg > 2:1176142/1‑90 (MQ=255) gTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGATCTCGGGTTTTACATCTCCTTTTTCCGCATTGTGGCCTTCCCTAATTCgg > 2:818207/1‑90 (MQ=255) cAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAAc < 2:838993/90‑1 (MQ=255) cAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAAc < 1:1322986/90‑1 (MQ=255) gagGACAGAGAAACGGCGGGGAAAATACTTGGTCTCCGGTTTTAAATCTCCTTTTTCCGCATTGTGACCTTCCCTAATGCGGAGGAAATg > 2:883356/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATACTGGGTCTCCGGTTTTAAATCTCCCTTTTCCGCATTGTGGCCTTCCCTAATACGGGGCAACTgcg > 2:1362438/1‑90 (MQ=255) ggACAGAGAAACGGCGGGGAAAATAATTGGTTTCGGGTTTTTTATCTTCTTTTTTGGCCTTGTGGCCTTTCCTAATTGGGGGCAAATgcg > 1:1406760/1‑90 (MQ=255) gagaAACGGCGGGGAAATTACTGGGTCTCGGGTTTTAAATCTCCTTTTTCGGCCTTTTGACCTTCCCTAATGCGGGGCAACTGCGCGGtg > 2:499149/1‑90 (MQ=255) gagaAACGGCGGGGAAAATAATTGATCTTGGGTTTTTAATCTCCCTTTTCCGCCTTTTGACCTTTCCTAATGCGGGGGAACTGCGCGGtg > 1:1347942/1‑90 (MQ=255) agaAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgc < 2:1406760/90‑1 (MQ=255) agaAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgc < 1:959213/90‑1 (MQ=255) aCGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcg < 1:1102441/90‑1 (MQ=255) aCGGCGGGGAAATTACTGGATCTCGGATTTTAAATCTCCTTTTTCCGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcg > 2:392500/1‑90 (MQ=255) ggTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTAt < 2:19601/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 2:1347942/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:420105/90‑1 (MQ=255) tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGt < 1:191150/90‑1 (MQ=255) | GCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGT > NZ_CP009273/1152646‑1152816 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |