Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,819,603 | 0 | C | G | 40.9% | 16.8 / 25.2 | 22 | F207L (TTC→TTG) | bamA‑II | outer membrane protein assembly factor |
Reads supporting (aligned to +/- strand): ref base C (6/7); new base G (1/8); total (7/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGCGCAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTA > NC_002947/3819454‑3819620 | ccGCGCAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGttcttc < 2:952151/150‑1 (MQ=32) ccGCGCAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGttcttc < 2:3776143/150‑1 (MQ=32) gcgcCACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCATCAGCGGTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa > 1:660193/1‑150 (MQ=255) gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa > 2:2401137/1‑150 (MQ=32) gcgcAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCaa > 2:1350446/1‑150 (MQ=32) cAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGaa < 2:2172178/150‑1 (MQ=32) aaCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATTAGACCCTGTGGCAGCTGTTCGAGCTGAAGACCAACAACTGGCTGTCGTTCTTCAAGAac < 2:1190906/150‑1 (MQ=12) aaCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAac < 2:5681488/150‑1 (MQ=17) aaCCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAac > 1:5084318/1‑150 (MQ=32) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacga > 1:1460681/1‑150 (MQ=32) ccGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAacga < 2:1144619/150‑1 (MQ=32) cGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgac < 1:1859622/150‑1 (MQ=17) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca < 2:2873237/150‑2 (MQ=17) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTGAAGAacgaca < 2:4572377/150‑2 (MQ=17) gtgtGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgaca < 1:3551280/150‑2 (MQ=17) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgacag < 1:1287612/150‑3 (MQ=12) tgtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACGGGCTGTCGTTCTTGAAGAacgacag < 2:5572790/150‑3 (MQ=12) gtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACgacgac < 1:3238439/150‑1 (MQ=32) gtgGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACgacgac < 2:968327/150‑1 (MQ=32) gggCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACa < 2:1621168/149‑1 (MQ=17) ggCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACaa < 1:5640741/150‑1 (MQ=32) gcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTa > 1:4400047/1‑150 (MQ=11) gcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgaca > 1:4697628/1‑142 (MQ=11) gcTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCCGTCGTTCTTGAAGAacgaca < 2:4697628/143‑2 (MQ=11) | CCGCGCAACCGTGTGGCGCTGAAGATCAAGATCAACGAAGGCACCGTCGCCGCCATCCAGCACATCAACATCGTTGGCAACAACGTATTCGATGATGAGACCCTGGGGCAGCTGTTCGAGCTGAAGACCACCAACTGGCTGTCGTTCTTCAAGAACGACGACAAGTA > NC_002947/3819454‑3819620 |
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: ‑ |