Hi,
I'm trying to extract a region from a VCF file using bcftools view, however no matter what I do it only extracts the whole chromosome, not the region. The chromosome name is farily complex (i.e. not just "chr1"). The command line I'm using is:
bcftools view -r 'gi|417531910|ref|NC_019471.1|,length=62722625':51709128-51711551 input.vcf.gz > output.vcf
Genome used for the VCF is the Oar_3.1 sheep reference.
I've put the chromosome ID in quotes as otherwise the pipes confuse the command line. Including the desired coordinates of the region (51709128-51711551) in the quotes doesn't help.
Any help appreciated!
Matt
•
written
by
ayh •
1 answer
You don't say what your chromosomes are actually called.
As noted in the bcftools documentation, -r takes a comma-separated list of regions. So your command line is asking for
- All of
gi|417531910|ref|NC_019471.1| - Bases 51709128 to 51711551 of
length=62722625(which probably is not the name of any of your contigs!)
If this extracts the whole chromosome as you say, then what you need to do is remove the extraneous guff from the comma up to (but not including) the colon.
If this returned no output and your chromosome really is named gi|417531910|ref|NC_019471.1|,length=62722625 then you would have a problem because comma is not a valid character in contig names in VCF files.
Log in to answer this question.
-
How to filter vcf files with bcftools and to check whether the filtering worked
written by Khaleesi95 •
-
bcftools csq creates no output
written by ommegang •
-
Extract the position from a huge vcf file
written by the_cowa •
-
BCFTOOLS view -S not working
written by hila.fridman •
-
bcftools filter doesn't filter variants based on expression
written by Adamo0_91 •
-
Extract only RNA-seq reads in a BAM file that extend into an intron
written by foxdie •
-
Conventions for sorting provenance (original URL, processing steps) in VCF file
written by hyanwong
-
Error in SAMTools mpileup
written by julia92796 •
-
error with vcfutils when building consensus sequence with samtools and bcftools 1.2
written by cheng •
-
Samtools Header From A Bacterial Genome Ba File
written by Sakti
try to create a BED file with the coordinates and use the option
-Rotherwise, try to rename the contigs: with
bcftools annotate --rename-chrs