Skip to content

Commit

Permalink
exclude chrX variant from inheritance
Browse files Browse the repository at this point in the history
  • Loading branch information
thomasguignard authored Nov 13, 2023
1 parent 9eb33ef commit d96faeb
Showing 1 changed file with 9 additions and 7 deletions.
16 changes: 9 additions & 7 deletions wwwachab.pl
Original file line number Diff line number Diff line change
Expand Up @@ -2224,11 +2224,14 @@
if ($line[0]!~/X/){

if (defined $skipCaseWT && $finalSortData[$dicoColumnNbr{"Genotype-".$case}] eq "0/0"){
switch ($familyGenotype){
#Check if case/dad and case/mum inheritance are consistent
case /^_0\/0_0\/1_0\/0_/ {$dadVariant ++ ;}
case /^_0\/0_0\/0_0\/1_/ {$mumVariant ++ ;}
}
# exclude chromosome X variants due to mother/son unbalance
if ($line[0]!~/X/){
switch ($familyGenotype){
#Check if case/dad and case/mum inheritance are consistent
case /^_0\/0_0\/1_0\/0_/ {$dadVariant ++ ;}
case /^_0\/0_0\/0_0\/1_/ {$mumVariant ++ ;}
}
}
next;
}

Expand All @@ -2240,8 +2243,7 @@
case /^_0\/1_0\/0_0\/1_/ {$caseMumVariant ++;}

}
}


if ($finalSortData[$dicoColumnNbr{"Genotype-".$case}] eq "0/0" or (! defined $hashAffected{$dad} and $finalSortData[$dicoColumnNbr{"Genotype-".$dad}] eq "1/1") or (! defined $hashAffected{$mum} and $finalSortData[$dicoColumnNbr{"Genotype-".$mum}] eq "1/1") ){
$finalSortData[$dicoColumnNbr{'MPA_ranking'}] += 100;
}
Expand Down

0 comments on commit d96faeb

Please sign in to comment.