Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,530,026 | 0 | T | A | 30.8% | 27.5 / 14.8 | 26 | intergenic (+32/‑136) | PP_4018/topB | acetyltransferase/DNA topoisomerase III |
Reads supporting (aligned to +/- strand): ref base T (13/5); new base A (2/6); total (15/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.84e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGC > NC_002947/4529879‑4530167 | gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTa > 1:135524/1‑149 (MQ=255) ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCcg > 4:134386/1‑149 (MQ=255) cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCcgc < 3:134386/149‑1 (MQ=255) gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa < 2:135524/149‑1 (MQ=255) gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa < 6:24366/149‑1 (MQ=255) tGCACGCCACGCCGTTCTATGAGCGGCTGGGTTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCGCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa < 4:189237/148‑1 (MQ=255) tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCCGCGAag < 2:217709/149‑1 (MQ=39) acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAAGAgg < 6:300364/149‑1 (MQ=255) acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCGCGAAGAgg > 4:298735/1‑149 (MQ=255) ttCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCACTTCGCCGGTTTACCCGCGAAGAGCCCGGCTCAGGc > 7:82138/1‑149 (MQ=255) gCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATcc > 1:283528/1‑149 (MQ=255) tCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCt > 3:294952/1‑149 (MQ=255) ttaaTTCCTCGAAGCTGGCATCCCGCATGTCGACTTGGTGCGCTACTCGCGCGCCTTATCCATCACCTGCACTGGCCTCTTCGCGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcat < 4:66126/146‑1 (MQ=255) ccTCGATGCTGCAATCCCGCTTGTCCACTTGGGGCGCTACCCGCGCGCCTTACCCATCCCCTCCCCTGGCTTCTTCGCGGTTTAACCGGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACatccatcc < 4:294952/148‑1 (MQ=255) ccTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc > 3:69123/1‑147 (MQ=255) gacgCTGCCATCCCGCTTGTCGACTTGGGGCGCTACTCGCGCGCCTTATCAACCACCTCCATTGGCCTTTTCGCGGGTTTACCCGCGATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGa < 4:69123/144‑1 (MQ=255) ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCctct > 2:31764/1‑149 (MQ=255) cccgACTCCCGCGCCTTATCCATCACCTGAACTGCCCTCTTCGCGGTTTAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCAt < 3:298735/147‑1 (MQ=255) gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATcc > 5:115674/1‑148 (MQ=255) gcgcgcCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATg < 5:308660/147‑1 (MQ=255) ccATCACCTGCACTGGCCTCTTAGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGCCAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCAcccccc > 5:221395/1‑147 (MQ=255) cATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGCAGCGGCCGGCTCCGGCGATATACATCCCCGCTGCCAAAAAACCGGTACATCCATCCAGATAAACCTTGCCTCTGCGCCCCCGCTTGCGCCTCCTATGGCCCCCCACCCCGa > 7:255449/1‑149 (MQ=255) cTGCACTGGCCTCTTAGCGGGTTTACCCGCGAAGACGACCGCTCAGGCGATATAAAAACCCCCTGACAAAAAACCAGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCCCATCCTATTGACCACCACCCCgaagaa > 4:23464/1‑147 (MQ=255) ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGTAAACACTTGACTCTGCGCCCCCGCTTGCGCATCCTAAGGACCACCACCCCGATGAAGCGTCACCGGCCaag > 1:103149/1‑147 (MQ=255) gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc > 5:14348/1‑149 (MQ=255) gcgGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATcc > 4:177586/1‑109 (MQ=255) | GCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGC > NC_002947/4529879‑4530167 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |