Question: Merge Samtools and Freebayes vcf files
0
gravatar for kspata
9 weeks ago by
kspata40
Chicago
kspata40 wrote:

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?

variants alignment exome • 157 views
ADD COMMENTlink modified 9 weeks ago by Pierre Lindenbaum116k • written 9 weeks ago by kspata40
0
gravatar for Pierre Lindenbaum
9 weeks ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum116k wrote:

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.
ADD COMMENTlink written 9 weeks ago by Pierre Lindenbaum116k
Please log in to add an answer.

Help
Access

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