Mercurial > repos > dereeper > sniplay
diff PedToFasta/Ped2Fasta.pl @ 3:345f88a8f483 draft
Uploaded
author | dereeper |
---|---|
date | Fri, 10 Jul 2015 10:38:43 -0400 |
parents | 420b57c3c185 |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/PedToFasta/Ped2Fasta.pl Fri Jul 10 10:38:43 2015 -0400 @@ -0,0 +1,70 @@ + +#!/usr/bin/perl + +use strict; +use Getopt::Long; + +my $usage = qq~Usage:$0 <args> [<opts>] + +where <args> are: + + -i, --in <PED input> + -o, --out <Fasta output> +~; +$usage .= "\n"; + +my ($input,$out); + + + +GetOptions( + "in=s" => \$input, + "out=s" => \$out, +); + + +die $usage + if ( !$input || !$out); + + +my %IUPAC = +( + '00'=> "?", + 'AA'=> "A", + 'CC'=> "C", + 'GG'=> "G", + 'TT'=> "T", + 'AG'=> "R", + 'GA'=> "R", + 'CT'=> "Y", + 'TC'=> "Y", + 'TG'=> "K", + 'GT'=> "K", + 'CG'=> "S", + 'GC'=> "S", + 'AT'=> "W", + 'TA'=> "W", + 'AC'=> "M", + 'CA'=> "M", +); + +open(my $O,">$out"); +open(my $P,$input) or die "File does not exist"; +while(<$P>) +{ + my $line = $_; + $line =~s/\r//g; + $line =~s/\n//g; + my @infos = split("\t",$_); + my $ind = $infos[0]; + print $O ">$ind\n"; + for (my $i = 6; $i <= $#infos; $i= $i+2) + { + my $code = $infos[$i].$infos[$i+1]; + my $letter = $IUPAC{$code}; + print $O $letter; + } + print $O "\n"; +} +close($P); +close($O);