Get 10% of triobra1 and calculate REAP /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.2 –noweb –geno 0 –recode –out trio1_bra.phasing.20   /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.1 –noweb –geno 0 –recode –out trio1_bra.phasing.10 Of these, 956761  snps are new, 5637 already exist Cálculo só com os brasileiros merge3.ped

 

Get 10% of triobra1 and calculate REAP /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.2 –noweb –geno 0 –recode –out trio1_bra.phasing.20   /projects/relatedness/plink-1.07-x86_64/plink –file trio1_bra.phasing –thin 0.1 –noweb –geno 0 –recode –out trio1_bra.phasing.10 Of these, 956761  snps are new, 5637 already exist Cálculo só com os brasileiros merge3.ped

 

Report 210213

Vcf::throw(‘Vcf4_1=HASH(0x1ce37c8)’, ‘The command “tabix -l /projects/1000genomes/integrated_call_s…’) called at /lgc/programs/vcftools_0.1.10/lib/perl5/site_perl/Vcf.pm line 2630
VcfReader::get_chromosomes(‘Vcf4_1=HASH(0x1ce37c8)’) called at /lgc/programs/vcftools_0.1.10/bin/vcf-merge line 195
main::init_cols(‘HASH(0x18d5c18)’, ‘Vcf4_1=HASH(0x18dcd30)’) called at /lgc/programs/vcftools_0.1.10/bin/vcf-merge line 277
main::merge_vcf_files(‘HASH(0x18d5c18)’) called at /lgc/programs/vcftools_0.1.10/bin/vcf-merge line 12
Traceback (most recent call last):
File “/lgc/scripts/reap.py”, line 51, in