SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-seq results interpretation - help needed rebrendi Bioinformatics 11 09-02-2012 03:36 AM
BWA produces odd alignment results dandyrilla Bioinformatics 2 11-27-2011 11:28 PM
SHRiMP alignment results for colorspace miRNA dsidote Bioinformatics 2 11-10-2011 09:32 AM
Help interpretation indel labmtl Illumina/Solexa 0 03-29-2011 07:01 AM
view alignment and assembling results with UCSC genome browser? dicty Bioinformatics 1 03-04-2011 07:38 AM

Reply
 
Thread Tools
Old 04-27-2013, 01:51 AM   #1
carolW
Senior Member
 
Location: US

Join Date: Apr 2013
Posts: 103
Default interpretation of alignment results

Hi,
How to interpret the results of a program alignment like BWA or bowtie? In general, is it the percentage of different reads, for ex less than 2% or 5%? How about in details?

Any tutorial, review and documentation is welcome.

Look forward to your reply,

Carol
carolW is offline   Reply With Quote
Old 04-27-2013, 12:51 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

No one can answer this without knowing what your experiment is. And only you know that.
swbarnes2 is offline   Reply With Quote
Old 04-28-2013, 12:23 AM   #3
carolW
Senior Member
 
Location: US

Join Date: Apr 2013
Posts: 103
Default

As a newbie, I have started with the public data before analysing private data. So how to interpret already on the public data like 1000 genomes and then, focus on the private data?
carolW is offline   Reply With Quote
Old 04-28-2013, 01:35 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

To reiterate what swbarnes2 said, you have to ask a coherent question to receive a useful answer. "Analyze" is a broad term, you need to specify what sorts of things you're interested in looking at. Given that you mentioned 1000 genomes data, are you interested in looking at SNPs?
dpryan is offline   Reply With Quote
Old 04-28-2013, 07:57 AM   #5
carolW
Senior Member
 
Location: US

Join Date: Apr 2013
Posts: 103
Default

yes, but already after the alignment, how can I evaluate that the results are good? for ex, I got the following result for one of the 1000 genomes after running bwa aln. I can think max of 9 mis-alignments could be considered as negligeable but how to evaluate in other cases? Is it 9 out of 225 bp read or number of seq 109811 should be considered?

[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_aln_core] 109811 sequences have been processed.
carolW is offline   Reply With Quote
Old 04-30-2013, 06:27 AM   #6
Inma
Member
 
Location: Spain

Join Date: Feb 2013
Posts: 25
Default

I think it's a good idea visualizing data from BWA in IGV (http://www.broadinstitute.org/igv/) or tablet (http://bioinf.scri.ac.uk/tablet/) and check the sequence coverage or the way your reads are mapping in the genome.
Inma is offline   Reply With Quote
Old 04-30-2013, 06:31 AM   #7
carolW
Senior Member
 
Location: US

Join Date: Apr 2013
Posts: 103
Default

Quote:
Originally Posted by Inma View Post
I think it's a good idea visualizing data from BWA in IGV (http://www.broadinstitute.org/igv/) or tablet (http://bioinf.scri.ac.uk/tablet/) and check the sequence coverage or the way your reads are mapping in the genome.
Thanks for your input.

I can do this as the data set is small but would it be a good idea to do so if the data set is large?
carolW is offline   Reply With Quote
Old 04-30-2013, 07:11 AM   #8
Inma
Member
 
Location: Spain

Join Date: Feb 2013
Posts: 25
Default

I did this in a 5 Mb bacterial genome, mapping 301453 reads against the sequence. I don't know if this programs are restricted to a small set of data.
Inma is offline   Reply With Quote
Old 04-30-2013, 08:48 AM   #9
carolW
Senior Member
 
Location: US

Join Date: Apr 2013
Posts: 103
Default

And then, how to evaluate quantitatively the sequence coverage checking and acceptable number of mis-alignments?
carolW is offline   Reply With Quote
Old 04-30-2013, 07:01 PM   #10
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quote:
Originally Posted by carolW View Post
And then, how to evaluate quantitatively the sequence coverage checking and acceptable number of mis-alignments?
You can get an estimate of coverage by looking at the # of aligned reads, length of reads, versus the length of your genome.

There is no magical value to tell you if reads mis-aligned. If reads mis-aligned, it's because your genome is repetitive, or it doesn't match what you actually sequenced. Neither one is the fault of the software, and neither one is something you really need to worry about much, let alone try to use as a general guide to how well your experimetn worked.

What is more helpful is to know the % of reads that aligned to your target. If the benchwork people screwed up the library prep, or misinformed you as to what the sample actually is, an abnormal alignment % will tell you that.

What % is expected depends on what the experiment is, and how good the people doing the benchwork are.
swbarnes2 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 07:29 AM.


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