Extract sequence by using Ids file (assuming one ids in a line): perl -ne 'if(/^>(\S+)/){$c=$i{$1}}$c?print:chomp;$i{$_}=1 if @ARGV' ids.file fasta.file #PerlTricks #PerlOneliner #Perl #Extract #Ids #
A perl script for mapping fasta formated sequences to multiple reference sequences using one of several alignment programs. http://seq.crg.es/download/software/Miro/doc/run_Mapping.html #Perl #Script
One code to find them all is a set of perl scripts to extract useful information from RepeatMasker http://doua.prabi.fr/software/one-code-to-find-them-all