Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,022,633 | 0 | C | T | 29.6% | 34.1 / 16.4 | 27 | intergenic (‑328/+340) | PP_3547/emrB | NAD(P)‑binding oxidoreductase/multidrug transporter permease |
Reads supporting (aligned to +/- strand): ref base C (12/7); new base T (4/4); total (16/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCC > NC_002947/4022490‑4022775 | cAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGTATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTcc < 1:145008/148‑1 (MQ=255) cGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTAcggcg > 4:96590/1‑148 (MQ=255) cAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTAcggcgg < 2:87751/147‑1 (MQ=255) cAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTAcggcgg < 3:45995/147‑1 (MQ=255) cAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGACACCCTTCCGCCCTTACGGCggg > 1:152408/1‑148 (MQ=255) cAAACGGACGGTTCGCTGAGGCAAAACAGACGAGTGCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGGGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTAcggcgg < 7:80714/147‑1 (MQ=255) gacggTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACtt > 3:88521/1‑148 (MQ=255) ggTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACttttt < 6:176079/148‑1 (MQ=255) ggTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACttttt < 8:146246/148‑1 (MQ=255) gCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCaa > 6:132009/1‑148 (MQ=255) gTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAAcgc < 7:7252/147‑1 (MQ=255) tCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGGATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCAACTTTACGCCCTTAAGGCGTGTCAATTTTTGTCAAAcgcg > 7:147215/1‑147 (MQ=14) tCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAAcgcg > 5:130963/1‑147 (MQ=255) tCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGa < 4:88521/148‑1 (MQ=255) aaCAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACaa > 3:76871/1‑147 (MQ=255) cAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGCGTAGGTATTCATTCGCGCAGGGGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACaaaaa < 7:134029/148‑1 (MQ=255) ttcttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGcc > 3:148146/1‑148 (MQ=39) ttgcttcttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCt > 7:32231/1‑148 (MQ=39) ttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTAc > 5:9079/1‑148 (MQ=39) gggCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCggg < 6:70928/148‑1 (MQ=39) cATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGtt > 7:83806/1‑148 (MQ=39) aTCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTACCTCGc > 5:106246/1‑148 (MQ=255) ggTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCt > 3:96857/1‑148 (MQ=255) aCGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGgga > 6:5517/1‑148 (MQ=40) aCGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGgga > 4:86407/1‑148 (MQ=40) tAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAg < 7:31682/148‑1 (MQ=255) tCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGGTCCCGGGAGGACCGCGCTGCCGGcc > 6:160814/1‑148 (MQ=255) | CAGCGTCGCCATCGCAAACGGACGGTTCGCTGAGTCAAAACAGACGAGTTCGCATCTTGCTTCTTGCTTCTTGCTTCTTGCTTCTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAGGTGGCGCATAGTCACCCTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCC > NC_002947/4022490‑4022775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |