Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,184,591 | 0 | G | C | 34.6% | 34.6 / 20.6 | 26 | intergenic (‑201/‑211) | mltF/purL | membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase |
Reads supporting (aligned to +/- strand): ref base G (7/10); new base C (6/3); total (13/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.11e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGC > NC_002947/1184444‑1184733 | ttGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCcgcgcg > 1:932705/1‑150 (MQ=255) ttGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGATTGAGCTCTCGATCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCcgcgcg < 1:1091285/150‑1 (MQ=38) gggTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGa < 2:1662827/150‑1 (MQ=255) gTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCgctcgctc < 1:1755902/150‑1 (MQ=255) gTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgc < 1:15013/150‑1 (MQ=255) ccATCAGGGCTGCCATGTAGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTTTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCg < 1:707769/149‑1 (MQ=255) tGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcg < 2:932705/150‑1 (MQ=255) gaTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcggg > 1:110640/1‑65 (MQ=255) gaTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcggg < 2:110640/65‑1 (MQ=255) gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGa > 2:920072/1‑71 (MQ=255) gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGa < 1:920072/71‑1 (MQ=255) aGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAc > 1:133179/1‑150 (MQ=255) gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcg < 2:456400/97‑1 (MQ=255) gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcg > 1:456400/1‑97 (MQ=255) tCCGCCTTAAGATATTCCGCGCCTATGGGGTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGTACAGTGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTAACCCCGAAGTCCGGtt < 1:1352579/150‑1 (MQ=38) cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCg > 1:1753288/1‑105 (MQ=255) cGCCTTAAGATACTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCg < 2:1753288/105‑1 (MQ=255) tGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTag > 2:836844/1‑150 (MQ=255) tGAGATCGAGCGCCGCCCGCGCGGCGCTCAATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa < 2:1623017/106‑1 (MQ=37) tGAGATCGAGCGCCGCCCGCGCGGCGCTCAATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa > 1:1623017/1‑106 (MQ=37) gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTaga > 1:1615008/1‑150 (MQ=255) tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa > 1:1441215/1‑101 (MQ=255) tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa < 2:1441215/101‑1 (MQ=255) tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGcc > 2:1521439/1‑86 (MQ=21) tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGcc < 1:1521439/86‑1 (MQ=21) gAGCGCCGCCCGCGCGGCGATCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt > 2:1034750/1‑150 (MQ=255) aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTg > 2:1613455/1‑150 (MQ=255) gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGc > 1:406443/1‑150 (MQ=255) gcgcCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt < 1:1512113/65‑1 (MQ=25) gcgcCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt > 2:1512113/1‑65 (MQ=25) | TTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGC > NC_002947/1184444‑1184733 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |