Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,781,385 | 0 | T | G | 53.3% | 25.5 / 40.5 | 30 | G396G (GGT→GGG) | envC | murein hydrolase activator EnvC |
Reads supporting (aligned to +/- strand): ref base T (5/9); new base G (16/0); total (21/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTTC > NZ_CP009273/3781297‑3781473 | gCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTc < 1:213896/90‑1 (MQ=255) ataATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGggcggc > 1:454878/1‑90 (MQ=255) gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAct > 2:733659/1‑90 (MQ=255) cACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTAtt > 2:831040/1‑90 (MQ=255) gTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCCTCTATTTCGaa > 2:728472/1‑90 (MQ=255) gAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCg < 2:933918/69‑1 (MQ=255) gAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCg > 1:933918/1‑69 (MQ=255) gAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCCGGGGCGGCCCTCACCCCATTTCGTAAt > 1:415862/1‑90 (MQ=255) gTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTcgc < 1:808241/90‑1 (MQ=255) ttGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCAc < 2:119430/72‑1 (MQ=255) ttGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAc > 1:119430/1‑72 (MQ=255) ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCCTCTATTTCGAAATTcgccgc > 1:126905/1‑90 (MQ=255) gCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCCGGGTCGGCCTTCTCTCTATTTCGAAAATCGCCGCCAggg > 2:419847/1‑90 (MQ=255) aGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTc < 2:93993/90‑1 (MQ=255) ttCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTTGCCGCCCGGGGCAgg > 1:218491/1‑90 (MQ=255) gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 2:454878/90‑1 (MQ=255) gcgGGCCACCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 1:898713/90‑1 (MQ=255) ttGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGt < 2:707827/90‑1 (MQ=255) cTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCActc > 2:604593/1‑39 (MQ=255) cTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCActc < 1:604593/39‑1 (MQ=255) ggtggGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGTCAGGCGGGCAATCCACAGCCGGGGTTggg > 2:602910/1‑90 (MQ=255) gcagcaGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGGGGGTGGGGAGAt > 2:785603/1‑90 (MQ=255) cagcagTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCCGGGGCAGGCGGTCAATCCCCAGCCGTGGGTGGGGAGATa > 1:278285/1‑90 (MQ=255) agTGGCGGTTACGGGTGGCCCTCACCCTATTTCGAAATTCGCCCCCAGGGGCAGGGGGGCAATCCCCAGCCGTTGGTGGGGAGGTaaagt > 1:513669/1‑87 (MQ=255) agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCCCAGCCGTGGGTGGGGAGATAAGtt > 1:226765/1‑90 (MQ=255) agTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGGGGGTGGGAAGATAAGtt > 2:191158/1‑90 (MQ=255) agTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGGGGTTGGGGAGATAAGtt > 1:150971/1‑90 (MQ=255) ggCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGGGGGTGGGGAGATAAGttttt > 2:678595/1‑89 (MQ=255) ggCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCCGGGGCAGGCGGGCAATCCCCAGCCGGGGGTGGGGAGATAAGttttt > 1:380923/1‑89 (MQ=255) ggTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGttt > 2:378346/1‑90 (MQ=255) gTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTtc > 2:451951/1‑90 (MQ=255) | GCTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTTC > NZ_CP009273/3781297‑3781473 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |