converting Bam to fastq while removing clipping(hard/soft clip bases)
1
0
Entering edit mode
5 weeks ago
eli_bayat ▴ 50

Hello,

I want to do some analysis and my raw data is paired-end reads fastq files. So far:

  1. I used BWA mem to convert them to Sam file

  2. then used samtools to convert to BAM file.

My next step is to convert bam file to fastq file and for that, first I need to exclude the clipping (hard/soft clip bases). Is there any tool that does this while converting to fastq file?

So far, I looked into few tools below:

  1. samtools: I couldn't find any option that remove clipping while converting to fastq.
  2. picard: --CLIPPING_ACTION seems to do this if I set it to 'X'. But I am not sure.
bam_to_fastq hard_clipping Soft_clipping tools remove_clipping • 309 views
ADD COMMENT
1
Entering edit mode
5 weeks ago
GenoMax 108k

Hard clipped bases are not in the BAM file.

As for soft-clipped reads use the answers here: Remove Soft Clipped Bases

ADD COMMENT
0
Entering edit mode

Great, that's what I wanted. Thanks,

ADD REPLY

Login before adding your answer.

Traffic: 1578 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6