Post date: May 12, 2020 9:56:6 PM
GOAL: study friends the same way we studied the Pando to build a control
1 - extract proba of being hets from filtered vcf file.
2 - try several thresholds
0.94
3 -
remove header from vcf file:
33499 pando_friends_variants_noHeader.vcf
80% filter - retains 33471 out of 33499 SNPs
filter vcf based on Boolean vector:
33471 pando_friends_variants_80_hets.vcf
4 - run extract proba hets
33471 friends_only_80.txt