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

Split out contig names from Reference .fai file on white space (to support the GATK resource bundle's file human_g1k_v37.fasta.fai.gz, which does not use tab delimiters)

This commit is contained in:
Menachem Fromer 2012-06-07 16:56:32 -04:00
parent a90482c772
commit a13d125ba1
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
public/perl/sortByRef.pl
View File

@ -50,7 +50,7 @@ my %ref_order;
my $n = 0; my $n = 0;
while ( <DICT> ) { while ( <DICT> ) {
chomp; chomp;
my ($contig, $rest) = split "\t"; my ($contig, $rest) = split '\s';
die("Dictionary file is probably corrupt: multiple instances of contig $contig") if ( defined $ref_order{$contig} ); die("Dictionary file is probably corrupt: multiple instances of contig $contig") if ( defined $ref_order{$contig} );
$ref_order{$contig} = $n; $ref_order{$contig} = $n;