COMANDOS DE PERL PARA SUBSTITUIR EL NUMERO AC POR EL ID

5' FILE:

$~/bin/ExtractAllSeqs.sh *5*dat  | gawk ' {if  ($1=="ID") {id=$2;  printf (">%s\n" , id) } else if (substr ($1,1,1) ~ / [actg] / ) {$NF=" "; printf ("%s%s%s%s%s%s\n",$1,$2,$3,$4,$5,$6)}}' > 5UTR.all.fa

3' FILE:

$~/bin/ExtractAllSeqs.sh *3*dat  | gawk ' {if  ($1=="ID") {id=$2;  printf (">%s\n" , id) } else if (substr ($1,1,1) ~ / [actg] / ) {$NF=" "; printf ("%s%s%s%s%s%s\n",$1,$2,$3,$4,$5,$6)}}' > 3UTR.all.fa