hi,
I have a long list of different miRNAs like below
hsa-miR-7641
bta-miR-2904
hhi-miR-7641
hsa-miR-4454
bta-miR-2478
efu-let-7c
mmu-miR-6240
hsa-miR-7704
is there any way to avoid searching one by one in mirbase and copy pasting their mature sequences and taking their sequences at the same time?
It looks like you can download sequence files: http://www.mirbase.org/ftp.shtml
"Fasta format sequences of all miRNA hairpins" & "Fasta format sequences of all mature miRNA sequences"
When you have those files, a simple grep -f <youridentifiers.txt> <mirbase.fasta> would do the trick
Hi ,
You can use faSomeRecords to extract multiple fasta records
Usage:
faSomeRecords mature.fasta your-list-file result.fa
in the grep -f/-F be careful of space and blank lines. you can check them with Notepad++ in windows.
Take care.
That is not cool.
Whole idea is to use "upvotes" judiciously so you can acknowledge posts that are offering significant help that either solve or go towards solving your problem. Since all communication on Biostars stays public, if every exchange is "upvoted" then that defeats the purpose of the voting system.
Please don't do this.
In addition all of your posts (Angel) are being automatically upvoted by
Who is also from Iran. That's why you have 2K votes when in fact you don't help other people.
I am very sympathetic for people that live in countries where there is a ... let's say "lack of democracy" but please do understand you are helping no one - and you are not serving any cause by gaming the system and pointlessly voting each other up.
Please remove the phony votes otherwise all these votes will be removed and your scores will be set to zero.
Dear @Ram, @Istevan and other nice people in Biostars; of course I do not feel any necessity to answer but as you "tagged me in" I must say that there is no "reciprocity" between me and any other people in here unless their response or answer is valuable (from scientific point or even teaching us some kind of correct appreciation or offering some helpful scripts).
In addition, I do not know how to automatically upvote some one else profile for 1K in a minute!
The scores in this (and any other) websites are not important for me (of course solving the people problems and my problems is valuable for me).
Dear @Angel, please do not act somehow that have negative feedback for my personality as a friend and do not take any action that the others call people of your country "poor in culture of using online voting system due to lack of democracy" !
I'm not sure under which comment my reaction is the most appropriate, but it's mainly targeted to Farbod and Angel.
Don't misunderstand the reactions above, we value every contribution to this website and are happy to have you here so frequently. But since all of us are volunteers, and votes are the only type of "recognition" that people can get for the time and effort that is spend here, it's a rather important system. And just like with any other thing, if you don't spend these votes "wisely", they lose their value. It's understandable that you want to upvote people who have been helpful to you in the past, but just do that for posts in which they actually are helpful. That way, other people in the community can also see which posts have a higher value than others.
"If every post is upvoted - no post is"
believe I did not that to disturb any one, I was only joking with farbod. I am sorry for doing like a childish act while I am too busy but I spent about an hour to joke might be that is why I did not learn data analysis after two years :( :( :( actually Iranian people is not all alike me :( :( they are good and affectionate people but exceptionally I usually do unreasonable things.
You should not have to do this, Istvan. There are rules and etiquettes, and our fear that we'd become SO is tying our hands from any kind of disciplining.
I've been silently monitoring their activity over the past few weeks and not intervening, and now you can see the thing Goutham Atla and I were fearing would happen.
But yes, there are a few holes we could plug - the trophy on answering own questions - it could be made less frequent (once in six months or so).
you know Ram, actually poor Farbod did not know anything about my ridiculous action and jivarajivaraj user is my previous profile I created two years ago when I was starting to learn data analysis. I did not know my action would irritate admins. now what is suffering me is whether biostars in future continue to help me as already by this illogical manner :( :( :(
actually I am not going to make intervention here, you, Devon, Goutham, genomax2 , WouterDeCoster and many others saved me many times from puzzling questions. for example newly @Persistent LABS solved my problem and last year @Antonio R. Franco helped me to run DEseq2 what even I can't do that myself now. today I was roaming in biostars, suddenly I thought what happen if my scores go higher like a moderator then I did so and I also made a joke with Farbod posts but he did not know anything.
now please excuse me and treat me as already and I promise be a good user might be have progress in my learning.
Sorry about the "lack of democracy" quip. I was in a hurry - all I wanted to say that I believe your government does not treat people equally - and that has the tendency to have people look out for one another more and help each other however they feel is necessary.
All positive sentiments - not trying to be judgmental or anything - one of the problems with the internet is the lack of metacommunication - it is easy to state something in writing that sounds very different than intended.
Your culture is not being targeted, not is your being active on the site. We are raising a valid point on your misuse of the voting system. We welcome everyone and are happy to have you as part of the community, until you start abusing the system. It demeans the amount of effort Istvan and the mods put into this.
Dear Ram, if you are with me,
I did not "misuse" or "abusing" your system.
The scores here is not important for me (not at all).
I have repeat it several times here but it seems that I was not so clear, So for the last time: I did not try to maximize mine or others scores! (I have upvoted "genomax2" or "Wouterdecouster" or "Denis" many times more than "Angel" - haven't you realized that in your silently monitoring our activity over the past few weeks)
And maybe I am not so expert in the field of bioinformatic, as it was not my field of interest in university, but I have learned in my life that it is not a good manner to anathematize all the people of a country and speaking about the democracy situation in their country !
I think we all made our point and can reach a mutual understanding of the voting system. We should avoid making this personal or discuss about the role of cultural elements in this. Perhaps it's time to end this discussion and move on to more pleasant bioinformatics.
No one's targeting your country, Farbod - I don't understand where you get that from. Who spoke of your nationality or ethnicity?
And I reported what I observed - threads between Angel and yourself with no content but one upvote each. As a mod, it is my job to monitor and check for odd usage.
I think it is fine if people respond to each other preferentially. Or to seek out each other if they wish so and contribute that way.
There is nothing wrong with that - this is a social site after all. As long as posts are on topic everyone is welcome to make use of it the way want to.
The only issue with misusing the voting system - it undermines its value.
Dear Angel Hi,
I really appreciate your kindness and generosity (to tell you the truth, it was somehow bizarre for me),
but I use Biostar website as a place to learn and enjoying one can find in solving other person problems, so as maybe you have realized, my SCORES was not so high as it was not one of my goals, gaining them. If you remember, I have suggest to run your miRNA data analysis on my server and be sure that it was not for maximizing my online scores and just for helping one of my compatriot.
I really do not want to be rude specially with a Fereshteh from Iran, but . . .
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
version 2.3.6
thank you,
grep -f <youridentifiers.txt> <mirbase.fasta> > result.txt
You need to redirect the output to a file, e.g.
grep -f youridentifiers.txt mirbase.fasta > mysequences.fasta. Without that redirection, grep just writes to stdout, which is the terminal.You write "As a txt file", what exactly do you have in mind as output? With standard linux tools (grep, sed, cut, awk) you can probably do all processing you would like.
Oh good, I see you already found out how to get it in a file ;)
yes exceptionally sometimes grey cells in my brain work :)
sorry I inspect the result file that contained all of mature miRNAs while my identifiers.txt contained about 80 miRNAs
:(
then what to do?
Is there by chance a newline (empty line) at the bottom of your identifiers.txt file? Additionally, might be better to go for
grep -w -f identifiers.txt mirbase.fastato only get full 'word' matches and not partial 'pattern' matches.thank youuuuuu
exactly the problem was some empty lines at the end of my id file
Yeah, and then everything matches to the newlines :p (I have this problem about 95% of the times that I do something like this)
I tried these codes all give me whole of fasta file
grep -A1 -w -f id.txt seqFile.fasta > output.fasta
grep -Fwf miRNA.txt -A 1 mature.fasta | grep -v '^--$' > out.fasta
grep -Fwf < ( sort -u miRNA.txt ) -A 1 mature.fasta | grep -v '^--$'
Could you try
grep -A 1 -Fwf <(tr -d '\n' miRNA.txt) mature.fasta > result.fasta? I see you do grep -v at the end to remove newlines, but you should remove newlines from the identifier file (although I have my doubt whether^--$would work).thank you
grep -A 1 -Fwf <(tr -d '\n' miRNA.txt) mature.fasta > result.fasta
told missing name for redirect
sorry,
supposing a file containing
dataset: 1
target: EAL01688
length: 507
miRNA : hsa-miR-26b-3p
length: 22
mfe: -18.2 kcal/mol
p-value: 1.000000e+00
position 190
target 5' A UACCU AAA A A 3'
miRNA 3' C UUCAUU 5'
dataset: 1
target: EAL01691
length: 2484
miRNA : hsa-miR-26b-3p
length: 22
mfe: -23.2 kcal/mol
p-value: 1.000000e+00
position 913
target 5' G CA GGA AAAGA U 3'
miRNA 3' A C 5
how to extract only unique target for example EAL01688
Which output do you want? I'm not sure I understand.
thank you
actually I was searching for targets of my miRNAs by rnahybrid as I presented some parts of the results above. now I want to extract only target from this output specially unique target. for example if EAL01691 repeated some times I only need that once.
If I understand it correctly, something like this would work:
grep "target:" yourinput.txt | sed 's/"target: //' | sort -u