More posts like this
-
Are there any aligners that enable me to align sense strands only to a sense strand reference seque…
written by b10hazard •Is there an aligner that allows me to ignore reverse complement alignments. For example, my reference sequence I'm aligning to will always be single stranded …
-
Just complement or reverse sequence fom Biopython, but not reverse-complement one!
written by natasha.sernovaDear all, I have a problem with Biopython. I give it a fasta-sequence and need to make either reversed fasta sequence out of it in …
-
How to make a frequent word counting program in python (record occurrence of every possible 7-mer)?
written by Kevin_Smith •I will like to discover the common 7-mer from all the sequences in a file (.txt , each line in the file contains a single …
-
How to find out the reverse complement of DNA from each FASTA formated sequence file in a directory…
written by Sumit •<p>I have total 2000+ genome sequence files in a directory. I need reverse complement sequence for each of the files and want to generate FASTA …
-
Merging Illumina Paired End Reads
written by vijay<p>Dear All,</p> <p>I have fastq a dataset containing forward and reverse sequences obtained through paired end module of Illumina platform. I am trying to merge …
-
Forum: Biostar 1.3.4: Post Closing And Poll Added
written by Istvan AlbertBased on the results of the polling and other feedback the post closing option has been re-enabled. One restriction is that a post can only …
-
Command Line Tool To Obtain Complementary And/Or Inverse Strand From Fasta Files
written by etienne.decarie •<p>From a fasta file, I would like to obtain three new fasta files containing:</p> <ol> <li>Inverse</li> <li>Complement</li> <li>Antiparallel</li> </ol> <p>I am looking for a command …
-
Printing Reverse Complement Of Dna In Single-Line Perl
written by soosus •<p>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, …
-
Computing The Reverse And Complement Of A Sequence With Pygr
written by Biostar User •Computing the reverse complement with the [Pygr](http://code.google.com/p/pygr/wiki/PygrDocumentation) bioinformatics framework: ```py # # Reverse complement example with pygr # from pygr.sequence import Sequence # needs a …
-
Computing The Reverse And Complement Of A Sequence With Biopython
written by Biostar User •An example that computes the reverse complement of a sequence with [BioPython][1] ```py # # Reverse complement example with BioPython # from Bio.Seq import Seq …
Please do a proper search within the forum site
http://whathaveyoutried.com/
Computing the reverse and complement of a sequence with Biopython
This is a great example of a lazy question with no evidence of an attempt to solve the problem yourself. If down-votes were still implemented at this site, it would get them.
In the absence of downvotes, any reason why we shouldn't just delete this question?
Good point. No.
testing the new post closing feature
same here, testing the comment-before-closing feature
So: this question would be/has been closed because in my opinion (1) the problem is common with a trivial, easily-found solution and (2) the user has made no attempt to find the solution themselves and is simply asking someone else to do their work for them.
These comments then serve as a record as to why the question was closed.
Simple hey :)
copy the following code to rcdna.pl
then use:
perl rcdna.pl ATCG
the script will return: TAGC
#!/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; }