This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bwa/Samtools Soft-Clipping To Retrieve Back Whole Read

Hi,

(I've been googling for this and couldn't find the answer)

Am I right in thinking that the default bwa/samtools instructions in the website will hard-clip the reads in the SAM/BAM files?

In the usual pipeline for single fragment reads: bwa aln ... bwa samse ... samtools view ... samtools sort ... samtools index ...

How can I end up with a final SAM/BAM file that contains the whole of the read/qualscore even if it's not fully matched?

Cheers

bwa samtools
• 5,374 views • link
updated by Ketil Ketil • written by 2184687-1231-83- 2184687-1231-83-

No, it is soft clip.

• 0 views • link
written by lh3 lh3

1 answer

I've been using bwa, and I see soft-clipped reads in the results. I doubt that samtools will change this. As far as I can tell, the BAM files also contain the whole reads.

• 1 views • link
written by Ketil Ketil

Log in to answer this question.

More posts like this