I want to do some analysis and my raw data is paired-end reads fastq files. So far:
I used BWA mem to convert them to Sam file
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:
- samtools: I couldn't find any option that remove clipping while converting to fastq.
- picard: --CLIPPING_ACTION seems to do this if I set it to 'X'. But I am not sure.