Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,822,630 | 0 | C | G | 25.0% | 79.2 / 29.8 | 44 | F207L (TTC→TTG) | bamA‑II | outer membrane protein assembly factor |
Reads supporting (aligned to +/- strand): ref base C (10/23); new base G (0/11); total (10/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.57e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGCAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGT > NC_002947/3822483‑3822646 | gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCa < 1:8701003/149‑1 (MQ=32) gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCa < 1:7510249/149‑1 (MQ=32) gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCa < 1:7510219/149‑1 (MQ=32) gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCa > 2:8162432/1‑149 (MQ=32) gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCa < 1:6216919/149‑1 (MQ=32) cAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACACCGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGa < 1:711476/149‑1 (MQ=17) cAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGa < 1:6750401/149‑1 (MQ=32) cAACCGTGTCGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGa < 2:3271639/149‑1 (MQ=17) aaCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGaa > 1:1543252/1‑149 (MQ=32) aCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac > 1:436006/1‑149 (MQ=32) aCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac > 1:435618/1‑149 (MQ=32) aCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac > 1:1122545/1‑149 (MQ=32) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacg < 2:4345532/149‑1 (MQ=17) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacg < 1:2401659/149‑1 (MQ=32) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacg < 2:4309511/149‑1 (MQ=32) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacg < 1:80456/149‑1 (MQ=32) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacg < 1:1906696/149‑1 (MQ=32) cGTGTTGCGCAGAAGATCAAGATCAACGAAGGCACCGTCGCAGCAATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCATGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGATTTCGATATTCAATAaaga > 2:7683730/1‑149 (MQ=11) cGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacga < 2:5125139/149‑1 (MQ=17) cGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacga > 2:5246715/1‑149 (MQ=32) cGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacga < 2:8522820/149‑1 (MQ=32) cGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacga < 2:8589945/149‑1 (MQ=32) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgac < 2:2441854/149‑1 (MQ=17) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgac < 2:2819192/149‑1 (MQ=17) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgac < 2:2819317/149‑1 (MQ=17) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgac < 2:3329044/149‑1 (MQ=32) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgac < 2:3329109/149‑1 (MQ=32) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca < 1:7644187/149‑2 (MQ=17) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca < 1:7644997/149‑2 (MQ=17) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca < 2:2280164/149‑2 (MQ=17) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca < 2:4195988/149‑2 (MQ=17) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgacg > 1:7806353/1‑149 (MQ=32) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacgacg < 1:6328818/149‑1 (MQ=32) gtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgacag < 1:3672243/149‑3 (MQ=17) gtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAAcgacga < 2:7621326/149‑1 (MQ=32) gtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAAcgacga < 1:4579741/149‑1 (MQ=32) gcgcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACaa < 2:1051147/149‑1 (MQ=32) gcgcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACaa > 2:542589/1‑149 (MQ=32) gcgcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACaa > 2:5458817/1‑149 (MQ=32) cgcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgacagccag < 1:4122662/149‑7 (MQ=9) cgcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAg < 2:3926041/149‑1 (MQ=31) gcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGt < 2:3042667/149‑1 (MQ=11) gcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGt < 2:2629001/149‑1 (MQ=11) gcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGt < 1:3246349/149‑1 (MQ=11) | GCGCAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGT > NC_002947/3822483‑3822646 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |