This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Difference in fastqc output after trimming with fastp vs afterqc

Hello, Upon trimming ddrad fastq files with fastp, the per base sequence content is still a 'red cross' but after using afterqc, the same is now yellow with an exclamatory mark, indicating it has improved. But upon using stacks, the RE site is not identified in the trimmed clean fastq file generated through afterqc. But the RAD tags are perfectly identified in the fastp cleaned fastq file. What might be the reason?

fastqc fastp ngs ddrad afterqc
• 3,354 views • link
written by ymberzal • ymberzal

How did you specify trimming with fastp? It seems the trimming may be done differently and afterqc has more default qc settings.

• 0 views • link
written by rfran010 • rfran010

Thank you for your reply.

$ fastp -i 1A_R1.fastq' -I 1A_R2.fastq --dont_overwrite -q 30 -u 40 --disable_length_filtering --cut_right -c -g -x -o fastp1A.fq -O fastp1B.fq

After I used STACKS to check for RAD cut sites with afterqc trimmed files, this is what I get: 4301656 total sequences 0 failed Illumina filtered reads (0.0%) 0 barcode not found drops (0.0%) 1518 low quality read drops (0.0%) 3839460 RAD cutsite not found drops (89.3%) 460678 retained reads (10.7%)

But fastp trimmed files pass through STACKS perfectly. I was wondering what might be the reason.

• 0 views • link
updated by GenoMax GenoMax • written by ymberzal • ymberzal

If you are following the STACKS pipeline you don't need to do separate trimming of the data. I see the following in manual

If you believe your reads may contain adapter contamination, process_radtags can filter it out.

Simply follow the STACKS pipeline.

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this