Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,906,007 | 0 | T | G | 52.4% | 16.5 / 27.1 | 23 | H276P (CAC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (2/8); new base G (11/0); total (15/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT > NZ_CP009273/1905919‑1906092 | gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 1:282771/90‑1 (MQ=255) gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:511979/90‑1 (MQ=255) gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:431174/90‑1 (MQ=255) gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:159372/90‑1 (MQ=255) aCGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATg < 1:471262/90‑1 (MQ=255) cgcgCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa < 1:672029/90‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGaaa > 1:232555/1‑89 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaaa > 1:649398/1‑89 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc > 1:828671/1‑90 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGGTGAGCAACTAACTATGCGGttttc > 2:1202503/1‑88 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:295925/90‑1 (MQ=255) tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:283992/90‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 1:978153/54‑1 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg > 2:978153/1‑54 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGttt > 2:1173454/1‑89 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTCCACTTACCGTTaa > 1:343535/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCCCCGAGCGGTTTTGACTTTCCCTTAATGCag > 2:932981/1‑90 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 2:803143/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc > 1:457383/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAAAAAATCCCCTAGCGGTTTTGACTTACCGTTTATGGAGGGAGcc > 2:905237/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTTACTTTCCGTTAATGCAGAGAGcc > 2:961748/1‑90 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTTAGCGATTTTGACTTACCGTTTATGCAGAGAGcc > 2:278036/1‑90 (MQ=255) tACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGAGGGGCATGAAAAACTCCCTGAGCGATTTTGACTTACCGGTAAAGCAGAGAGCca > 1:597309/1‑90 (MQ=255) cGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatc < 2:1020311/90‑1 (MQ=255) gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatca < 1:891641/90‑1 (MQ=255) gTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGGTGcttctt > 1:674193/1‑90 (MQ=255) gTTTATCCAGCGGCGGGGGGGGAATCAACAAATCACTGGGCGATTTTGGCTTACCGTTAATGGAGGGGGGCAATATGCTGGTTGgttctt > 1:1167533/1‑90 (MQ=255) gcggcgGGAGGGGCATGGAAAACACACCGAGCTATTTTGACTTCACCTTAAACCAGAGAGCCATCATGTTGGGTGGGTATTGGGGTTctt > 2:995087/1‑88 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:995087/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 < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |