Question: (Closed) Remove sequence from fasta file by samtools
0
gravatar for Janey
18 days ago by
Janey10
USA
Janey10 wrote:

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 • 119 views
ADD COMMENTlink modified 18 days ago by Ram15k • written 18 days ago by Janey10

searching on this site would give you many solutions.

ADD REPLYlink written 18 days ago by shenwei3563.8k

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

Does no one have any easy solution.

ADD REPLYlink written 18 days ago by Janey10
seqkit grep -v -f ID-list.txt seqs.fa

ID-list.txt: one ID per line.

ADD REPLYlink written 18 days ago by shenwei3563.8k

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!

ADD REPLYlink written 18 days ago by Ram15k
1

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

ADD REPLYlink modified 18 days ago • written 18 days ago by finswimmer2.8k
Please log in to add an answer.
The thread is closed. No new answers may be added.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 961 users visited in the last hour