Beefy Boxes and Bandwidth Generously Provided by pair Networks
The stupid question is the question not asked

comment on

( #3333=superdoc: print w/replies, xml ) Need Help??

"I would like to pull out all of the full length CDS with some 5' and 3' UTR info for the proteins I have found in the tissue by mass spec"

If I understand you correctly, I will look at the gff/gtf coordinates for the CDs feature and try to obtain that for the specific protein annotation in the gff file, A fastA file may not have these features handy so it maybe not the way to go unless you have an annotated genome somewhere for you to extract this stuff from.

"I have tried simply blasting the peptides against the assembly but I get hits that are not in open reading frames"

Maybe because BLAST is a local alignment and search tool, that behavior is totally expected.

" I would suppose the major difference is that I need to return all the ORFs from tens of thousands of entries as a new fasta file."

Read the file one line at a time. Maybe something like what follows ?... If you don't have BioPerl installed then certainly there are many other examples on reading fastA file that don't use BioPerl around here

#UNTESTED CODE TO DEMONSTRATE A WAY AROUND use strict; use warnings; use Bio::SeqIO; my $in = Bio::SeqIO->new(-file =>"FastA.fa", format = "FASTA"); while(my $seq = $in->next_seq){ my $sequence = $seq->seq; #Do something with $sequence my $ORFpattern = "foo"; if($sequence =~ /$ORFpattern/){ #report or anything } }

"Hi Pearl Monks"

we are Perl Monks and not Pearl Monks...

Welcome to the Monastery and good luck....

David R. Gergen said "We know that second terms have historically been marred by hubris and by scandal." and I am a two y.o. monk today :D, June,12th, 2011...

In reply to Re: How do I extract ORFs from a fasta file into a new fasta file by biohisham
in thread How do I extract ORFs from a fasta file into a new fasta file by Wasp_Guy

Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post; it's "PerlMonks-approved HTML":

  • Are you posting in the right place? Check out Where do I post X? to know for sure.
  • Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
    <code> <a> <b> <big> <blockquote> <br /> <dd> <dl> <dt> <em> <font> <h1> <h2> <h3> <h4> <h5> <h6> <hr /> <i> <li> <nbsp> <ol> <p> <small> <strike> <strong> <sub> <sup> <table> <td> <th> <tr> <tt> <u> <ul>
  • Snippets of code should be wrapped in <code> tags not <pre> tags. In fact, <pre> tags should generally be avoided. If they must be used, extreme care should be taken to ensure that their contents do not have long lines (<70 chars), in order to prevent horizontal scrolling (and possible janitor intervention).
  • Want more info? How to link or or How to display code and escape characters are good places to start.
Log In?

What's my password?
Create A New User
Domain Nodelet?
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others browsing the Monastery: (3)
As of 2021-12-06 07:47 GMT
Find Nodes?
    Voting Booth?
    R or B?

    Results (32 votes). Check out past polls.