...
Merging multiple VCF files
We will use 'bcftools' to merge VCF files. Initially, we need to create a list of VCF files that we plan to merge. For example: while in the above directory do the following:
Code Block |
---|
for i in `ls --color=never *.vcf`; do echo $i; echo $i >> file_list.txt; done |
if you do not yet have 'bcftools' installed then run the following command to install it using conda
Code Block |
---|
conda install -c bioconda bcftools |
Now let’s create an index file for each VCF file:
...
Code Block |
---|
bcftools merge -l file_list.txt -0 -Oz -o merged_ALL.vcf.gz |
uncompress version of merged vcf
Code Block |
---|
bcftools merge -l file_list.txt -0 -o merged_ALL.vcf |
Visualising the merged VCF file in IGV
...