Question: how to relate BWA input parameters to get minimum coverage and identity
0
gravatar for gm
3 months ago by
gm0
gm0 wrote:

Hi all and thanks in advance for any help you will kindly provide.

I'm aligning some short sequences against a few databases by means of BWA. I need to recover bot identity and coverage in a fast way at the end of the alignemnet. I read that I can use the M field in CIGAR form SAM and further parameters from SAM to distinguish between match/mismatch. Can you give more precise information on this point?

Is there easy way to modify BWA to have identity and coverage directly (as it is done in DeconSeq, whose code I tried to analyse but with little success)?

Is there any settings for input parameters to BWA to obtain some desired coverage and identity?

Thanks a lot, Giovanna

alignment • 115 views
ADD COMMENTlink modified 3 months ago • written 3 months ago by gm0
0
gravatar for Jennifer Hillman Jackson
3 months ago by
United States
Jennifer Hillman Jackson23k wrote:

Hello,

There isn't a public tool (that I am aware of) that takes an input BAM/SAM dataset and then outputs per-query-read-alignment percent identity and coverage.

BWA has many options for filtering alignments. Open up the advanced settings to review. These options correspond to the same options found in the BWA manual.

Another choice is to post-alignment filter your current results. Try these two first: Filter mapped reads on MD tag string or BAM filter?

Thanks! Jen, Galaxy team

ADD COMMENTlink written 3 months ago by Jennifer Hillman Jackson23k
0
gravatar for gm
3 months ago by
gm0
gm0 wrote:

Thanks for your kind support. I'll have a try.

ADD COMMENTlink written 3 months ago by gm0
Please log in to add an answer.

Help
Access

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