This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract sequence ID from Fasta for a given sequence

I have a fasta file myfasta.fasta like this:

>aat.2.2344.a
ATTGCCGGTTTAATATTA
>aat.2.d2344.acc
ATTGCCGGTTTAATAAA
>aat.2.2bb344.a
ATTGCCGGTTTAATAGGAGAGAATT
>aat.2.2ccc344.a
ATTGCCGGTTTAATAGGGAG
>aat.2.2344.acc
ATTGCCGGTTTAATAAA

I also have a text file my.txt which contains the sequence that matches the sequence in fasta file above:

ATTGCCGGTTTAATAAA

Based on this sequence, I want to extract all matched IDs for this sequence. Can someone please help me with this? Thanks!

The result I want is:

>aat.2.2344.acc
>aat.2.d2344.acc
fasta
• 2,187 views • link
updated by michael.ante michael.ante • written by MAPK MAPK

Are the sequences all one line? If so you can just use grep -B 1 ...

• 1 views • link
written by Joe Joe

Yes they are 50 bps reads.

• 1 views • link
written by MAPK MAPK

Dear MAPK, if you usually work with FASTA files you may find SEDA (http://www.sing-group.org/seda/) an useful tool. It has a great variety of operations to manipulate, filter, and transform FASTA files (check out the manual to see all of them: https://www.sing-group.org/seda/manual/index.html). It also allows you to explore a set of FASTA files and extract only the information you need, such as the sequence identifiers (see https://www.sing-group.org/seda/manual/graphical-user-interface.html#the-input-area).

With best regards, Hugo.

• 1 views • link
written by Hugo Hugo

2 answers

This works:

grep --no-group-separator -B 1 -F -f my.txt 344194.fasta | grep -v -f my.txt

Only downside is reading the my.txt file twice. There are other non-grep approaches that could avoid this but this is simple.

• 1 views • link
written by Joe Joe

Hi MAPK,

I find fasgrep from the FAST suit quite handy:

fasgrep -s ATTGCCGGTTTAATAAA myfasta.fasta

It reports the full entry, thus you can just grep for ^\>.

[EDIT] since you have multiple sequences (as I have read right now) you can provide these as regex: "ATTGCCGGTTTAATAAA|CCCCGCGC|ATATATATA"

Cheers,

Michael

• 1 views • link
written by michael.ante michael.ante

Log in to answer this question.

More posts like this