This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to exclude a region with bcftools

Heys,

Quick question here. I want to exclude a region with bcftools, how can I do it? I am aware that I can check specific regions with a file.txt and

bcftools view -R file.txt file.vcf

How can I exclude the very same regions?

Thanks in advance!

bcftools
• 3,063 views • link
written by gubrins gubrins

1 answer 

bcftools view -T ^file.txt file.vcf

Sorry, I've just found it, just in case it is useful for my future me or others. Important to add the ^ symbol before the coordinates!

• 0 views • link
written by gubrins gubrins

Log in to answer this question.

More posts like this
  • Problem creating new conda environment
    written by gubrins gubrins
    Heys, So I am trying to install a new conda environment with an .yml file but I encountered this error: ``` ResolvePackageNotFound: - r-base==3.2.2=0 ``` …
  • Counting variants with SnpEff
    written by gubrins gubrins
    Heys, I am working with snpEff and I want to calculate the number of variants for specific categories. I am following this code: cat output.ann.vcf …
  • create vcf from bed coordinates
    written by gubrins gubrins
    Heys, I have a silly question for you, that I am not managing to accomplish. I have a vcf file and a bed file with …
  • print next line IF statement
    written by gubrins gubrins
    Heys, I'm working with a fasta file with several individuals and I want to split them by a particular coordinate in the genome. I wrote …
  • How much weights a million reads
    written by gubrins gubrins
    Heys, I have a quick question: I was wondering how much weights a million Illumina reads of let's say 100 bp long in Gb. May …
  • bcftools multiple commands
    written by Filago • Filago
    Hello, I want to process a VCF-file with bcftools. Amongst other things I want to: 1. use multiple commands in one line of code e.g. …
  • bcftool, find a row with type BND
    written by gdm01 • gdm01
    Hello, I am trying to learn how to use bcftools --targets and --regions options. By specifiying some start and end positions in my linux command, …
  • VEP repeated region flag
    written by cocchi.e89 cocchi.e89
    Quick question, I googled around but found no answer for it, is there any VEP functionality that allows to know if a variant falls in …
  • Target capture sequencing
    written by gubrins gubrins
    Heys! First of all sorry if the question is to basic, but I'm new to bioinformatics. I'm working with target capture sequencing and I'm having …
  • BCFTOOLS view -S not working
    written by hila.fridman • hila.fridman
    Hello, I'm trying to exclude some samples from my vcf. My command looks like this: bcftools view -S file.txt vcf_file.vcf >new_vcf.vcf file.txt looks like this …