Re^10: how to read input from a file, one section at a time?

Hi,

My apologies for not being clear. Just to let you know, multiple proteins can have different header sequences but identical sequence information. When I say duplicate entries, I mean the actual sequence (and not the header). I want the script to read the input file and identify if there are more than one entries with the same sequence information and print them. Does that help? Again, sorry for the confusion and thank you for your help.

Comment on Re^10: how to read input from a file, one section at a time?

Replies are listed 'Best First'.
Re^11: how to read input from a file, one section at a time? by poj (Abbot) on Apr 02, 2019 at 15:43 UTC
Try `my %fasta_seen; FASTA_RECORD: while ( my $para = <$PROTFILE> ) { # Remove fasta header line if ( $para =~ s/^>(.)//m ){ $name = $1; }; # Remove comment line(s) $para =~ s/^\s#.*//mg; next FASTA_RECORD if $fasta_seen{ $para }++; …` [download] This may not be a sensible solution if your sequences are very long in which case consider using a message digest like Digest::MD5 poj	[reply] [d/l]
Re^12: how to read input from a file, one section at a time? by davi54 (Sexton) on Apr 02, 2019 at 16:02 UTC
And how do I print the duplicate entries?	[reply]
Re^13: how to read input from a file, one section at a time? by poj (Abbot) on Apr 02, 2019 at 16:18 UTC
`#next FASTA_RECORD if $fasta_seen{ $para }++; if ( $fasta_seen{ $para }++ ){ print "DUPLICATE : $name \n $para\n"; next FASTA_RECORD; }` [download] poj	[reply] [d/l]
Re^14: how to read input from a file, one section at a time? by davi54 (Sexton) on Apr 02, 2019 at 17:40 UTC
Re^15: how to read input from a file, one section at a time? by poj (Abbot) on Apr 02, 2019 at 17:43 UTC
Re^15: how to read input from a file, one section at a time? by davi54 (Sexton) on Apr 02, 2019 at 17:54 UTC
Some notes below your chosen depth have not been shown here
Re^14: how to read input from a file, one section at a time? by davi54 (Sexton) on Apr 02, 2019 at 16:38 UTC
Re^15: how to read input from a file, one section at a time? by poj (Abbot) on Apr 02, 2019 at 17:21 UTC
Some notes below your chosen depth have not been shown here
Re^15: how to read input from a file, one section at a time? by davi54 (Sexton) on Apr 02, 2019 at 17:30 UTC


Think about Loose Coupling
	PerlMonks