Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
BBMap (aligner for DNA/RNAseq) is now open-source and available for download. Brian Bushnell Bioinformatics 678 12-12-2019 08:36 AM
Please help my BBMap cannot remove Illumina adapter TofuKaj Bioinformatics 4 04-28-2015 08:53 AM
BBMap Error Phage Hunter Bioinformatics 5 01-14-2015 04:34 AM
Introducing BBMap, a new short-read aligner for DNA and RNA Brian Bushnell Bioinformatics 24 07-07-2014 09:37 AM
BitSeq narges Bioinformatics 1 11-19-2012 01:53 PM

Thread Tools
Old 04-29-2015, 09:48 PM   #1
Senior Member
Location: Vienna

Join Date: Mar 2010
Posts: 107
Default BBMap for BitSeq


I want use BBmap for estimating transcript expression level from RNA-seq reads. Mapped will be against transcripts, not the genome.

The Manual for BitSeq uses Bowtie with following conditions:

bowtie -q -v 3 -3 0 -p 4 -a -m 100 --sam
3 mismatches allowed
no 3' trimming (and no 5' trimming) - both default
report reads only if < 100 possible mapping positions

how should the BBmap parameters look like?

Code: ambiguous=all maxsites2=100 ( secondary=TRUE sssr=0.95 maxsites=100 )
is there a parameter for allowed mismatches per read?
will "secondary=TRUE sssr=0.95 maxsites=100" be a good idea for mapping reads to the human transcriptome. sssr=0.95 means approximate how many mismatches per 75 bases reads ~ 2?

should the reads be trimmed and quality filtered before, and if yes with which parameters?

dietmar13 is offline   Reply With Quote
Old 04-30-2015, 08:40 AM   #2
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

I imagine the main reason for specifying 3 mismatches with Bowtie is that 3 is the maximum it allows.

BBMap normally adjusts sensitivity with the mind/minratio parameter, though you can optionally set "subfilter=3" to ban alignments with more than 3 substitutions. I don't see how that would be beneficial to transcriptome mapping, though.

I would suggest:

" (file parameters) ambig=all maxsites=100 maxindel=100"

...and if you really want, add "subfilter=3". For transcriptome mapping there's no reason to allow the default maxindel=16000. And there is not much need to apply any quality-trimming or filtering unless you add the "subfilter" flag and have low-quality data, though you can do that if you want with "qtrim=rl trimq=10" to trim both ends to Q10. I do, however, always recommend adapter-trimming, particularly when requiring high-identity alignments.
Brian Bushnell is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:42 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO