Tag: readDNAStringSet

How to add a refseq() slot to a MAG based phyloseq object ?

How to add a refseq() slot to a MAG based phyloseq object ? 0 Hello, I was asked to build a phyloseq object from some metagenome assembled genomes that I built. By far I have created a tax_table() slot from the classification of the MAGs using GTDB-Tk, a out_table() that…

Continue Reading How to add a refseq() slot to a MAG based phyloseq object ?

Rarefaction curve of microbiome study

I don’t know why my rarefaction curve for most samples does not reach saturation. Please how can I improve te curve? The curve does not look like most rarefaction curves I see in publications, Please help and I am relatively new to microbiome study. Below is the code I used:…

Continue Reading Rarefaction curve of microbiome study

What is Bioconductor in R ?

Bioconductor is an open-source and open-development software project. It provides tools, packages, and resources for the analysis and comprehension of genomic data. Focuses on the statistical analysis and interpretation of high-throughput biological data. These packages include preprocessing, quality control, normalization, differential expression analysis, pathway analysis, genomic annotation, visualization, and machine…

Continue Reading What is Bioconductor in R ?

Query regarding phyloseq object construct with QIIME output

Query regarding phyloseq object construct with QIIME output 0 @6d5973d2 Last seen 15 hours ago India I used qiime pipeline to analyze 16s amplicon sequencing data now I want to take those outputs into r studio with the help of phyloseq package and want to create phyloseq object that’s how…

Continue Reading Query regarding phyloseq object construct with QIIME output

Introduction to Computational Evolutionary Biology

R is a very flexible programming language, and it allows developers to create their own data structures (called classes) for their packages. Over the years, some packages have become so popular that the classes they use to store data are now used the “standard” representations for particular types of data….

Continue Reading Introduction to Computational Evolutionary Biology

bioconductor – Getting a consensus sequence per species from a fasta with sequences from many species in R

In R, I have a fasta (ex: Andrenidae.FASTA) with a few hundred nucleotide sequences from a dozen species. I want to get 1 consensus sequence for each species. Each sequence is made up of A/C/T/G/N, with N representing unknown nucleotides, and the sequence names are species names. All sequences are…

Continue Reading bioconductor – Getting a consensus sequence per species from a fasta with sequences from many species in R

r – Displaying BIC and AICs from MSA nucleotide multisequence fasta file

I am trying to replicate MEGA (models/find best DNA protein models) using R. After reading R AICs and BIC documentation I can’t understand how I can implement it. How can I implement AICs and BICs without having to especify the number of sequences in the fasta file (in case that…

Continue Reading r – Displaying BIC and AICs from MSA nucleotide multisequence fasta file

error reading fastq-files with readDNAStringset

I am trying to read a fastq-file with readDNAStringSet and having quite some trouble doing so. I need the names, aswell as the quality-scores. Right now I am using: readDNAStringSet(myFastqFile, format=”fastq”, use.names= TRUE, with.qualities = TRUE) But here i get the Error: “@” expected at beginning of line 1 I…

Continue Reading error reading fastq-files with readDNAStringset

Trimming DNAStringSet

Trimming DNAStringSet 1 Hello, I am currently dealing with the problem of reading in a Fastq-File with “readDNAStringset”, trimming the Sequences and then writing them in to a new fastq-file. The reading of the fastq-file with “readDNAStringSet” is working just fine. I am then trying to trim a fixed length…

Continue Reading Trimming DNAStringSet