This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Remove sequence from fasta file by samtools

Hi

I use the following script to extract sequence from fasta file by samtools. How can I write similar command to remove the sequences by sequence IDs list from fasta file by samtools.

cut -c 2- ID.text | xargs -n 1 samtools faidx in.fasta > out.fasta

Thanks for your help

rna-seq

searching on this site would give you many solutions.

I searched this site but i didn't find suitable solution.

Does no one have any easy solution.

seqkit grep -v -f ID-list.txt seqs.fa

ID-list.txt: one ID per line.

Hello Janey!

We believe that this post does not fit the main topic of this site.

Please search the site - this is one of the most common questions here.

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

It is also a duplicate post from the same OP. Please stay in your original post and give reply to my question and the answers of the user users.

fin swimmer

0 answers

No answers yet.

Log in to answer this question.