Question: Merge Samtools and Freebayes vcf files
gravatar for kspata
2.1 years ago by
kspata70 wrote:


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?

variants alignment exome • 1.2k views
ADD COMMENTlink modified 2.1 years ago by Pierre Lindenbaum133k • written 2.1 years ago by kspata70
gravatar for Pierre Lindenbaum
2.1 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum133k wrote:

gatk CombineVariants with the option genotypemergeoption

    Make all sample genotypes unique by file. Each sample shared across RODs gets named sample.ROD.
    Take genotypes in priority order (see the priority argument).
    Take the genotypes in any order.
    Require that all samples/genotypes be unique between all inputs.
ADD COMMENTlink written 2.1 years ago by Pierre Lindenbaum133k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2564 users visited in the last hour