How to match fasta header list of name?
Hi,
I have a multi fasta file that contains 453 fasta sequences that looks like the following:
>M.Bce12308ORF4755P GTWWAC
ATGCGTGACCTGATCGAAGAGCCGGGCGGCGGCGCCGCGAGCGAGGCGGAGGCGGTTCAGCCCGCCGCTGCCGTGCCGCGCGCGCTGCCGTCCGGTATCG
>M.Bce1254ORF9725P GTWWAC
ATGCGTGACCTGATCGAAGACCCGGGCGGCGGCGCCGCGAGCGAGGCGGAGGCGGTTCAGCCCGCCGCTGCCGTGCCGCGCGCGCTGCCGTCCGGTATCG
And I have sequence name list that contains 461 name as below:
M.Bce12308ORF4755P
M.Bce122ORF1082P
M.Bce12308ORF4755P
M.Bce1254ORF9725P
May I know how to match the name list to the fasta file, so that I can know which of the sequence from the name list is missing in the fasta file?
Thank you!
Felix
• 4,230 views
•
link
updated
by
shenwei356
•
written
by
fec2 •
•
written
by
fec2 •
2 answers
# IDs in seqs.fa
$ grep '^>' seqs.fa | awk '{print $1}' | sed 's/^>//'
M.Bce12308ORF4755P
M.Bce1254ORF9725P
# IDs not in list.txt
$ grep -w -v -f <(grep '^>' seqs.fa | awk '{print $1}' | sed 's/^>//') list.txt
M.Bce122ORF1082P
• 0 views
•
link
written
by
shenwei356
library(Biostrings)
library(tidyverse)
## get target seq names
target_seq_names <- c("M.Bce12308ORF4755P", "M.Bce122ORF1082P", "M.Bce12308ORF4755P", "M.Bce1254ORF9725P")
## get seq names from fasta file
fasta_seq_names <- Biostrings::readDNAStringSet("input.fasta") %>%
names() %>% ## get names
gsub(pattern = "\\s.*" ,replacement = "" ,x = .) ## clean headers. remove stuff after first space
fasta_seq_names
[1] "M.Bce12308ORF4755P" "M.Bce1254ORF9725P"
## present in both
intersect(fasta_seq_names , target_seq_names)
[1] "M.Bce12308ORF4755P" "M.Bce1254ORF9725P"
• 0 views
•
link
written
by
Chirag Parsania
Log in to answer this question.
More posts like this
-
Replace multiple text with corresponding text
written by fec2 •
-
How to reduce the length of Ns in a fasta sequence
written by hemu.csb •
-
annotate unknown peptides by fasta reference
written by robinycfang •
-
Pairwise alignment of multi-FASTA file sequences
written by aurora •
-
Split multifasta file using awk command
written by fec2 •
-
Blast locally on multiple fasta files with multiple database
written by fec2 •
-
Edit header of multifasta file
written by fec2 •
-
How to concatenate multiple fasta file
written by fec2 •
-
Regarding Whole Genome Visualization
written by namgyalwangstuklama •
-
Glimmer 3 multi-extract input problem
written by thepunkoklos •
Did you try searching the forum at all? This is one of the most widely addressed problems here.
Hi, I have tried searching it, but couldn't found exactly same issue, as what I want is a list of name of the missing fasta file. Do you have any idea what is the key word should I use to search for this issue? Thanks.
Extract reads from fasta file (specific read_names)and make a new fasta file
Extract fasta sequences from a file using a list in another file.
Filtering sequences from a Fasta file using Biopython
Pulling out gene names from a FASTA file with a list of IDs
Extracting specific IDs + sequence from multifasta
How to extract specific genes from a fasta file
Printing specific sequence and ID from combined fasta file using bash commands
Extract Sequence From Fasta File Using Ids From a separate txt File in linux
....
Although I highly doubt you could not find the exact same case (I recall addressing multi-part sequence identifiers and how to deal with them a few years ago), this approach is not helpful. What you need is not something that you can copy-and-paste and "just works" - such solutions are rare and don't teach us anything. You need a "pointer", which is a hint that takes you one step closer to a solution than you are right now. That way, you get to solve the problem yourself while overcoming an obstance that might have taken quite some time to solve on your own.