What is default paprameter of BWA mem mapping
Please Discuss the default parameter setting of BWA mem.
I used bbmap and BWA mem with default setting to map raw reads to CDS transcripts and I found that BWA is able to map more number of reads in comparison to bbmap.
I know in bbmap default setting is 76% identity but not aware about BWA mem.
Thanks
• 9,343 views
•
link
updated
by
Jenez
•
written
by
Bioinfonext
•
written
by
Bioinfonext
2 answers
I would take a guess and say that the default settings for each parameter can be found on the documentation page for bwa mem here.
See the values within [ ] for the default parameters.
Log in to answer this question.
More posts like this
-
Bwa mem different alignment results for the same reference genome
written by EarthShaker-EAC •
-
Calling variants after using different read alignment programs
written by ManuelDB
-
mismatch setting in bowtie2
written by Bioinfonext
-
Cufflinks produces empty transcript.gtf files
written by william.mcentegart •
-
comparison of BWA MEM and hisat2
written by guido.leoni •
-
De novo assembled transcript mapping to genome scaffold
written by Bioinfonext
-
To extract unmapped pair end reads
written by Bioinfonext
-
How to extract unmapped Paired end reads from BWA output sam file
written by Bioinfonext
-
Error while running BWA
written by Bioinfonext
-
GSNAP cammand for mapping pair end read to CDS transcripts and to extract unmapped pair end reads
written by Bioinfonext
By default, BWA produces multiple alignments per read, while BBMap does not unless you add the "ambig=all" flag. Be sure that you are counting correctly - the number of alignments produced is not the same as the number of reads that map. The difference is particularly large for transcriptome mapping because of all the duplication due to isoforms.