Per base GC content with a shoulder after adapter removal and trimming
1
0
Entering edit mode
3.2 years ago
serpalma.v ▴ 70

Hello!

The following graph shows a fastq file (WGS) after trimming (sample comes from mouse). Adapters and low quality bases were removed with trimmomatic. The module "Overrepresented sequences" was flagged with a "PASS".

I have gone through other threads. This one in particular help me link1, but I am not sure what it means when a curve has "Small shoulders".

This other thread Link2 also gives useful hints, but again the meaning of "a broad secondary peak" can be ambiguous.

Based on the figure below, could you help me decide if the sample is acceptable for downstream analyses, and if I should apply some other tool to correct this issue?

EDIT: links fixed and sample nature added.

enter image description here

fastqc GC content • 1.1k views
ADD COMMENT
0
Entering edit mode

What is the nature of the sample? That is, is it RNA-seq, ChIP-seq, etc.?

ADD REPLY
0
Entering edit mode

I added that to the post. It is WGS data

ADD REPLY
0
Entering edit mode

A bit of PCR bias could produce this. As others have mentioned, you might as well continue the analysis.

ADD REPLY
0
Entering edit mode

Your two links are the same one.

I don't see a second peak on your graph (at least not a clear one). The parabol is not perfect but I would personally not worry about this for the next step

You work on mouse, if you got human contaminations, reads from human will not map the mouse reference genome and will stay unmapped

ADD REPLY
0
Entering edit mode
3.2 years ago

Shoulders in this graphics are usually compatible with at least two issues

  1. Presence of contaminant DNA

  2. Poor quality in the sequencing

The point2 can be solved at least in part after a quality trimming. If the shoulder persist after the quality trimming a contamination fo your samples can be likely to be happening. A rapid metagenomic analysis with programs like Kramer can maybe help you in getting an answer

ADD COMMENT
0
Entering edit mode

The sample was already processed with trimmomatic, including adapter removal.

FastQC did not detect adapters.

The quality Phred score was above 28 for all bases.

ADD REPLY

Login before adding your answer.

Traffic: 2916 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