Extract fasta sequnces not matching a list
Hey guys,
I know it's possible to use seqtk subseq to extract fasta files matching a list. Now, is there a way to extract those fasta files that don't match the exact fasta headers on the list? Thanks!
• 3,677 views
•
link
updated
by
GenoMax
•
written
by
genomes_and_MGEs •
•
written
by
genomes_and_MGEs •
1 answer
You can use faSomeRecords from Jim Kent's utils for that.
faSomeRecords - Extract multiple fa records
usage:
faSomeRecords in.fa listFile out.fa
options:
-exclude - output sequences not in the list file.
Log in to answer this question.
More posts like this
-
Parsing genbank or gff files based on specific nucleotide coordinates
written by genomes_and_MGEs •
-
Remove whitespaces on fasta files, except on fasta-header
written by genomes_and_MGEs •
-
Only copy multifasta files from one directory to another
written by genomes_and_MGEs •
-
Extracting fasta sequence from CDS matching with sequence ID
written by metallicasster •
-
Batch rename protein fasta headers
written by genomes_and_MGEs •
-
List names in column that match a given word in another column
written by genomes_and_MGEs •
-
Add strain name directly on fasta header
written by genomes_and_MGEs •
-
Extract fasta-headers that have the same word
written by genomes_and_MGEs •
-
How to remove a list of reads from fastq file?
written by Tao
-
Problems in use of seqtk subseq
written by donger1103 •
Looks like you already have a solution. In general, this boils down to list manipulation. If you have a list of all sequences and a list of those that you don't want (those to exclude), grepping the two files using a
-voption will give you their difference that can be used in normal fashion.