Question: How to compare the alignment quality of two or more mapping software
gravatar for ravi.uhdnis
4.4 years ago by
United States
ravi.uhdnis180 wrote:


I need to compare the performance of 2 or more alignment software (e.g. BWA, a new software X) while aligning a particular gene (e.g. HLA-A) against hg19 as human reference genome. Basically, i need to know sensitivity, Specificity, Accuracy, False Positive, False Negative of these software taking the same set of gene as input and against same reference human genome. Although BWA gives basic alignment summary and we can have a overall look via evaluating mapping file (.bam) using IGV but is there any other better way to look the overall alignment summary by different software of same genomic region of human genome after mapping by different software for their performance comparison ?. Looking for suggestions, Many thanks. 

ADD COMMENTlink modified 4.4 years ago by harold.smith.tarheel4.5k • written 4.4 years ago by ravi.uhdnis180
gravatar for Devon Ryan
4.4 years ago by
Devon Ryan95k
Freiburg, Germany
Devon Ryan95k wrote:

The only way to do this is to generate synthetic (i.e., fake) data with known results, align it with all of the tools, and compare the results. You can do this with tools like wgsim, art, etc. Here's an example from Heng Li using wgsim to compare a number of common aligners (N.B., the comparison is very very outdated).

ADD COMMENTlink written 4.4 years ago by Devon Ryan95k
gravatar for harold.smith.tarheel
4.4 years ago by
United States
harold.smith.tarheel4.5k wrote:

BBMap 'randomreads' command lets you generate synthetic data of specified read lengths, error types and rates, etc for benchmarking. And you can use 'samtoroc' from the same package to generate ROCs for comparing aligner performance.

ADD COMMENTlink modified 4.4 years ago • written 4.4 years ago by harold.smith.tarheel4.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1011 users visited in the last hour