Beefy Boxes and Bandwidth Generously Provided by pair Networks
laziness, impatience, and hubris
 
PerlMonks  

comment on

( #3333=superdoc: print w/replies, xml ) Need Help??
So you basically have some coordinates and want to check if anyone of the alleles matches the location given by the coordinates right ? But your coordinates don't tell about which chromosome these alleles are expected but what you can probably do is download the chromosomes in say FastA format and then use Bio::SeqIO using the subseq function.

here is the documentation to the module Bio::SeqIO

The subseq function requires a beginning and an end to what it is going to extract so probably you want to reformat your coordinates as follows

$seq->subseq(coordinate-1, coordinate+1)

to extract three bases with the middle one being the actual nucleotide you want to compare your alleles against.

Note that BioPerl can be a bit slow and therefore consider using something like sfetch or esl-sfetch which will again require providing beginning and end positions to what it is to extract


David R. Gergen said "We know that second terms have historically been marred by hubris and by scandal." and I am a four y.o. monk...

In reply to Re: bioperl module to extract specific nucleotides given chromosome and exact location of that nucleotide by biohisham
in thread bioperl module to extract specific nucleotides given chromosome and exact location of that nucleotide by xxArwaxx

Title:
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?
Username:
Password:

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

How do I use this? | Other CB clients
Other Users?
Others drinking their drinks and smoking their pipes about the Monastery: (5)
As of 2021-12-01 21:11 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    R or B?



    Results (14 votes). Check out past polls.

    Notices?