Dear all Can someone tell me how to retrieve multi-line fasta file from for instance below mentioned data using a list of locus_tag like locus_tag=ABUW_RS00010, locus_tag=ABUW_RS00015. I have 400 entries of this.
>lcl|NZ_CP008706.1_cds_WP_001237350.1_2 [locus_tag=ABUW_RS00010] [protein=DNA polymerase III subunit beta] [protein_id=WP_001237350.1] [location=1590..2738] [gbkey=CDS]
GTGCGTTTGAAAATCGCTAAAGAAAGTTTACTCAATGTTTTATCGCATGTGGTCGGTGCGGTAGAACGTC
GTCATACCTTGAATATTCTTTCAAACGTCAAAATTCAGACTAATGCTCAAGCGCTCACTATTACGGGTTC
>lcl|NZ_CP008706.1_cds_WP_000550807.1_3 [locus_tag=ABUW_RS00015] [protein=DNA replication/repair protein RecF] [protein_id=WP_000550807.1] [location=2753..3835] [gbkey=CDS]
ATGCATCTTACGCGCTTAAATATTGAACGTGTGCGTAATTTAAAAACGGTTGCACTCCATGGGTTGCAAC
CGTTTAATGTCTTTTATGGCGCAAACGGTTCGGGTAAAACATCAATTTTAGAAGCGATCCATTTGCTTGC
AACAGGTCGCTCATTCCGTACACATATACCCAAAAATTATATTCAATATGAAGCTGACGACGCGATTGTA
TTTGCTCAGTCAGCAACTGAAAAAATTGGTATGCAAAAACTCGCTTCTGGTGAGCAACTCATGAAAGTGA
•
written
by
sharmatina189059 •
2 answers
Hello cpad0112,
I understood the question in that way, that the OP wants to fetch a limited list of sequences out of one fasta file containing many sequences, based on the locus_tag. But maybe I'm wrong here.
If I'm right, than seqkit could help (again):
$ seqkit grep -n -r -f tags.txt input.fasta > output.fasta
Where tags.txt contains the locus_tags.
fin swimmer
Thankyou so much!!! It works. This is exactly I want.
I have moved the comment of finswimmer to an answer since it solves your question.
If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Hello Finswimmer, can your answer above apply to extracting the sequences from the genbank file using the locus tags?
Hello sharmatina189059,
for this type of operations you can use SEDA (http://www.sing-group.org/seda/).
Specifically, the "Pattern filtering" operation (http://www.sing-group.org/seda/manual/operations.html#pattern-filtering) allows to do what you are asking for. You only need to load your FASTA file(s) (or take the sequences from the clipboard) and then apply your text pattern to the sequence headers.
Please, do not hesitate to contact us in case you need further assistant.
Regards,
Hugo.
Log in to answer this question.
-
Renaming fasta files with their headers
written by sebabiokr •
-
Compare two columns from two files and print the content of both the files
written by Utkarsha •
-
unexpected output from efetch
written by andres.firrincieli
-
Error when using featurecounts
written by m.radz •
-
How to extract two genomic location numbers within the following fasta header?
written by mrj •
-
Protein fasta file header shorten
written by wang-yanfang •
-
Filtering bacterial genes by annotations and finding common genes between genomes
written by dch014 •
-
Filtering sequences with multiple headers by length
written by dod •
-
how to change fasta headers?
written by sharmatina189059 •
-
Get taxonomy into a CDS from genomic custom database
written by Shred
Hello sharmatina189059,
could you please reformat your post using the code button (the one with 101010) so that one can see how your input file realy looks like.
Does it look like this?
fin swimmer
retrieve from where?
I added markup to your post for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below: