SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
[Bowtie2] CIGAR string calculation. Coryza Bioinformatics 6 03-14-2014 02:12 AM
Amazon Web Services: Tophat cannot find Bowtie2 even though Bowtie2 is in path Dr_Funkenstein RNA Sequencing 4 03-06-2014 04:55 AM
how to calculate MAPQ value in the tophat feixue1039 Bioinformatics 5 09-05-2013 03:05 PM
Interpreting MAPQ values biobee07 Bioinformatics 6 08-13-2013 12:42 PM
How is the p-value for MAPQ calculated in Bowtie2? Cameron MacPherson Bioinformatics 0 03-22-2013 07:17 AM

Reply
 
Thread Tools
Old 03-08-2015, 05:40 AM   #1
KyuzoSeq
Junior Member
 
Location: Turkey

Join Date: Mar 2015
Posts: 9
Default Bowtie2 and MapQ Calculation

Hi

I'm going to ask about bowtie2's mapping quality calculation for ChIP-Seq data. For sequences which align multiple locations Bowtie2 gives the MAPQ 30. It says for multiple aligned sequences, mapq is not meaningful, so why it gives 30 for all? Other thing it gives 30 for completely aligned sequence and sequences which has indels too. What is that 30 and why?

I attached a png file about my ask.
Attached Images
File Type: png file.png (36.7 KB, 8 views)

Last edited by KyuzoSeq; 03-08-2015 at 05:47 AM.
KyuzoSeq is offline   Reply With Quote
Old 03-08-2015, 09:56 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,882
Default

http://biofinysics.blogspot.com.au/2...pq-scores.html

Here is added discussion on biostars: https://www.biostars.org/p/110958/
GenoMax is offline   Reply With Quote
Old 11-02-2015, 08:55 AM   #3
Gonza
Member
 
Location: Ithaca, NY

Join Date: Mar 2013
Posts: 78
Default Chip-seq mappers

Hello Geno Max and all,

I am inquiring about Chip-seq mapping software. I am using the Arabidopsis TAIR10 genome as a reference.

I have used for RNA-seq tophat2 to account for splicing, which is something we are not really concerned when doing Chip-seq (right ?).

I have read that people uses Bowtie, and SOAP, SOAP2 when mapping against TAIR10. Before i get my Illumina reads, I'd like to know if you have an idea as to which mapper would be best?

Many thanks in advance.

Cheers

G

Last edited by Gonza; 11-02-2015 at 09:01 AM.
Gonza is offline   Reply With Quote
Old 11-02-2015, 09:05 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,882
Default

@Gonza: You could use any modern aligner for the alignment. You don't need to worry about splicing with ChIP-seq. I like @blancha's advice in a recent thread. Pick an aligner and become familiar with its myriad options (most have many) rather than hopping from one aligner to next.
GenoMax is offline   Reply With Quote
Old 01-31-2016, 12:51 PM   #5
Gonza
Member
 
Location: Ithaca, NY

Join Date: Mar 2013
Posts: 78
Default

Hi GenoMax

I am mapping Illumina reads to a indexed genome with soap Version: 2.21 http://soap.genomics.org.cn/soapaligner.html

It seems to work fine, the only problem is that i cannot find the % of mapped reads. I am used to tophat that tells you % of mapped/unmapped. Is there a way to see % of mapped reads when using soap?

Script:
$soap -r 0 -M 1 -a sample1.fq -D genome.fa.index -o SOAP.sample1

Thanks
Gonza is offline   Reply With Quote
Old 01-31-2016, 06:11 PM   #6
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,882
Default

@Gonza: I am not familiar with soap but have you looked to see if soap writes any log files? Those may contain the information you are looking for.

Otherwise you could use "samtools idxstats" command or a package like Qualimap to find this information.
GenoMax is offline   Reply With Quote
Old 02-01-2016, 01:08 AM   #7
Gonza
Member
 
Location: Ithaca, NY

Join Date: Mar 2013
Posts: 78
Default

thanks, will look into Qualimap
Gonza is offline   Reply With Quote
Reply

Tags
alignment, bowtie2, chip-seq, mapq, sequence

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 03:10 PM.


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