SEQanswers

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 586 11-07-2017 10:06 AM
Please help my BBMap cannot remove Illumina adapter TofuKaj Bioinformatics 4 04-28-2015 09:53 AM
BBMap Error Phage Hunter Bioinformatics 5 01-14-2015 05:34 AM
Introducing BBMap, a new short-read aligner for DNA and RNA Brian Bushnell Bioinformatics 24 07-07-2014 10:37 AM
BitSeq narges Bioinformatics 1 11-19-2012 02:53 PM

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

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

Hello,

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:

Code:
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:
bbmap.sh 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?

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

Join Date: Jan 2014
Posts: 2,695
Default

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:

"bbmap.sh (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
Reply

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 12:28 PM.


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