Question: Bwa Alignment Length
0
gravatar for Newvin
8.6 years ago by
Newvin340
Newvin340 wrote:

After performing a transcriptome assembly, I am using BWA to align the reads to my assembled contigs. From the resulting SAM file, I am able to calculate the # of reads that map to each contig, but I would like to normalize this abundance measure by using the length of each alignment.

My naive question: is the length of each alignment equal to the length of the read that is being aligned, or does BWA sometimes align only part of a sequence query to the target?

alignment bwa • 2.7k views
ADD COMMENTlink modified 4.2 years ago by Biostar ♦♦ 20 • written 8.6 years ago by Newvin340

What sequencing platform are you using? Your read lengths vary that much?

ADD REPLYlink written 8.6 years ago by Damian Kao15k

These are Illumina short reads. The variation in read length occurs because we used an adaptive algorithm to trim the reads based on quality score prior to assembly.

ADD REPLYlink written 8.6 years ago by Newvin340
2
gravatar for Jeremy Leipzig
8.6 years ago by
Philadelphia, PA
Jeremy Leipzig19k wrote:

BWA is performing a local alignment. If you see an 'S' in the CIGAR string of the resulting SAM file (e.g 40M36S), that 3' soft clipped end is not part of the alignment.

ADD COMMENTlink written 8.6 years ago by Jeremy Leipzig19k
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 2.3.0
Traffic: 1587 users visited in the last hour