More posts like this
-
Paired-End Reads and Variant Calling
written by ocortes •Dear colleagues, I have the fastq files of paired end sequences of 42 animals but when I run bwa mem to generate the .sam files …
-
Read showing Insertion of 41 bases when aligned with BWA
written by kspataHi All, I aligned reads of a viral vector using both Bowtie2 local and BWA-MEM. 1300 reads are mapped to this position. I observed that …
-
Alternative for unsupervised clustering in scRNA Data
written by JohnHi dear one, As far I have seen in literatures all the data of single cell RNA seq have been subgrouped by unsupervised clustering by …
-
alternative for molquest, for inter-operonic region finding
written by JohnHi dear one, I want to find the interoperonic regions with any gene finder, my guide suggested molquest GUI software. I treied all the way …
-
BWA MEM vs BOWTIE2 , which is best?
written by JohnHello dear one, for my new project I would need sequence alignment tool to analyse RNA-seq data. I found that bowtie2 is a tool used …
-
MAPQ values are 60 in bwa mem
written by JohnHi dear one, I am getting lots of MAPQ values as 60 in bwa mem? I have ran the same read (paired end) in bowtie2 …
-
bwa mem has more MAPQ values as 0 and 60
written by JohnHi dear one, Why maximum no. of the MAPQ values are 0 and 60 in bwa mem? I have ran the same read (paired end) …
-
How to calculate bowtie2 run time using python
written by JohnHello dear one, I am new to programming! I want to calculate the run duration of bowtie2 and bwa. 1. First I don't know (the …
-
bwa mem hard clipping
written by ARB •Hi all, I apologize if this question has been posted or answered earlier. But I did not seem to find a reasonable explanation and hence …
-
bwa mem mismatch option
written by Varun GuptaHi Everyone, I am using bwa mem to align my paired end read data. My reads are 150bp in length(paired end). I am using bwa …
for bwa : http://bio-bwa.sourceforge.net/bwa.shtml
Thank you! But In this page, it is only given that about the penalty values! I am intended to compare bowtie2 and bwa-mem so I want to compare the penalty model of both! Is there anything related to how penalty model build?
The BWA-MEM algorithm has not been published yet. (bottom of the link http://bio-bwa.sourceforge.net/bwa.shtml you didn't read )
But your question is bad because penalty is penalty, you are looking to alignment algorithm/model maybe ...
It's published on ArXiv but not in a peer-reviewed journal (and it probably never will).
oh didn't know that i told that on their web site information, sorry for the mistake. Do you know why it's not reviewed ?
AFAIK it was submitted, and Heng Li did not agree with (unreasonable) reviewers comments and then decided not to proceed.
Also doesn't really matter. You can cite and read the preprint.
Is it only me who is surprised that we are discussing whether BWA-mem is published and peer-reviewed or not? It looks absolutely ironical given the usage of BWA-mem in the early years before the explosion in the numbers of different kind of aligners, and not to discount the fact that there are still die-hard fans of BWA-mem (including me).
that's not changing my point view about this tools at all , i was just curious about it :) It show once again that publication evaluation system got bad points :)
Searcing for the word "penalty" in the respective manuals gives you all the details
http://bio-bwa.sourceforge.net/bwa.shtml
http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml
thanks, i know them, but i want penalty model to compare them
bwa mem and bowtie2 are two tools which have been used extensively in the field. Are you sure your comparison is still necessary? It might, of course, be interesting for your own understanding of the tools.
thanks for all your answers and discussions, Even though these two are best tools I just want compare and understand both the tools !
Still my question is,
ex: Are Gap open penalty in bowtie2 and bwa-mem meant to be same (in terms of value)? my guide suggested to have a look at penalty models of both the tools then set (change accordingly) a default parameters for both of them!
I didn't get the answer!
Thanks in advance for the future replies!
Do you compare the ripeness of bananas and oranges by comparing their colors, because both are fruits?
BWA-mem and bowtei2 are different algorithms, and for sure there will be differences in how the gap penalty affects one algorithm vs other. Moreover, gap penalty can be set depending upon context too (and bowtie2 provides more parameters to tweak various penalties). For routine exercises, it is recommended to use the default settings. However, if you have specific needs (like your PI wants that), then please read the manuals and relevant materials first, and ask specific Qs, like: "I want to choose X in bwa vs Y in bowtie in context Z because I thought about A, B and C -- is is justified or what are your views?"
Hi, Even I'm looking for a optimal aligner and variant caller for analysing human exome data. I performed alignment with different aligners BWA-mem, Bowtie2, Cushaw3 and Novoalign I don't have a accurate answer to say which one is the best.
Have you tried increasing the gap penality to reduce the gaps ?