How to split fasta file by ID, making sure to keep two entrys in each file
Hi, I have a fasta file which looks like:
>reference
AGCT
>entry1
AGCT
>entry2
AGCT
I want to split the fasta file by ID, but also include the reference in each file.
So far I am using seqkit split --by-id file.fasta and then cat reference.fasta file1.fasta > results.fasta
However this is quite slow, is there a way of doing this in one command?
• 1,168 views
•
link
updated
by
Istvan Albert
•
written
by
SaltedPork
•
written
by
SaltedPork
1 answer
If you have already split your files then it is just a question of looping through the list, you can loop in bash or more elegantly with parallel:
ls -1 *.fa | parallel 'cat reference.fa {}.fa > {}.merged.fasta'
• 0 views
•
link
written
by
Istvan Albert
Log in to answer this question.
More posts like this
-
Remove duplicates in multifasta, where entries are paired
written by SaltedPork
-
Split multi-fasta file and keep structure
written by genomes_and_MGEs •
-
Having trouble with Seqkit Replace, trying to rename fasta headers from file
written by SaltedPork
-
Use Read ID File to parse out reads from a Fasta File
written by kalanir •
-
Problem with FASTA filtering
written by Janey •
-
How to split fasta by '>' into a file each containing one sequence, and have the name of that file …
written by SaltedPork
-
Possible bug in seqkit?
written by lakhujanivijay
-
merge large amount of fastq files into a single one
written by catherine
-
compare and trim fasta files by header ID
written by joie.cannon •
-
Fastq File Split Using Corresponding Id Identifier In Tag File, Not Barcode Sequence
written by Tonyzeng •
