Beefy Boxes and Bandwidth Generously Provided by pair Networks
Think about Loose Coupling
 
PerlMonks  

Re: unique sequences

by Dallaylaen (Chaplain)
on Dec 11, 2017 at 01:42 UTC ( [id://1205274]=note: print w/replies, xml ) Need Help??


in reply to unique sequences

It's hard to tell from the source what the input looks like, and what should be done...

Do I understand it right that
1) Input has [ACGT] sequences + some comment lines starting with >?
2) You need to find strings of 12 nucleotides ending in a GG that occur in the input only once?

It would be nice if you posted a short input sample (although I guess it's may not generate the needed discrepancy).

Replies are listed 'Best First'.
Re^2: unique sequences
by erix (Prior) on Dec 11, 2017 at 08:58 UTC

    These kind of dmel-all-chromosome files can be found here:

    ftp://ftp.flybase.net/releases/current/dmel_r6.18/fasta/

    dmel-all-chromosome-r6.18.fasta.gz

    That's a few releases later but I suppose that doesn't matter for the processing solution. (The older releases will be somewhere there as well.)

[reply]
      That's a few releases later...

      So basically species have releases, and are versioned these days. I for one would love to git revert a mammoth!

[reply]

In Section Seekers of Perl Wisdom

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?

www.com | www.net | www.org
Node Status?
node history
Node Type: note [id://1205274]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this? • Last hour • Other CB clients
Other Users?
Others taking refuge in the Monastery: (6)
As of 2024-04-19 08:59 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found