Alignment of mirna on specific organism data from mirbase
2
0
Entering edit mode
13 months ago
otieno43 ▴ 30

I have micro-RNA sequences that I need to align specifically to Drosophilar melanogaster mature microRNA from mirbase using bowtie. Have realized the fasta sequence of mature miRNA in mirbase contain all the organisms that whose microRNAs have been sequenced and submitted. I was wondering how I can extract/or specifically use D. melanogaster as reference as opposed to using all the mature miRNAs.

mirbase miRNA • 1.5k views
ADD COMMENT
0
Entering edit mode

Do not type in a sentence to use as tags. Use relevant subject matter keywords. I'll try and fix it for you this time but please apply yourself a little more in the future.

ADD REPLY
0
Entering edit mode
13 months ago
shelkmike ★ 1.4k

A simple solution is to download the file https://mirbase.org/download/mature.fa , then open it in a text editor (for example, Notepad++) and manually copy all sequences that have "Drosophila melanogaster" in their titles to another file.

ADD COMMENT
0
Entering edit mode

A simple solution (...) then open it in a text editor (for example, Notepad++) and manually copy all sequences that have "Drosophila melanogaster" in their titles to another file.

no, this is not a simple solution.

ADD REPLY
0
Entering edit mode

There are different definitions of simplicity.

ADD REPLY
0
Entering edit mode

Sure, simple as in easy and simple as in foolish (like in simpleton). Your solution is along the latter lines.

ADD REPLY
0
Entering edit mode
13 months ago

assuming https://mirbase.org/download/mature.fa there is only one line per sequence and there is 'Drosophila melanogaster' in each header:

wget -q -O - "https://mirbase.org/download/mature.fa" |\
grep -F 'Drosophila melanogaster' -A1 --no-group-separator > melanogaster.fa
ADD COMMENT
0
Entering edit mode

Thank you, sorted this out.

ADD REPLY
0
Entering edit mode

Please accept the answer so the question is marked solved on the website. To do that, click on the green check mark on the left side of the answer.

ADD REPLY

Login before adding your answer.

Traffic: 2375 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6