Skip to content

Commit

Permalink
Support buggers who dunno how to convert hg19 to grch37
Browse files Browse the repository at this point in the history
  • Loading branch information
ckandoth committed Nov 1, 2015
1 parent e32ef37 commit 25161f8
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion maf2vcf.pl
Original file line number Diff line number Diff line change
Expand Up @@ -142,8 +142,8 @@

# Make a version of chrom without chr-prefixes, and chrM renamed to MT
my $no_prefix_chr = $chr;
$no_prefix_chr =~ s/chrM/MT/;
$no_prefix_chr =~ s/^chr//;
$no_prefix_chr =~ s/^M$/MT/;

# To represent indels in VCF format, we need to fetch the preceding bp from a reference FASTA
my ( $ref_len, $al1_len, $al2_len ) = map{( $_=~m/^(\?|-|0)+$/ ? 0 : length( $_ )) } ( $ref, $al1, $al2 );
Expand Down

0 comments on commit 25161f8

Please sign in to comment.