Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,396,513 | 0 | C | G | 21.9% | 117.5 / 12.3 | 41 | S32S (TCG→TCC) | yjfQ | DNA‑binding transcriptional dual regulator |
Reads supporting (aligned to +/- strand): ref base C (1/31); new base G (9/0); total (10/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.85e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGATTTGA > REL606/4396439‑4396588 | tCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCg < 1:882209/76‑1 (MQ=255) tCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCg < 1:1239338/76‑1 (MQ=255) ttGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCa < 2:733429/76‑1 (MQ=255) tGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAg < 1:183597/76‑1 (MQ=255) aCTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCt < 2:1382772/76‑1 (MQ=255) ttttttCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCa < 1:1749926/76‑1 (MQ=255) ttttCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAAc < 1:1687103/76‑1 (MQ=255) ttCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGa < 2:730847/76‑1 (MQ=255) tCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGAc < 1:454484/76‑1 (MQ=255) cAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACt < 1:1684490/76‑1 (MQ=255) aGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACtt < 2:1730321/76‑1 (MQ=255) ttGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAATTCCCAGACGCTCACAGACCTtct > 2:1590180/1‑76 (MQ=255) tGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAATTCCCAGAACCTCAACGACCTtctc > 2:186622/1‑76 (MQ=255) tGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATCCCCCACGCGCCACGACTTtctc > 2:1335328/1‑76 (MQ=255) tGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTtctc < 1:1208059/76‑1 (MQ=255) gCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATCCCCGACGCTCCACGACCTTCTCa > 1:386258/1‑76 (MQ=255) ccGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCaa < 2:158969/76‑1 (MQ=255) cGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATGCCCGACGCTCAACGCCCTTCTCAAc > 2:1236430/1‑76 (MQ=255) cTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATCCCCCACGCCCCACGACCTTCTCAAccg > 1:414877/1‑74 (MQ=255) tCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACaa < 1:55290/76‑1 (MQ=255) tCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACaa < 1:74123/76‑1 (MQ=255) tCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACaa < 1:1189439/76‑1 (MQ=255) gTTTATTGATATCGCGTCGCGCAGTGGCAGGGGGAAATCCCAGACGCCCCACGACCTTCTTAAACGGCACAAAGcc > 2:538752/1‑76 (MQ=255) tATTGATATCGCGTCGCGCAGTGGCAGGGGGAATTCCCCGACGCCCCACCACTTTTTCCACGGGCACCAAAcccca > 2:1722251/1‑74 (MQ=255) aTTGATATCGCGTCGCGCAGTGGCAGGGGGAAATCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAACCCCAt > 1:527915/1‑76 (MQ=255) ttGATATCGCGTCGCGCAGTGGCAGGGGGAAATTCCAGACGCTCAACGACTTTCTCAACCGTCACACAACCCCAtt > 1:1581299/1‑76 (MQ=255) tGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATtg < 1:458878/76‑1 (MQ=255) atatCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATtgtg < 2:1260306/76‑1 (MQ=255) tatCGCGTCGCGCAGTGGCAGGGACAATTCCCCGACACTCCACCACTTTCTTCACGGGCACCAAACCCAAATTTGc > 2:1684044/1‑76 (MQ=255) atCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCg < 2:1043339/76‑1 (MQ=255) atCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCg < 1:948631/76‑1 (MQ=255) cgcgtcgcgCAGTGGCAGGGGGAAATTCCAGACGCGCAACGACCTTCTCATCGCGCACCAAACCCAACTTTGCggg > 2:875424/1‑74 (MQ=255) gcgtcgcgCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGCGCGAGc < 2:919558/76‑1 (MQ=255) cgtcgcgCAGTGGCAGGGGGAAATTCCAGACGCCCCACCACTTTTTCTACGGTCACCAAGCCCACATGTGCGAGCa > 2:970348/1‑76 (MQ=255) tcgcgCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCAtt < 2:1534789/76‑1 (MQ=255) tcgcgCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCAtt < 2:42906/76‑1 (MQ=255) cgcgCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCAttt < 1:600639/76‑1 (MQ=255) gcgcAGTGGCAGGCGAAATTCNCAGACGCTCAACGACTATCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTc < 2:283331/76‑1 (MQ=255) gcgcAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTc < 2:260160/76‑1 (MQ=255) gcgcAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTc < 2:353211/76‑1 (MQ=255) gcgcAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTc < 2:87264/76‑1 (MQ=255) gcgcAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTc < 1:1141194/76‑1 (MQ=255) aGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCagg < 1:97290/76‑1 (MQ=255) tGGCAGGCGAAATTCCCAGACGCTCAACGGCTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCaggag < 1:248975/76‑1 (MQ=255) cAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGAt < 1:1711046/76‑1 (MQ=255) aGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGAtt > 1:228276/1‑76 (MQ=255) cGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGATTTGa < 1:1718276/76‑1 (MQ=255) | TCTGCGCCATTGCGCACTTTTTTCAGTTTGCCGCTTTCGTCAAGTTTATTGATATCGCGTCGCGCAGTGGCAGGCGAAATTCCCAGACGCTCAACGACTTTCTCAACGGTCACAAAGCCCAATTGTGCGAGCATTTCCAGGAGGATTTGA > REL606/4396439‑4396588 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |