This is a test version of Biostars. For the public version, visit https://www.biostars.org.
retrieving sequences of a upstream and downstream of a coordinate for hg19

I have location of peak summit from my chipseq data and want to extract sequences 500 bp upstream and downstream of this coordinate from hg19. What should be easieast approach. Apology for a naive question, but I could not recollect or find on net how I used to do this.

Thanks

chip-seq
• 2,100 views • link
updated by ATpoint ATpoint • written by kanwarjag kanwarjag

1 answer 

awk 'OFS="\t" {print $1, $2-500, $3+500}' summits.bed | samtools faidx hg19.fa /dev/stdin > sequences.fa
• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this