Monthly Archives: May 2016

ReadMapper: a program that maps sequencing reads to a nucleotide backbone, and creates a publication quality figure.

Introduction

ReadMapper is based on R code that was written to analyze our inhouse data.  The code has been modified to make it user friendly and allow other people, who are not as versed in the R language, to be able to map sequencing reads to a nucleotide backbone, and create publication quality figures.  In our project we were mapping translatome data to the genome of the phage Lambda (this is what the included sample data is), as well as the currently annotated gene open reading frames.  Which looks just like the figure below.

Continue reading

crAssphage poops up again!

Yet again, analysis of a metagenomic sample shows that crAssphage is the most abundant phage anywhere. It also shows what a dis-service NCBI did to science by deleting the crAssphage record. We used meta-spades to reconstruct the entire crAssphage genome from someone else’s data set, but in their paper, the largest contig was ~3 kb. This analysis suggests that crAssphage is present in ulcerative colitis samples but the abundance goes down after treatment!

Continue reading