Extract reads from specific taxa in blob tools?
0
0
Entering edit mode
2.3 years ago
YOUSEUFS ▴ 30

My blob tools plot looks like:

blobtools viewer output

I want to extract all reads associated to Escherichia (or any taxa), how would I go about doing this?

genomics metagenomics blobtools • 615 views
ADD COMMENT
1
Entering edit mode

has been a while since I used blobTools, but if I remember correctly you should have a file called blast.out that should look like this:

contig_1    979556  200
contig_2    979556  500
contig_2    979556  1000
contig_2    979556  500
contig_2    979556  300
contig_3    979556  10000
contig_4    979556  1000
contig_5    6252    2000
contig_6    232323  2000
contig_6    6252    2000
contig_6    979556  2000
contig_6    232323  2000
contig_7    6252    2000
contig_8    6252    2000
contig_8    979556  2000
contig_9    6252    200

The second column is the NCBI taxonomy ID. Use this information to extract the reads mapping to a specific set of contigs by using the blobtools module bamfilter

ADD REPLY

Login before adding your answer.

Traffic: 1721 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