I have mapped reads to a reference genome of a related species. I want to remove reads which map to a specific region (chromosome) of the reference, but I don't know what the best way to go about it is. I'm doing variant calling, and I've been using the -t ^ option in bcftools call to exclude the relevant chromosomes, but does this target regions in the reference, or in the mapped reads? Is it the same thing? Is it better to do this during mpileup phase?
I feel completely lost and would appreciate an ELI5.
•
written
by
clinnaeus •
1 answer
use bedtools complement to get the complement of your blacklisted.bed and then use samtools view -L complement.bed -O BAM -o out.bam in.bam
or use samtools view -M -L blacklisted.bed -O BAM --unoutput out.bam in.bam > /dev/null
Log in to answer this question.
-
bcftools mpileup before bcftools call
written by Egelbets •
-
Error: RSEM Bam file for variant calling
written by kartikayprasad •
-
Indel realignment: GATK vs. BCFtools
written by martytriggers •
-
Amplicon sequencing and unwanted amplicons
written by misterie
-
Generate VCF from different .bam files with different chromosome names
written by evoclive •
-
Detect Transgene Insertion Site using SAM format
written by kspata
-
Nanopolish variant calling. Need to specify input region of the reference genome
written by SUDOsundu •
-
Does samtools view only show reads that completely span the region?
written by Noam Teyssier
-
Phasing alleles, indel gap filled with reference.
written by Vincent Manzanilla •
-
priors for calculting genotype likelihood using RNAseq data
written by prasundutta87
Not sure if this would help but
samtools viewnow has the optionYou could create intervals for parts you want to keep.