Merge Samtools and Freebayes vcf files · Biostar

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

Merge Samtools and Freebayes vcf files

Hi,

I am analyzing presence of variants in exonic regions of 10 genes. I have sequencing data for samples obtained from different individuals. I aligned the NextSeq PE 150 sequencing data to hg38 human genome assembly using bowtie2 and called variants using samtools and freebayes and filtered variants for exonic regions using bedtools intersect.

I observed some common variants between samtools and freebayes vcf and unique variants in both files. I visualized these variants from both files using IGV which confirms their presence.

Can I merge these two variant files obtained from different tools? If so how can I merge the vcf file?

alignment exome variants

• 2,886 views • link

updated 7.7 years ago by Pierre Lindenbaum

• written 7.7 years ago by kspata

1 answer

gatk CombineVariants https://software.broadinstitute.org/gatk/documentation/tooldocs/3.8-0/org_broadinstitute_gatk_tools_walkers_variantutils_CombineVariants.php with the option genotypemergeoption

UNIQUIFY
    Make all sample genotypes unique by file. Each sample shared across RODs gets named sample.ROD.
PRIORITIZE
    Take genotypes in priority order (see the priority argument).
UNSORTED
    Take the genotypes in any order.
REQUIRE_UNIQUE
    Require that all samples/genotypes be unique between all inputs.

• 1 views • link

written 7.7 years ago by Pierre Lindenbaum

Log in to answer this question.

More posts like this

Merge bam files and get VCF file
written 6.9 years ago by khatami.mahshid •

I have 2 sorted Bam files and I want to merge them and get a VCF file from them using samtools pileup. first i merge …
Variant discrepancy between technical replicates
written 7.3 years ago by kspata

I am performing alignment to reference and variant calling for same sample using freebayes. Both the technical replicates were sequenced on MiSeq PE 250. The …
Set minimum variant frequency while calling variants with samtools
written 7.5 years ago by kspata

Hi All, Currently I am using following command to call variants using samtools/bcftools. samtools mpileup -d 1000000 -L 1000000 -uf ref.fasta sample.sorted.bam | bcftools view …
Remove low quality base calls from reads
written 7.5 years ago by kspata

Hi All, I was performing alignment of viral vector sample sequenced on MiSeq PE 250. I trimmed the raw reads using trim galore command as …
Missing sequence from a cosmid de novo assembly
written 7.7 years ago by kspata

Hi All, I performed de novo assembly for a cosmid sequenced on NextSeq PE 300 using SPADES. The pipeline i used is as follows: 1.Trim …
Obtaining depth for all positions for different exonic ranges from mpileup file
written 7.7 years ago by kspata

Hi, I am analyzing 96 samples sequenced using NextSeq PE 150 reads. The goal is to study variants and per exon coverage of 10 genes. …
Strategy to detect deletions in the PE 300 reads
written 7.9 years ago by kspata

Hi All, I have PE 300 sequencing data of a circular plasmid sample. I want to detect deletions which occur at low frequency with deletion …
Low coverage for bases near the end of target reference sequence
written 8.0 years ago by kspata

Hi All, I have viral samples sequenced using NextSeq PE 150 with over 90 million reads, I aligned these reads after trimming adapters and low-quality …
Call missing variants in VCF as reference allele
written 8.8 years ago by olavur

I have many VCFs from different samples. If I merge these into a single VCF using `vcftools` (`vcf-merge`), the samples where a variant wasn't called …
Combine single sample vcfs together
written 9.7 years ago by sumudu_rangika •

Hi, I have 8 bam files ( obtained from 8 isolates from same species) and need to call variants. In order to get a final …