Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,396,513 | 0 | C | G | 22.2% | 93.1 / 11.2 | 36 | S32S (TCG→TCC) | yjfQ | DNA‑binding transcriptional dual regulator |
Reads supporting (aligned to +/- strand): ref base C (4/24); new base G (8/0); total (12/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.64e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGATTTG > REL606/4396441‑4396587 | tGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGaa < 1:1098140/76‑1 (MQ=255) tGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGaa < 2:1658287/76‑1 (MQ=255) cgcCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAAt < 2:237357/76‑1 (MQ=255) tGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAg < 1:29390/76‑1 (MQ=255) cgcACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGAc < 1:1385557/76‑1 (MQ=255) tttttCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCaa < 2:1731370/76‑1 (MQ=255) tttCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGGAATCCCAGCAG‑CTCCACg > 1:1317749/1‑76 (MQ=255) ttCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGa < 2:1356049/76‑1 (MQ=255) tCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATTCCC‑GACG‑CCCAACGAc > 1:18618/1‑76 (MQ=255) cAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑t < 2:1433435/76‑1 (MQ=255) aGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑tt < 2:285957/76‑1 (MQ=255) gTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑ttt < 2:228310/76‑1 (MQ=255) tttGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGAATTCCCACACG‑CTCAACGAC‑TTtc > 1:60418/1‑76 (MQ=255) tGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGAAATCCCCCACGCCCCAACAAC‑ttttt > 1:1080557/1‑74 (MQ=255) ccGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAAATCCCCGACC‑CTCAACGAC‑CTtctcca > 2:1811044/1‑74 (MQ=255) cGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAAc < 2:18618/76‑1 (MQ=255) cTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGGAAATCCCCGACG‑CCCCACGAC‑TTTTTCTAccg > 2:693861/1‑74 (MQ=255) tCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGAATTCCCAGACG‑CCCAACAAC‑TTTCTCCACgggcc > 2:1823942/1‑73 (MQ=255) cGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATTCCC‑GACG‑CCCCACGAC‑TTTTTCCACAGCcac > 1:795978/1‑76 (MQ=255) cGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGAAATCCCCCACC‑CCCAACGAC‑CTTCTTTACCGTcac > 1:1371940/1‑76 (MQ=255) tCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGGAATTCCAGACG‑CTCCACCAC‑TTTCTCCACCGCCACaa > 2:1045692/1‑76 (MQ=255) cAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACaaa < 1:776960/76‑1 (MQ=255) aGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGGGGAATTCCCCGACG‑CCCCACGAC‑CTTCTCAACCGTCACaaaac > 2:1457434/1‑74 (MQ=255) tttATTGATATCGCGTCGCGCAGTGGCA‑GGGGGAAATTCCAGACG‑CCCCACCAC‑TTTTTTCACAGTGACAAAAccc > 2:797198/1‑76 (MQ=255) tttATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGccc < 1:1691717/76‑1 (MQ=255) ttATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATTCCC‑GACG‑CCCCACCAC‑CTTCTCAACAGCCACAAAACCCa > 2:463212/1‑76 (MQ=255) aTTGATATCGCGTCGCGCAGTGGCA‑GGGGGAAATTCCACACG‑CCCCACCACTTTTTTCACCGGCCCCAAACCCCaa > 2:669679/1‑76 (MQ=255) ttGATATCGCGTCGCGCAGTGGCA‑GGCGGAAATCCCAGAAG‑CCCCACGAC‑CTTTTTCACAGTCACCAAACCCaaat > 2:480807/1‑74 (MQ=255) tGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATtg < 2:1052918/76‑1 (MQ=255) gATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATtgt < 2:513073/76‑1 (MQ=255) atCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCg < 1:381458/76‑1 (MQ=255) tcgcgtcgcgCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGa < 2:1527906/76‑1 (MQ=255) cgtcgcgCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCa < 2:60418/76‑1 (MQ=255) tcgcgCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCAtt < 1:1653781/76‑1 (MQ=255) cgcAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTcc < 2:130104/76‑1 (MQ=255) gcAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCa < 2:952925/76‑1 (MQ=255) ggCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCaggagg < 1:1139940/76‑1 (MQ=255) cA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGAt > 2:1047581/1‑76 (MQ=255) a‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGAtt > 1:591547/1‑76 (MQ=255) ggCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGAttt < 1:904834/76‑1 (MQ=255) gCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGATTTg < 1:222416/76‑1 (MQ=255) | TGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCA‑GGCGAAATTCCCAGACG‑CTCAACGAC‑TTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGATTTG > REL606/4396441‑4396587 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |