This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract genotype from hg19

I have a large list of chr and positions, and would like to obtain the genotype of these positions in hg19. I saw this earlier post, but is there any easier way besides using samtools for each of the positions please?

Thank you.

hg19
• 2,003 views • link
updated by Chris Fields Chris Fields • written by DVA DVA

1 answer

seqtk has the 'subseq' command which has a bed option:

seqtk subseq [options] <in.fa> <in.bed>|<name.list>

This can generate tab-delim output as well using -t.

• 0 views • link
written by Chris Fields Chris Fields

Thanks Chris. I will give that a shot.

• 0 views • link
written by DVA DVA

Log in to answer this question.

More posts like this