Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA Soft Clipping Bio.X2Y Bioinformatics 11 03-09-2015 09:08 AM
Bismark and unique matches mixter Bioinformatics 3 12-12-2012 04:06 AM
unmapped vs soft-clipping CNVboy Bioinformatics 2 04-08-2012 10:41 PM
Reference that matches vcf? brofallon Bioinformatics 5 02-24-2012 06:36 PM
bwasw Bgansw Bioinformatics 2 02-20-2012 04:00 PM

Thread Tools
Old 11-07-2013, 03:17 PM   #1
Junior Member
Location: California

Join Date: Dec 2009
Posts: 5
Default bwasw soft clipping of exact matches

I am using bwasw to align reads from a MiSeq run to a very short reference sequence. For some reason, bwa is soft-clipping (assigning S in the cigar string) the leading and trailing bases of many or most of my reads. (All of my reads seem to get the first base soft-clipped, after that there is a distribution of soft clipping that drops off toward the center of the read).

This is occurring despite the fact that the bases which are getting soft clipped (at least the ones I checked by scanning the sam) are EXACT MATCHES to the reference. They are not consistently surrounded by mismatches (they occur even in reads with an edit distance of 0), nor are they low-quality (prior to assembly, I ran a very stringent quality filter, so every base in every one of my reads is phred 36 or higher).

Possibly pertinent: all of my reads are the same length as my reference +/- 2 nts. (This is a ribozyme population deep sequencing project).

Does anyone know why bwa might be acting this way or how to get around it? Here's my bwa command:

bwa bwasw -M -NM myreference.fasta myfiltered.fastq > myoutput.sam

Alternatively, is there a way to tell samtools to ignore soft-clipping? Or a tool to convert the cigar strings to M/I/D only? The soft-clipped bases are artificially reducing the coverage in my mpileup, and I am worried how they will affect other downstream modules I want to use.

Thanks for any insight!
petrie is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:15 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO