Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,906,007 | 0 | T | G | 51.6% | 35.8 / 45.2 | 31 | H276P (CAC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (2/13); new base G (16/0); total (18/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905919‑1906092 | gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:2030976/90‑1 (MQ=255) gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:324271/90‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaac > 1:1985981/1‑90 (MQ=255) tttttAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaacaac < 2:1058148/90‑1 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAAGAACAAATCACTGAGCgt > 1:1875179/1‑89 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 1:1448281/61‑1 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAAACCCCTCAGCGATTTct > 2:1345748/1‑89 (MQ=255) tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAAAAACTCACTGAGCGAttttttc > 2:1180113/1‑86 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGACCAACTCACTGAGCGATTTCGACTTa < 2:1129695/90‑1 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 1:778058/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTAcc < 1:734219/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAAAAAATCACTCAGCGATTTTGACTTAcc > 2:2181282/1‑90 (MQ=255) cAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGtt < 2:479260/90‑1 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAAAAAATCACTGAGCGATTTTGACTTACCGTTaa > 2:39535/1‑90 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTGAGCGGTTTTGACTTACCGTTAATGCa > 1:1931308/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 1:1180113/90‑1 (MQ=255) cagcccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 1:631235/87‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 1:807865/90‑1 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGc < 2:727644/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATTAAGAACTCACTGAGCGGTTTTGGCTTACCGTTAATGCagagaggc > 1:1939647/1‑88 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACCCCCCTGGCGGTTTTGGCTTACCCTTAAAGCGGAGGGcc > 2:1075854/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTTAGCGATTTCGACTTACCGTTAATGCAGAGAGCCAACATGc > 1:331841/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCACAGAGCGATTTCGGCTTACCGGTAATGGAGAGAGCCATTATGc > 1:1228909/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACCGAGCGGTTTTGACTTTCCCTTTATGCAGAGAGCCCTCATGc > 2:760151/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCCCTGAGCGATTTTGACTTTCCCTTAATGCAGAGAGCCATCATGc > 2:1877174/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGGTTTTGACTTTCCGGTTATGCAGAGAGCcatcactg > 2:631235/1‑87 (MQ=255) tCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGtt > 1:876310/1‑90 (MQ=255) tCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAAACCCTGAGCGGTTTTGACTTTCCCTTAATGCAGAGAGCCATCATGCTgggt > 2:1977721/1‑88 (MQ=255) gTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt > 1:1587810/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:212967/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:616430/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:1197893/90‑1 (MQ=255) gcggcgGGGGGGTCATGAACAAATCACTGAGCGATTTTGACTTTCCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAt > 2:1378815/1‑90 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 2:1931308/90‑1 (MQ=255) | GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905919‑1906092 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |