Extract reads from specific taxa in blob tools?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Phyloseq returning black bars for plot_bar
written by Jacob 0Hi all, I am new to phyloseq and trying to plot a dataset with the plot_bar function. I read in my taxa, otu, and metadata …
-
Extract cis and trans mutations from bam file
written by RD 3Hi 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 15Hello 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 0I would like to make a plot similar to the one here in R. I have a phyloseq object that contains the abundance of all …
-
Collapse multifasta file by specific chromosome names
written by YOUSEUFS 3I have a multicast file with unique identifiers ('SUBJECT.1', 'SUBJECT.2' etc) like this: >SUBJECT.1.1:1203-2742(+) AAATTT >SUBJECT.1:354-700(+) CCCGGG >SUBJECT.2:789-2000(+) GGGCCC >SUBJECT.2:2012-2742(+) TTTAAA how would I extract …
-
How to batch download multiple mitochondrial FASTA & GFF files?
written by YOUSEUFS 3Downloading 1 by 1 is taking an eternity. I have a refseq_mito.txt file with 115 RefSeq accession numbers for mitochondrial genomes (examples below) refseq_mito.txt contents …
-
AFD plot for each VCF file
written by anasjamshed 14I 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 5Hi, 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 2I 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 11Hello 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 …
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