Skip to content

Commit

Permalink
subroutine of decide_h_a_b
Browse files Browse the repository at this point in the history
If a genotype can't be decided as h, a, or b, then assign it as missing
  • Loading branch information
lileiting committed Jun 3, 2015
1 parent de1dd52 commit a0c86ed
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion pileup2singledose.pl
Original file line number Diff line number Diff line change
Expand Up @@ -111,7 +111,7 @@ sub decide_h_a_b{
return 'a';
}elsif(&judge($bases,$mutant,$threshold) and !&judge($bases,$main,1)){
return 'b';
}else{die}
}else{return ".."}
}

sub variation_type {
Expand Down

0 comments on commit a0c86ed

Please sign in to comment.