Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 720,464 | 0 | C | G | 35.7% | 31.5 / 21.2 | 28 | A22G (GCC→GGC) | phrB | deoxyribodipyrimidine photolyase, FAD‑binding |
Reads supporting (aligned to +/- strand): ref base C (6/12); new base G (10/0); total (16/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.83e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT > REL606/720389‑720539 | tCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gc < 1:1702153/76‑1 (MQ=255) tCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gc > 1:5352157/1‑76 (MQ=255) aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcg > 1:6368165/1‑76 (MQ=255) aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcg > 1:2719428/1‑76 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gccgcct < 2:3164485/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gccgcct < 2:40613/76‑1 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gccgcctg > 1:3366906/1‑76 (MQ=255) ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gccgcccgc > 1:3255626/1‑76 (MQ=255) tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggggcccgcc > 1:751865/1‑76 (MQ=255) tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgcc > 1:4002430/1‑76 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggggcccgccc > 2:686806/1‑75 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgccc > 1:4657933/1‑75 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑gccgcctgccg < 1:3433109/76‑1 (MQ=255) tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcggctgccg‑c > 1:6518542/1‑76 (MQ=255) tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgcctggcg‑c > 1:3375856/1‑76 (MQ=255) tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgcctgccg‑c > 1:5870034/1‑76 (MQ=255) tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgcccc > 1:6486091/1‑74 (MQ=255) aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgcccgccccc > 2:3428197/1‑73 (MQ=255) gATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGct‑ggcgccccccg‑ccaa > 2:2547457/1‑73 (MQ=255) gATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑Caaa > 1:3258353/1‑75 (MQ=255) tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGGCTGGCC‑CAAATct > 2:1983135/1‑75 (MQ=255) cccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGGCCGCCC‑CCACTCGTc > 2:6807901/1‑76 (MQ=255) aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGGGCCTGCCG‑CACATCGCCTGc > 2:5384130/1‑76 (MQ=255) aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGGCG‑CCCATCGTCtcc > 2:3089781/1‑74 (MQ=255) tctggtctggTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCCGCC‑GCCGCCAATTC‑TCTGcc > 1:6528923/1‑75 (MQ=255) tgggtttgttTTCGCCAGGATTTACGTCGCCACGATAACCTCGCACTGGCT‑GCCCCCTGCCG‑CAATTCGTCTGCAc < 2:2012248/74‑1 (MQ=255) tggtttggttTCGCCGGGATTTCCGTCTCCAAGATACTCTCGCACGGGCT‑GCCCCCTGCCG‑CAATTCGTCTGCAcg < 1:3150951/76‑1 (MQ=255) gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CCATTCGTCTGGAcccg > 1:2413350/1‑76 (MQ=255) tctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCC‑CCACTCCTCTGCCcgccg > 2:3381910/1‑74 (MQ=255) tgtgtTCCGCCAGTTTTCACGTCGCCAGGATATTCTCCCACGGGCT‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTg < 2:5040985/73‑1 (MQ=255) ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CCAATCGTCTGCACGCGTg > 1:821112/1‑76 (MQ=255) ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCTGCCC‑CCAGTCGTCTCCAAGCGTGcc > 1:2162677/1‑75 (MQ=255) ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CCATTCGTCTGCACGCGTGCt > 1:470125/1‑76 (MQ=255) ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCCGCCG‑CCATTCGTCGGCGCCCGCGCt > 1:6919666/1‑76 (MQ=255) ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCCGCCG‑CAATTTGTCTGCACGCGCGCt > 1:1676719/1‑76 (MQ=255) gtTTCGCCAGGTTTTCCGTCTGCACGATAATCCCGCCCGGGCT‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTg < 1:6485096/76‑1 (MQ=255) ttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGGGCCCGCCG‑CCATTCGTCTGCACGcgggcgggg > 2:979495/1‑75 (MQ=255) gcGCCAGTTTTTGCGTCTGCCCGATATCCGCCCCCTGGCG‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCg < 1:3281857/75‑1 (MQ=255) cgcgAGGTTTTCCGTCTGCACGATAATCCGGCGCGGGCG‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGt < 2:2880593/76‑1 (MQ=255) cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCCCCC‑CCAATCGTCTTGCCCCCTGCTGGCgg > 2:3365907/1‑75 (MQ=255) gcgaGGTTTTCCGTCTCAACTATAATCCCCCCGTGGCG‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGtt < 1:5735316/73‑1 (MQ=255) gCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCC‑CCAATCGTCTGCCCGCCCGCTGGCGtt > 2:4999835/1‑76 (MQ=255) cAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CCATTCGTCTGCACCCGTGCCGGCGTtat > 1:4577348/1‑76 (MQ=255) cAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCCCGCCCCG‑CCCTTCGTCTTCCCGCCTGGTGGGGttct > 2:774770/1‑74 (MQ=255) aGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCC‑CCATTCGTCCTCCCGCCCGCTGGGGTtata > 2:1993987/1‑76 (MQ=255) ggttttACTTCTGCAGGATACTCTCGCGCGGGCT‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTtatat < 2:3758516/73‑1 (MQ=255) gATTTACGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCCCCCC‑GAAATCTCCTGCCCCCGTGGTGGCGGgatata > 1:4346013/1‑76 (MQ=255) gATGTACGTCTCCACAATAATCTCGCACTGCCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTtatata < 2:4838402/76‑1 (MQ=255) aTTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CCATTTCTCCCCCCGCGGGCTGGGGTTatatat > 1:1784109/1‑76 (MQ=255) aTTTACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CAATTTGTCTGCACGCGTGCGGGGGGTttatat > 1:2604087/1‑76 (MQ=255) tttACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCTCCC‑CCATTCGTTTGCCCGCGCGCTTGCGGtatattta > 1:2662071/1‑73 (MQ=255) tttACGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCC‑CCATTCGTCTTCCCGCCGGCTCGGGGTTTATATc > 1:2116265/1‑76 (MQ=255) ttACGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCG‑CAATTCGTCTGGACGCGGGCTGGCGTTATATATCg > 1:4289956/1‑76 (MQ=255) tgccGTCCGCACGATAATCCCCCGCTGCCC‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCg < 1:1991376/73‑1 (MQ=255) tACGTCTGCACGATAATCTCGCACTGGCT‑GCGGCCTGCCG‑CCAATTGTCTGCACGCGCGCTGGCGTTATATATAGc > 1:3248194/1‑76 (MQ=255) aCGTCTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCC‑CACATTGTCTTCCCCCGTGCTCGCGTTTTTTATTGCt > 2:76577/1‑76 (MQ=255) ggTCGGCACGATACTCGCCCCGTGGCG‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTa < 2:1289838/75‑1 (MQ=255) gTCTGCAAGATAATCTCGCACGGGCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTac < 2:6011454/76‑1 (MQ=255) tctcCACGATAATCCCCCACGGGCT‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTaca < 2:4650986/76‑1 (MQ=255) tCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCG‑CCATTCGTCTTCACGCGGGCTGGGGTTATATATCGCTacc > 1:5168460/1‑75 (MQ=255) cTGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCG‑CAAATCGTCTGCCCGCGCGCCGGGGGTATATATCTCTacac > 2:4948785/1‑76 (MQ=255) tGCACGATAATCTCGCACTGGCT‑GGGGGCCGCCG‑CCAATCGTCTTCACCCGGGCTGGGGCTTTATATTTCTTcacc > 2:3039536/1‑76 (MQ=255) tGCACGATAATCTCGCACTGGCT‑GGCGCCCGCCC‑CCATTTCTCCGCCCGCGTGCTGGCGGTATATATTGCTAcacc > 1:6110081/1‑76 (MQ=255) tGCACGATAATCTCGCACTGGCT‑GCCGCCTTCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc < 2:1374817/76‑1 (MQ=255) cgCACGATAATCCCGCACGGGCT‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc < 2:3213154/75‑1 (MQ=255) cACGATAATCTCGCACTGGCT‑GGGGCCCGCCG‑CCATTCGTCTGCACCCGTGCTGGCGTTATATATAGCTAcaacac > 2:2008968/1‑76 (MQ=255) cACGATAACCCCCCCCGGGCC‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcaccac < 1:4818342/76‑1 (MQ=255) aCGATAATCTCGCACTGGCT‑GGCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGGTATATATCGCTACACCACg > 1:866409/1‑76 (MQ=255) aCGATAACCCCGCCCTGCCC‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg < 1:8206/76‑1 (MQ=255) aTAATCTCGCACTGGCT‑GGCGCCCGCCC‑CCAATCGTCTGCCCCCGTGCTGGCGGTATATATTGCTTCAACACCCCa > 2:4455274/1‑76 (MQ=255) aaTCTCGCACTGGCT‑GGGGCCTCCCG‑CCCATCACCTTTCCACGTCCTGGCGTTATATATCGCTTAACAACTCCAgc > 2:4913919/1‑75 (MQ=255) aaTCTCGCACTGGCT‑GGCGCCTGCCG‑CCATTCGTCTGCACGCGGGCTGGCGTTATATATCGCTACACCAccccccg > 2:5056880/1‑73 (MQ=255) aaTCTCGCACTGGCT‑GGCGCCCGCCC‑CCACTCGTCTGCACCCGGGCCGGCGGTATATATTGCTAAACCAccccccg > 2:4552703/1‑73 (MQ=255) cccTGGCC‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGAc < 2:6600382/74‑1 (MQ=255) cTGGCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACCCCACGCCAGTGGGCGgcgc > 2:1454927/1‑76 (MQ=255) cTGGCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc > 2:2903767/1‑76 (MQ=255) cTGGCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGAAGc > 1:6044389/1‑76 (MQ=255) cTGGCT‑GCCGCCCGCCG‑CAATTCGTCTGCACGCGTGCTGGCGGTATATATCGCTACACCACGCCAGTGGGCGACGc > 2:1510522/1‑76 (MQ=255) cTGGCT‑CCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc < 2:5352157/76‑1 (MQ=255) gCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATa > 2:176335/1‑76 (MQ=255) gCT‑GCCCCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATa < 2:5481310/76‑1 (MQ=255) ggcgccccctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATa < 2:3972773/73‑1 (MQ=255) ct‑gccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa > 1:6355224/1‑76 (MQ=255) tgcccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACCCATaa < 2:4062741/73‑1 (MQ=255) t‑gccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 1:886991/1‑76 (MQ=255) t‑gccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 1:2857894/1‑76 (MQ=255) cccccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa < 2:2153175/73‑1 (MQ=255) cgccccctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc < 2:1502009/75‑1 (MQ=255) gccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 2:3289017/1‑76 (MQ=255) gccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 1:2266882/1‑76 (MQ=255) ccgcctgccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt > 1:1473478/1‑76 (MQ=255) | TCAGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCT‑GCCGCCTGCCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT > REL606/720389‑720539 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |