Skip to content

Commit

Permalink
bug fix: program halt when nothing is masked in truncated candidates
Browse files Browse the repository at this point in the history
  • Loading branch information
oushujun committed Jan 7, 2018
1 parent 662f1b6 commit 80840a6
Showing 1 changed file with 2 additions and 1 deletion.
3 changes: 2 additions & 1 deletion LTR_retriever
Original file line number Diff line number Diff line change
Expand Up @@ -682,7 +682,8 @@ if ($trunc>0){
`awk '{print \$1"\\t"\$1}\' $index.ltrTE.veryfalse.list > $index.ltrTE.veryfalse`;
`perl $script_path/bin/call_seq_by_list.pl $index.ltrTE.veryfalse -C $genome > $index.ltrTE.veryfalse.fa`;
`cat $index.ltrTE.stg2 $index.ltrTE.veryfalse.fa > $index.ltrTE.mask.lib`;
`${repeatmasker}RepeatMasker -e ncbi -q -pa $threads -no_is -norna -nolow -div 40 -lib $index.ltrTE.mask.lib -cutoff 225 $index.ltrTE.trunc > /dev/null 2>&1`;
my $info=`${repeatmasker}RepeatMasker -e ncbi -q -pa $threads -no_is -norna -nolow -div 40 -lib $index.ltrTE.mask.lib -cutoff 225 $index.ltrTE.trunc`;
`cp $index.ltrTE.trunc $index.ltrTE.trunc.masked` if $info=~/No repetitive sequences were detected/;

#file checking
die "ERROR: RepeatMasker is not running properly!
Expand Down

0 comments on commit 80840a6

Please sign in to comment.