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,
• 222 views
•
link
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.
• 0 views
•
link
Log in to answer this question.