SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics
Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie beginner question... milesgr General 6 03-14-2012 06:21 AM
Simple Tophat question kwatts59 RNA Sequencing 2 06-02-2011 07:31 PM
bowtie-inspect question dicty Bioinformatics 0 12-02-2010 04:45 PM
A simple question on cDNA and protein sequence cheng Bioinformatics 1 10-06-2010 08:36 PM
Bowtie Index Question Wei-HD Bioinformatics 5 07-29-2010 08:15 AM

Reply
 
Thread Tools
Old 02-07-2012, 08:26 AM   #1
nilmot13
Member
 
Location: Leeds, UK

Join Date: Jan 2011
Posts: 19
Default Simple Bowtie question

Hi,

I've been playing around with Bowtie on Galaxy platform and I would like to ask some very simple question really to confirm my understanding to the Bowtie parameters.

My aim for the mapping is to get the best sequencing alignment to a custome genome.

These are the parameter I've used:

Single-end
-n = 2
-e = 70
--nomaqround = round to nearest 10
-y
-a
--best
--maxbts = 800
--strata

Will the parameters I've set gave me a single alignment per sequence read and also ensure it's 'uniqueness/stringency'?

If anyone could help that would be fantastic.

Many thanks
nilmot13 is offline   Reply With Quote
Old 02-07-2012, 04:54 PM   #2
tianyub836
Member
 
Location: China

Join Date: Apr 2011
Posts: 23
Default

try option "-m 1"
tianyub836 is offline   Reply With Quote
Old 02-07-2012, 05:04 PM   #3
770349688
Junior Member
 
Location: china

Join Date: Feb 2012
Posts: 1
Default TMM normalization

Hello:
it is said that TMM is a normalization in RNA-seq,and I have learn it for a long time。But I donĎt understand it so far。Can you tell me how does it perform?And I would be gratefull if you could give me an example ?
770349688 is offline   Reply With Quote
Old 02-08-2012, 01:19 AM   #4
nilmot13
Member
 
Location: Leeds, UK

Join Date: Jan 2011
Posts: 19
Default

Quote:
Originally Posted by tianyub836 View Post
try option "-m 1"
Hi,

I didn't use -m = 1 beause it will reduce the mapping coverage of some reads. E.g in the ribisomal RNA region there gene cluster are very highly conserved and our lab is still interested about its coverage.

What I was really looking for is, would the parameter I've set gave me a stringent mapping result? i.e. less mismatches and poor quality reads being mapped.

Cheers
nilmot13 is offline   Reply With Quote
Old 03-07-2012, 01:56 AM   #5
oxydeepu
Member
 
Location: bangalore,india

Join Date: Jul 2011
Posts: 41
Default A simple question about bowtie output

Hi all..

i have a small doubt. In the output of bowtie, how to figure out mismatches. +ve is fine., but in negative there is reverse complementation thats what is giving me a confusion.
Please help as soon as possible.
Thank you in advance.

Deepak
oxydeepu is offline   Reply With Quote
Old 03-07-2012, 02:03 AM   #6
nilmot13
Member
 
Location: Leeds, UK

Join Date: Jan 2011
Posts: 19
Default

Hi Deepak, what exactly are the problems that you are encounting?
nilmot13 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 09:40 AM.


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