| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 720,464 | 0 | C | G | 42.9% | 5.0 / 12.8 | 14 | A22G (GCC→GGC) | phrB | deoxyribodipyrimidine photolyase, FAD‑binding |
| Reads supporting (aligned to +/- strand): ref base C (3/5); new base G (6/0); total (9/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.10e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACA > REL606/720389‑720538 | tCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctgc < 2:4631545/76‑1 (MQ=255) cAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctgcc < 2:3900330/76‑1 (MQ=255) aGGAGTTTTATGATTCCCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCCCGCACTGGctgccg < 2:3534175/76‑1 (MQ=255) aGGAGTTTTATGATATCCCATCTGGTCTGGTTTCGCCAGGATTT‑TCGACTGCACGATAATCTCGCACTGGctgccg < 2:4133659/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctgccg‑cct < 2:873054/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctgccg‑cct < 2:51718/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctgccg‑cct < 2:4050376/76‑1 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctggcg‑cccg > 1:112203/1‑74 (MQ=255) ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctgccg‑cctgc < 1:1472281/76‑1 (MQ=255) aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGctggcg‑gccgcccgc > 1:3113338/1‑73 (MQ=255) ttccccTTCTGTTTTGTTTCCGCAAGATTTT‑ACGTCTGCACGATAACCCCCCCCGGGCTCCCC‑CCTGCCGC‑AATTc < 2:2038084/73‑1 (MQ=255) tccccGTCTGGTGTGGTTTCCCCGGAATTC‑ACGTGTGCACGATAACCCCGCCCTGGCTGCCC‑CCTGCCGC‑AATTCg < 2:3757292/74‑1 (MQ=255) tACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCTGCCGC‑AATTCg > 1:1270993/1‑76 (MQ=255) aCCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGGCCCC‑AATTct > 1:1761830/1‑75 (MQ=255) aCCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑CAATCGt > 1:2190457/1‑76 (MQ=255) ccATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGGG‑CACGCCGC‑TGTTCGTCt > 1:3130503/1‑76 (MQ=255) ccATCGGGTTGGGTTTCGCCAGATTTT‑GCCTCTGCACGATACTCTCGCCCTGGCTCCCC‑CCTGCCGC‑AATTCGTCt < 2:3467143/76‑1 (MQ=255) cATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑GCCGCCGC‑CCATCGTCtc > 2:3566699/1‑75 (MQ=255) aTCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTGCCGCGAATTC‑TCTGc > 1:4303251/1‑76 (MQ=255) aTCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑AAATCGGCTGc > 1:2560469/1‑76 (MQ=255) ttggtctgttTCCGCAGGTATTT‑ACTCCGGCACGATAATCTCGCCCGGGCTGCCC‑CCTGCCGC‑AATTCGTCTGCAc < 2:139109/75‑1 (MQ=255) tgggtcttgtTTCGCCTGGTTTT‑CCGCCGGCACGATACTCTCCCAGTGGTTGCCG‑CCTGCCGC‑AATTCGTCTGCAc < 2:2980500/74‑1 (MQ=255) ggtTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑AATTCGTCTGCACGCGTGcg > 1:4348914/1‑75 (MQ=255) ggtTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCCCGCC‑AAATCCTCCGCACGCGTggc > 1:3408206/1‑74 (MQ=255) ttCGCCGGTTTTT‑CCGTCTCCACGATAATCTCGCCCTGGCTGCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGc < 1:1134294/76‑1 (MQ=255) cGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑CATTCGTCTTCCCGCGGGCTGGGGt > 1:3949547/1‑76 (MQ=255) gCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGGG‑CCCGCCCC‑CAATCGTCCTCCCCCGTGCTGGCggt > 2:3158461/1‑74 (MQ=255) ccAGGATTTCACGTCTGCAGGATAACCCGGCACTGGCGGCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGtt < 2:4904104/76‑1 (MQ=255) ccAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑CAATCGTCTTCACGCGTGCTGGCggtt > 1:3922735/1‑73 (MQ=255) agtaTTG‑ACGTCTCCACGATATTCGCGCCCTGGCGCCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTtata < 1:1056155/73‑1 (MQ=255) aGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑AATTCGTCTGCACGCGTGCTGGCGGtatt > 1:156073/1‑75 (MQ=255) gttttt‑GCGTCTCCACAATACTCTGCCAGTGGCTGCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTtatat < 2:3728770/73‑1 (MQ=255) gtaTTT‑CCGTCGGCACATTAATCTCGCACGGGCGCCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTtatat < 2:1029754/74‑1 (MQ=255) ggATTG‑CCGTCTCCCCGATATCCTCCCCCGGGCCCCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTtatat < 1:4119510/76‑1 (MQ=255) gATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTGCCGC‑CATTCGTCTGCACGCGCGCGGGCGGtatata > 1:294473/1‑76 (MQ=255) tttt‑CCGTCGGCACAAAAATCTCCCCCTGGCTGCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTatatat < 2:2310953/75‑1 (MQ=255) t‑ACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCGCCGC‑CATTCGTCTCCACGCGAGCTGGCGTTCTATATctc > 1:3268521/1‑74 (MQ=255) gCACGATAATCTCGCACTGGCTGGCG‑CCCGCCCC‑AATTCGTCCTCACGCGTGGCGGCGGTTTTTATTCCTAcaccc > 2:1333480/1‑75 (MQ=255) gCACGATAATCGCGCCCGGGCCCCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacca < 1:2790638/76‑1 (MQ=255) gCACGATAACCCCGCCCTGGCCGCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacca < 1:674840/76‑1 (MQ=255) cGATAATCTCGCACTGGCTGGCG‑CCCGCCCC‑CATTCGTCTGCCCGCGGGCTGGCGGTATTTATCGCTAcaccaccc > 2:304346/1‑74 (MQ=255) aTAACCCCGCACTGGCCGCCG‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCa < 2:2518658/76‑1 (MQ=255) aaTCTCGCACTGGCTGGCG‑CCCGCCGC‑CAATCGTCTGCACGCGGGCTGGGGGTTTTTATTGGTACACCCCGcccgg > 2:4218084/1‑73 (MQ=255) ccTGGCCCCCC‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACg < 1:3876927/75‑1 (MQ=255) cTGGCTGCCG‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc > 1:4109458/1‑76 (MQ=255) cTGGCTGCCG‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc > 1:2341546/1‑76 (MQ=255) tGGCGGCCG‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCa < 2:1253404/76‑1 (MQ=255) cccccg‑cctgccgc‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa < 1:682115/73‑1 (MQ=255) tgccg‑cctgccgc‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 2:3421465/1‑76 (MQ=255) tgccg‑cctgccgc‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 2:3238839/1‑76 (MQ=255) tgccg‑cctgccgc‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 1:3338017/1‑76 (MQ=255) gccg‑cctgccgc‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 2:1443114/1‑76 (MQ=255) | TCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTT‑ACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCTGCCGC‑AATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACA > REL606/720389‑720538 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |