How to remove aligned sequnces?
0
0
Entering edit mode
4.8 years ago

Hello, I have been working on a metagenomic data, and now trying to remove ribosomal protein using the blastp command in BLAST+, I want to retain the sequences which doesn't map to the reference. Can anyone tell me how can I do that so that all the mapped sequences Bacterial gets removed and I am left with unmapped sequences which are possible Viral sequences.

The command I used for blastp is;

blastp -db cog -query ../../../Prodigal/F_V_0_29_R.faa -out ../../../blast/output_1 -evalue 1e-10

Thanks in advance, Vishal

alignment next-gen sequencing • 675 views
ADD COMMENT
0
Entering edit mode

Parse entries that don't have hits. You can then fish those sequences out of your file. Here are some options for that: Extract Sequence From Fasta File Using Ids From a separate txt File in linux

ADD REPLY
0
Entering edit mode

Thank You so much for the help.

ADD REPLY

Login before adding your answer.

Traffic: 1813 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6