Question: How to retrieve hard clipped reads
0
gravatar for teresayahu
9 weeks ago by
teresayahu0
teresayahu0 wrote:

Hi,

Is there a way to retrieve the part of reads that has been hard clipped and discarded in bam files?

Thanks!

alignment • 179 views
ADD COMMENTlink modified 9 weeks ago by d-cameron880 • written 9 weeks ago by teresayahu0

If you have the original fastq files you could fish them out of there. I assume you probably don't have them and hence the question?

ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by genomax37k

Exactly. I don't have the fastq files and only have bam file.

ADD REPLYlink written 9 weeks ago by teresayahu0

Then you have answered your own question. If they are not in the BAM then there is no way to get them back.

ADD REPLYlink modified 9 weeks ago • written 9 weeks ago by genomax37k
0
gravatar for d-cameron
9 weeks ago by
d-cameron880
Australia
d-cameron880 wrote:

Hard clipping removes the read base and quality information from that record in the BAM file. It must be recovered from another source (e.g. the original FASTQ file) as the information is not present in the BAM file.

Exception: if the BAM file contains multiple record for a single read, the hard clipped record can be rehydrated from the other record. For example, if your hard clips are due to split read alignments generated by bwa, the primary read alignment mapping location record will be soft clipped and will containing the required information. The hard clipping on the supplementary alignments can be converted to soft clips using the information from the primary read alignment.

My SV caller GRIDSS includes a (somewhat poorly named) utility program gridss.ComputeSamTags that, using the SOFTEN_HARD_CLIPS parameter, can do exactly that (you'll need to sort by read name first).

ADD COMMENTlink modified 9 weeks ago • written 9 weeks ago by d-cameron880
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: 809 users visited in the last hour