thank you, but this is whole of IDs, for example if I want only 100 IDs from which what should I do?
getting sequence for a list of miRNAs
hi,
I have a list of miRNAs IDs from different organisms how I can take sequence of mature miRNAs for all my IDs?
for example
lgi-mir-34
mdo-mir-193
mdo-mir-21
mdo-mir-302d
mml-mir-18
• 2,673 views
•
link
2 answers
From miRBase using this file.
• 1 views
•
link
• 1 views
•
link
grep -A 1 your_ID them out of the miRBase file using a loop OR use filterbyname.sh from BBMap OR use samtools faidx as described here: A: How To Extract A Sequence From A Big (6Gb) Multifasta File ?.
• 1 views
•
link
another way was mirdeep2
fastaselect.pl reads.fa reads_select.ids > reads_select.fa
This script only prints out the fasta entries that match an id in the id file.
• 0 views
•
link
Log in to answer this question.