Question: Why ExAC used bwa aln instead of mem?
gravatar for abascalfederico
3 months ago by
abascalfederico970 wrote:

Anyone knows why in ExAC they used bwa aln instead of mem? I am not sure of the read lengths but somewhere said something about 150 bps x 2. Does it make sense to use bwa aln here? I though for long reads mem was the best choice

The actual parameters used in ExAC were:

bwa aln q 5 -l 32 -k 2  -o 1 ...
bwa sampe -T -P ...
bwa exac • 226 views
ADD COMMENTlink written 3 months ago by abascalfederico970

If you can cite your source for this, that'd be great. I'm sure Daniel MacArthur's team would be glad to answer your question by the way - you could just send them an email!

EDIT: Found the source: (supplementary in the paper

ADD REPLYlink modified 3 months ago • written 3 months ago by Ram12k

I'll send them an email, thanks for the suggestion. Regarding the use of bwa aln, I found it in the supplementary material of: Lek, M., Karczewski, K. J., Minikel, E. V., Samocha, K. E., Banks, E., Fennell, T., ... & Tukiainen, T. (2016). Analysis of protein-coding genetic variation in 60,706 humans. Nature, 536(7616), 285-291.

ADD REPLYlink written 3 months ago by abascalfederico970
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: 781 users visited in the last hour