This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Printing Reverse Complement Of Dna In Single-Line Perl

I want to write a quick single-line perl script to produce the reverse complement of a sequence of DNA. The following isn't working for me, however:

$ cat sample.dna.sequence.txt | perl -ne '{while (<>) {$seq = $_; $seq =~ tr /atcgATCG/tagcTAGC/; $revComp = reverse($seq); print $revComp;}}'

Any suggestions?

perl bash sequence
• 18,618 views • link
updated by winter_li • winter_li • written by soosus • soosus

The simpler version is: cat test.txt | perl -pe 'chomp;tr/ACGTacgt/TGCAtgca/;$_=reverse."\n"'

• 65 views • link
written by lh3 lh3

an even simpler version is cat test.txt | perl -pe 'tr/ACGTacgt/TGCAtgca/;reverse'

• 0 views • link
written by Jorge Amigo Jorge Amigo

6 answers

Why use single line. keep the code to your script fold and use these routine script everyday.

Like: perl /home/usr/bin/rcdna.pl ATCGATGCG

#!/usr/bin/perl
use strict;
my $dna=shift @ARGV;
my $rcdna= & reverse_complement_IUPAC($dna);
print "$rcdna\n";

sub reverse_complement_IUPAC {
        my $dna = shift;

        # reverse the DNA sequence
        my $revcomp = reverse($dna);

        # complement the reversed DNA sequence
        $revcomp =~ tr/ABCDGHMNRSTUVWXYabcdghmnrstuvwxy/TVGHCDKNYSAABWXRtvghcdknysaabwxr/;
        return $revcomp;
}
• 0 views • link
updated by Ram Ram • written by Shicheng Guo Shicheng Guo

shell script? you don't need perl:

tr -d "\n " < input.txt | tr "[ATGCatgcNn]" "[TACGtacgNn]" | rev
• 106 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Say I want to do it with perl, for illustrative purposes. How do I go about that? And why doesn't my code work?

• 0 views • link
written by soosus • soosus

you need to concatenate all the strings in the while lool. Once everything is read, call tr , reverse and print.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Here is the modified solution that can handle IUPAC degeneracies:

tr -d "\n" < input.txt | tr "[ATGCUatgcuNnYyRrSsWwKkMmBbDdHhVv]" "[TACGAtacgaNnRrYySsWwMmKkVvHhDdBb]" | rev
• 0 views • link
updated by Ram Ram • written by Vladimir Mikryukov • Vladimir Mikryukov

If you want to use Perl:

$ echo "CAAT" | perl -le '{$seq = <STDIN>; chomp $seq; $seq =~ tr /atcgATCG/tagcTAGC/; $seq = reverse($seq); print $seq; }'
ATTG

STDIN is a default handle, so you could just take that out for brevity. I like being explicit if it isn't too wordy, though.

• 118 views • link
updated by Ram Ram • written by Alex Reynolds Alex Reynolds

for this matter I usually go for this compact yet self-explanatory code inside my scripts

echo -n "CAAT" | perl -pe '{tr/atcgATCG/tagcTAGC/; $_=reverse}'

although I would agree that if reverse complementing is the only thing you want to do, you'll probably find shell scripting faster.

• 0 views • link
written by Jorge Amigo Jorge Amigo

This will give you the reverse complement of each line, try something like:

$tail -n +2 seq.fa | tr -d "\n" |tr "acgtACGT" "tgcaTGCA" |rev
• 0 views • link
updated by Ram Ram • written by Asaf Asaf 
perl -0777pe's/\n //g; tr/ATGCatgcNn/TACGtacgNn/; $_ = reverse $_;' sample.sequence.txt
• 0 views • link
updated by Ram Ram • written by soosus • soosus

Hi , maybe you can write a perl script like the following ,

  use Bio::Seq;

  $seqobj = Bio::Seq->new(-seq => $seq);
  $Query = $seqobj->revcom()->seq;

Good luck !

• 0 views • link
written by winter_li • winter_li

Log in to answer this question.

More posts like this