This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get The Dna Sequence From The Human Genome

I have hundreds of genomic coordinates that I want to retrieve their dna sequences in the human genome. e.g. chr1 11355 11391 +

Is there a standard way to get the dna sequence in the human genome. Note that I have the hg19 pre-built index but I'm not sure how to use it.

Thanks,

deleted-post

1 answer

Perhaps there is a more slick human-genome specific way to do this. However, using the bedtools function, fastaFromBED

fastaFromBED -fi hg19.fa -bed input.bed -fo output.fa

It's quite fast.

Log in to answer this question.