This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using tabix to subset a region from a VCF file

I've read quite a few posts that suggest I can subset my VCF file for a region of interest using Tabix.

However, when I try as below the output file is empty

VCF="/Volumes/Seagate Expansion Drive/temp/130iPSC_061118.snp.vcf.gz"
tabix -p vcf "$VCF"
tabix "$VCF" 15:31196055-31235311 > "$DIR"/vcf/sliced.vcf
tabix "$VCF" -R "$DIR"/source/regions.bed > "$DIR"/vcf/tabix.vcf

I've subsequently managed to get it working using vcftools, as follows:

vcftools --gzvcf "$VCF" --chr chr15 --from-bp 31196055 --to-bp 31235311 --recode --recode-INFO-all --out "$DIR"/vcf/sliced

But I'd still like to know how to use Tabix

tabix conda vcf
• 3,057 views • link
written by michael.flower.14 • michael.flower.14

you are also selecting a region 0 bases long with 1:17375-17375. you may also need to use chr1 instead of 1

• 0 views • link
written by cmdcolin cmdcolin

You have disagreement between the chromosome identifiers in your examples: 15:31196055-31235311 in the example that doesn't work vs chr15 in the example that does work. Tabix works to select regions just like in your example, if the region identifiers are correctly specified.

• 0 views • link
written by seidel seidel

0 answers

No answers yet.

Log in to answer this question.

More posts like this