How to extract sub fasta sequences from a multiline fasta file ?

Question

[sunnykevin97]

Hi

Multiline fasta 
>data1
AADFGHJKASDFG
ASDFGHJKLFGHJK
ZXCVBNMMMMM
>data2
POIUYTREWPOIU
LKJHGFDSSALKJH
MNBVCXZLKJHGF
>data3
QWERTYUIOQWE
ERTYUIRTYUITYUI
ASDFGHJSDFGHS

Interested in extracting only >data2 sequence into text file.

>data2
POIUYTREWPOIU
LKJHGFDSSALKJH
MNBVCXZLKJHGF

Using file option that I'd great. I had a 1000's of ID's in a text file. I tried using grep commands able to extract only 1 line.

 >data2
    POIUYTREWPOIU

Help!

Answer 1

[GenoMax]

how to filter a fasta file with seqs ID

[sunnykevin97]

Its not working, generates an empty file.

faSomeRecords fsd.fa mfsd.txt out.fa

[GenoMax]

You need to list plain fasta header names one per line in mfsd.txt. There can be no > in listing.

[sunnykevin97]

Other than faSomeRecords is their any tools that can do the same ? Actually it doesn't have append(>>) option. That doesn't solve my problem yet. Help!

[GenoMax]

What do you mean it does not have append option? All selected sequences should be in out.fa. If you are doing this for different datasets, just cat output files together. cat out1.fa out2.fa > final.fa.

You need to improvise (especially if you are not a programmer yourself) when needed. You will rarely find tools that do exactly what you need since you may be the only person needing a tool that works like that.

[sunnykevin97]

I tried the same, it worked.

Thanks