Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,781,385 | 0 | T | G | 64.3% | 14.1 / 28.3 | 28 | G396G (GGT→GGG) | envC | murein hydrolase activator EnvC |
Reads supporting (aligned to +/- strand): ref base T (3/7); new base G (18/0); total (21/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTTC > NZ_CP009273/3781300‑3781473 | ataATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGggcggc > 1:199011/1‑90 (MQ=255) aaTCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCt < 2:414549/90‑1 (MQ=255) aaTCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGCCt > 2:20015/1‑90 (MQ=255) gcACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCActctct > 1:347753/1‑88 (MQ=255) gcACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGCCCTCCCTCTAt > 2:291304/1‑90 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTc < 1:158470/90‑1 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTTTTTTGAAATTc > 2:254964/1‑90 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCCTCCCTCTATTTTGGAAttt > 1:341864/1‑89 (MQ=255) gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGGGGCCTTCACTCTATTTTGAAATTc > 2:158470/1‑90 (MQ=255) ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTcgccgc > 2:269005/1‑90 (MQ=255) ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTcgccgc > 1:138412/1‑90 (MQ=255) ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGCCTTCACTCTATTTCTAAATTcgccgc > 1:434971/1‑90 (MQ=255) gCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGGCTTCACTCTATTTCGAAATTCGCCGCCCggg > 1:262377/1‑90 (MQ=255) gCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAggg > 1:477075/1‑90 (MQ=255) gCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCGGGGGCGGGCCTCCCCCTATTTTGGGATTCGGCGGCgggg > 2:216294/1‑90 (MQ=255) ttCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGGCAgg > 1:81918/1‑90 (MQ=255) gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc < 2:83141/90‑1 (MQ=255) aGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCcac < 2:363838/90‑1 (MQ=255) ttGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGGAATTTGCCGCCAGGGGCGGGGGGGCAAACCCCAGCCGt > 2:231708/1‑90 (MQ=255) cTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCCTCTATTTTGAAATTCGCCGCCCGGGGCAGGGGGGCAATCCACAGCCGTGGTTg > 2:454456/1‑90 (MQ=255) ggtggGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCCATTTCGGAATTCTCCCCCCGGGGCACGGGGGCCCTCCCCACCCGCGGGTggg > 2:278737/1‑90 (MQ=255) agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGtt > 1:22029/1‑90 (MQ=255) gTGGCGGTCAGGGGCGGCCCTCACTCTATTTTGAAATTCGCCCCCCGGGGCAGGCGGGCCATCCCCAGCCGGGGGTGGGGAGAGAAGttt > 2:233412/1‑90 (MQ=255) ggTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTg < 2:138782/72‑1 (MQ=255) ggTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTg > 1:138782/1‑72 (MQ=255) tCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGTGGTTGGGGAGATAAGTTTTGTTTcc > 1:106219/1‑90 (MQ=255) gggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAtt > 2:400960/1‑90 (MQ=255) ggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAttt < 1:454456/90‑1 (MQ=255) gTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTtc < 1:269005/90‑1 (MQ=255) | ATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTTC > NZ_CP009273/3781300‑3781473 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |