/
vcftools
vcftools
A great set of tools for manipulating and analyzing VCF files
This exercise should be run on your personal computer
Installing vcftools:
Instructions for installing VCFtools are here:
https://vcftools.github.io/downloads.html
Exercise:
Download this file to your desktop, unzip it, then follow instructions below
Or copy over on TACC:
vcftools --weir-fst-pop
#Navigate to the exercise directory cd ~/Desktop/vcftools #Look at the input files: ls *.vcf ls *.txt #look at the genotype data for chromosome 1 less chrI.vcf #look at the genotype data for chromosome 12 less chrXII.vcf #a list of the male samples in the VCF less male_pun.txt #a list of the female samples in the VCF less female_pun.txt #Use VCF tools to get Fst between males and females for the two VCF files vcftools --vcf chrI.vcf --weir-fst-pop male_pun.txt --weir-fst-pop female_pun.txt --out chr1_Fst_results vcftools --vcf chrXII.vcf --weir-fst-pop male_pun.txt --weir-fst-pop female_pun.txt --out chr12_Fst_results #look at the results ls -l *.fst less chr1_Fst_results.weir.fst less chr12_Fst_results.weir.fst #Use plot_Fst.R to plot the results ./plot_Fst.R #look at the resulting plots ls *Fst_plot.pdf open *.pdf #What is biological explanation for these results?
Related content
Pungitius exercise
Pungitius exercise
More like this
2020 Analysis using BEDTools
2020 Analysis using BEDTools
More like this
2021 Analysis using BEDTools
2021 Analysis using BEDTools
More like this
Visualization using Integrated Genome Viewer (IGV)
Visualization using Integrated Genome Viewer (IGV)
More like this
Bedtools tutorial -- GVA2020
Bedtools tutorial -- GVA2020
More like this
genotyping with mpileup
genotyping with mpileup
More like this