Working with multi-biallelic rows of VCF in plink
0
0
Entering edit mode
3.3 years ago
MAPK ★ 2.1k

Hi All, I have a VCF file that I generated after merging INDELS.vcf and SNPs.vcf. I used picard tool to merge these into one VCF as java -jar picard.jar MergeVcfs I=${SNPS} I=${INDELS} O=${BOTH}.

This generates multi-row biallelic VCF (not multiallelic). I was wondering if plink has issues with this file. I preliminary analysis shows that plink gives warnings for variants at the same site, but I am not sure if this is a problem.

Also, can someone please suggest a better way to merge INDELs with SNPs ?

plink vcf gatk • 787 views
ADD COMMENT
1
Entering edit mode

Maybe pipe the output to bcftools -m +any - that should merge multi-allelic records into single records irrespective of if they're SNP records or Indel records.

ADD REPLY

Login before adding your answer.

Traffic: 1575 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6