Extract reads from specific taxa in blob tools?
My blob tools plot looks like:
I want to extract all reads associated to Escherichia (or any taxa), how would I go about doing this?
genomics
metagenomics
blobtools
• 875 views
•
link
updated
by
Ram
•
written
by
YOUSEUFS •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Is there too little variation to detect meaningful gene changes?
written by Megan •Hi everyone, I wanted to get some feedback on some DESEQ analysis I am currently doing. I have tried to run my pipeline and my …
-
Extract cis and trans mutations from bam file
written by RD •Hi All, I have annotated SNPs from a VCF file, and at position 110294564, there is a primary mutation, let's call it P260R. I want …
-
how to plot distribution of SNPs across a set of genes based on distance between snp and gene
written by rheab1230 •Hello everyone, I have a gene set file. I want to plot a distribution plot for SNP distance across different gene sets. I want to …
-
Plot in R: Presence of bacteria type A vs bacteria type B in different groups
written by flourde •I would like to make a plot similar to the one here in R. I have a phyloseq object that contains the abundance of all …
-
How to manage GO terms for metagenomics
written by young_bioinformaticianHello, I have been trying to retrieve GO terms for metagenomics for a while. I already retrieved all GO IDs, however I can't collect in …
-
AFD plot for each VCF file
written by anasjamshedI have 2 test directories: ![dir][1] and one tsv file that contains information about the population: ![pop][2] [1]: /media/images/453029d7-d7c4-4c3a-8334-e117fd3d [2]: /media/images/00802d58-b802-4145-a81a-82d03c99 I want to write …
-
Help with R code to get a unique list
written by Space_Life •Hi, I have gene ontology data that looks like this (sample data). I manually downloaded the gene ontology data from the Uniprot website. ![Data view][1] …
-
Polygenic Risk Score Plot
written by Muhammad •I have a dataset of about 6000 people, Chromosome 21. I calculated PRS using plink and PRSice. The plot is shown above. I have two …
-
extracting Genus abundance info from Kraken 2 output file
written by arshad1292 •Hello I have a kraken 2 output (k2report.txt file that looks like this: ![enter image description here][1] [1]: /media/images/000a25cb-600d-4b62-a41e-23a5ecef What I want to do is …
-
Visualization of GO results using enrichplot
written by christine.a.persen •Hello I am doing gene set analysis on DNA methylation microarray data (Illumina 450k and EPIC), and use the missMethyl package to do the analysis, …
has been a while since I used blobTools, but if I remember correctly you should have a file called
blast.outthat should look like this: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