using BBmap to map paire end reads to de-novo assembled transcriptome · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

using BBmap to map paire end reads to de-novo assembled transcriptome

Dear all

I received some RNA seq reads containing Lexogen primers. My reads do not have any refrence genome. Because of existing primers in the reads, mapping my reads by TopHat to assembled transcriptome will reduce the mapping rate significantly. Regarding having no reference genome, Can I map my reads by BBmap to assembled transcriptome provided by Trinity? Dose BBmap need to reference genome for mapping?

rna-seq

• 2,342 views • link

updated 8.7 years ago by GenoMax

• written 8.7 years ago by rahmati.razieh83 •

1 answer

If you are going to align then you will need a reference no matter what genome you are working with. You can use trinity assembled transcriptome as a reference with BBMap. There are other aligners like GMAP that may be useful for this as well.

• 0 views • link

written 8.7 years ago by GenoMax

Thank you so much for your comment

• 0 views • link

written 8.7 years ago by rahmati.razieh83 •

Log in to answer this question.

More posts like this

BBMap not mapping reads to sequences known to exist in reference
written 5.4 years ago by robert.murphy •

I am trying to use BBMap to map some degenerate primers to an assembled reference genome. The primers are for the adenylation domain of NRPS …
Mapping both PE reads and Single-end reads into transcriptome
written 7.8 years ago by pbigbig

Hi everybody, I built several de novo transcriptomes from non-model organism without genome reference by Trinity. Sequencing libraries were prepared with stranded method (dUTP), therefore …
mapping and quantification of paired end reads using BBmap
written 8.7 years ago by rahmati.razieh83 •

Hi everyone I have some paired end RNA -seq reads without reference genome. I assembled them by using trinity. because my reads containing specific Lexogen …
trimming PE RNA Reads
written 8.8 years ago by rahmati.razieh83 •

dear all; I have problem in the trimming my RNA seq reads. The problem is that although the length of all reads are the same, …
GMAP error - Problem Sequence: TRINITY_DN422143_c0_g1_i1 (301 bp)
written 8.9 years ago by TrentGenomics •

Hello, I'm using GMAP-2017-09-11 as part of the RECOT (read coordinate transformer; http://sesejun.github.io/recot/ ; https://scfbm.biomedcentral.com/articles/10.1186/1751-0473-8-6) pipeline to transform the coordinates of a transcriptome assembly to …
low mapping rate, finding possible source of contamination
written 8.9 years ago by vaslanzadeh •

Hello, I have received small RNA-seq data from a pathogenic bacteria for bioinformatic analysis. After trimming the adaptors with trim_galore, fastq file was mapped (almost …
low number mapped reads in ribosome profiling
written 10.5 years ago by laura.fancello •

Dear all, I performed both genome and transcriptome mapping on **ribosome profiling data **and I obtain** low percentages of mapped reads**: 14-20% for genome and …
troubleshooting bowtie2 low mapping rates
written 10.5 years ago by zizigolu

Hi, I am mapping my illumina reads on the transcriptom by bowtie2 but the rate of alignment is too LOW(0.9). When I tried my reads …
Align long nanopore minion cDNA reads
written 11.0 years ago by Nicolas Rosewick

Hi, I received some nanopore minion cDNA samples. My idea was to use STARlong to align them against my reference genome. Any other idea / …
How to deal with high percentage of non-specific matching during read mapping on Trinity assembly?
written 11.1 years ago by seta

<p>Hi everybody,</p> <p>Regarding Trinity tool for transcriptome assembly, I noticed that it has high percentage of non-specific matches during mapping reads back to assembled transcriptome. …