Question: Bwa Aln 5'-End Soft-Clipping
gravatar for GPR
7.5 years ago by
GPR320 wrote:

Hello, I am working my way through an RNA-editing flowthrough and am now trying to clip-off the Illumina 5'-end adapters. I have aligned my paired-end reads with bwa samse -n 4, and next have added headers, sorted and removed duplicates, to find out that I need to remove the Illumina adapter sequences. I am about to use the Picard tool CleanSam.jar to perform a soft-clip of my data and would like to ask the community is such "cleaning" is appropriate for my purposes: again trimming of 5'-end Illumina adapters. Thanks in advance!

ADD COMMENTlink written 7.5 years ago by GPR320
gravatar for John St. John
7.5 years ago by
John St. John1.2k
San Francisco, CA, Cancer Therapeutics Innovation Group
John St. John1.2k wrote:

I usually remove adapters prior to mapping. I wrote a program called SeqPrep to do this. More reads will map when you trim before mapping, also my program uses the overlap between the pair of reads as well as matches to the adapter (you are doing pairs, right?). That said if you are not doing pairs of reads, do not use SeqPrep, it is only for paired reads. If you are using single ended reads any tool that does adapter trimming (think there is one in GATK as well as Picard's CleanSam). The lab I used to be in used Trimmomatic for adapter trimming. There are a ton of programs though that will work reasonably well on single-ended reads if that is what you have.

ADD COMMENTlink written 7.5 years ago by John St. John1.2k

Thanks so much, useful recommendation

ADD REPLYlink written 7.5 years ago by GPR320
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1929 users visited in the last hour