|
|
@@ -27,14 +27,14 @@ task sister { |
|
|
|
${family_name} ${LCL7_name} 0 0 1 -9 |
|
|
|
${family_name} ${LCL5_name} ${LCL7_name} ${LCL8_name} 2 -9" > ${family_name}.D5.ped |
|
|
|
|
|
|
|
mkdir D5 |
|
|
|
mkdir VBT_D5 |
|
|
|
/opt/VBT-TrioAnalysis/vbt mendelian -ref ${ref_dir}/${fasta} -mother ${family_name}.vcf -father ${family_name}.vcf -child ${family_name}.vcf -pedigree ${family_name}.D5.ped -outDir D5 -out-prefix ${family_name}.D5 --output-violation-regions -thread-count $nt |
|
|
|
# prepare ped file, D6 |
|
|
|
echo "${family_name} ${LCL8_name} 0 0 2 -9 |
|
|
|
${family_name} ${LCL7_name} 0 0 1 -9 |
|
|
|
${family_name} ${LCL6_name} ${LCL7_name} ${LCL8_name} 2 -9" > ${family_name}.D6.ped |
|
|
|
|
|
|
|
mkdir D6 |
|
|
|
mkdir VBT_D6 |
|
|
|
/opt/VBT-TrioAnalysis/vbt mendelian -ref ${ref_dir}/${fasta} -mother ${family_name}.vcf -father ${family_name}.vcf -child ${family_name}.vcf -pedigree ${family_name}.D6.ped -outDir D6 -out-prefix ${family_name}.D6 --output-violation-regions -thread-count $nt |
|
|
|
>>> |
|
|
|
|