zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

bug fix: now contigs not present in the sequence dictionary are registered properly and do not cause the script to break 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3623 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
asivache 2010-06-23 20:59:38 +00:00
parent 611d834092
commit 17d2043354
1 changed files with 1 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
perl/sortByRef.pl
View File

@ -86,6 +86,7 @@ while ( <$INPUT> ) {
my $order; my $order;
if ( defined $ref_order{$contig} ) { $order = $ref_order{$contig}; } if ( defined $ref_order{$contig} ) { $order = $ref_order{$contig}; }
else { else {
$ref_order{$contig} = $n;
$order = $n; # input line has contig that was not in the dict; $order = $n; # input line has contig that was not in the dict;
$n++; # this contig will go at the end of the output, $n++; # this contig will go at the end of the output,
# after all known contigs # after all known contigs