Question: primer trimming tools
0
gravatar for J.F.Jiang
2.9 years ago by
J.F.Jiang710
China
J.F.Jiang710 wrote:

Hi all,

Due to the PCR based amplification, we need to trim to primers off the PE reads.

However, it seems that the trimming process might bring in low quality in tails, especially that some reads might even fall into Q15-Q20, which might lead to low Quality using GATK variant calling.

`$cutadapt -q 10 -g file:$fprimer -G file:$rprimer --minimum-length 20 $read1 $read2 -o $read1_trim -p $read2_trim`;


The command listed above is the one I used in my calling pipeline, although I can adjust the q to 30, there would be quite a lot reads to be removed.

Is there any smarter method to trim the primers?

Thanks.

Junfeng

primer trimming • 1.8k views
ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by J.F.Jiang710
0
gravatar for Brian Bushnell
2.9 years ago by
Walnut Creek, USA
Brian Bushnell15k wrote:

Doing adapter-trimming will not somehow reduce the quality of your reads.  Variant callers are expected to take quality into consideration when making calls, and good mapping programs are robust against low quality tails on reads; the extra length improves mapping confidence.

I don't recommend trimming to a quality threshold above 12 for mapping.  If cutadapt has trouble finding adapters in low-quality reads without aggressive quality-trimming, I suggest you use BBDuk instead, which does quality-trimming after adapter-trimming, as it is does not need quality-trimming to increase sensitivity.  Trimming for quality prior to trimming adapters reduces the number of adapter bases present, ultimately making them harder to detect.

Edit - it looks like the OP was talking about targeted amplicon sequencing and trimming the 3' primer, in which case my answer is irrelevant.

ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by Brian Bushnell15k
0
gravatar for J.F.Jiang
2.9 years ago by
J.F.Jiang710
China
J.F.Jiang710 wrote:

After trying and searching references, I found the most efficient way is to align the sequences using aligners first, e.g., bwa-mem, then using GATK ClipReads later to remove the matched primer sequences, since primers can help to align the reads to ref genomes, while trimming the reads first will generate a lot of "bad" quality reads that can result in lower depth of coverage or waste of inputs.

And another question is that I calculated the Ti/Tv ratio for target sequencing, e.g., BRCA1 and BRCA2 genes, and obtained a varied value from 2.X to 8.X for different samples. Is this normal or not?

ADD COMMENTlink written 2.9 years ago by J.F.Jiang710

We also adopted the workflow of BWA-MEM alignment first and then soft-clip primer afterwards. Since ours are nested amplicons, GATK ClipReads cannot properly handle the primer clipping. We then developed and use BAMClipper (Scientific Reports 7:1567). The bonus is that errors in primer sequences (synthesis and/or sequencing) are tolerated.

ADD REPLYlink written 16 months ago by Tommy Au70
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: 1182 users visited in the last hour