Perl deduplication fasta sequence

normal method

#! usr/bin/perl -w
use strict;
my $input=shift;
my %hash;
open IN,"<$input";
$/=">";
while(<IN>){
    chomp;
    $hash{$_}=1;
}
foreach my $key(keys %hash){
    print ">$key";
}
close IN;

Bioseq modular approach

#!/usr/bin/perl
use Bio::SeqIO;

my $fas=shift @ARGV;
my $IN=Bio::SeqIO->new(-file=>"$fas",-format=>'fasta');
my $OUT=Bio::SeqIO->new(-file=>">New_$fas",-format=>'fasta');
my $check={};
while (my $seq=$IN->next_seq()){
    my $id=$seq->id;
    unless($check->{$id}){
       $check->{$id}=1;
       $OUT->write_seq($seq);
    }
}
$IN->close();
$OUT->close();
print "Finished!\n";

One-line command

cat cat_allsample.fa |perl -076 -ne 'chomp; print ">$_" unless $c{$_}++ '|grep -c '>'

Guess you like

Origin www.cnblogs.com/jessepeng/p/11221893.html