Per sequence GC content still not good after trimming
1
0
Entering edit mode
28 days ago
zouini6 • 0

I am analyzing NGS data for a prokaryotic genome and FastQC initially failed "Per sequence GC content". I then ran Trimmomatic on my Fastq files and after running FastQC on the trimmed paired dataset, the GC content is still not good. could someone tell me what is the reason for this?

The first .jnp shows the result I got before trimming and the second one is for after trimming.

Thank you

GC content. trimmomatic. • 314 views
ADD COMMENT
1
Entering edit mode
28 days ago
GenoMax 107k

Is the organism you are working on GC rich? At this point in process I suggest that you move along with the analysis. If concerns crop up in dowstream analysis then you could revisit this to dig deeper.

There are no absolute failures with FastQC. Their appearence (red X) indicate that your data lies outside the defined bounds (which are for normal genomic sequencing).

ADD COMMENT
0
Entering edit mode

I don't know if it's a rich GC organism. I am moving along with the analyses but I just wanna know what are the potential interpretation for this because I'm supposed to comment and each step (with the results) for my assignment.

Thank you

ADD REPLY
0
Entering edit mode

A simple naive explanation can be that you have fragments in your data that contain more GC that most other fragments. That could simply be a by product of creation of libraries. At an extreme it could represent contamination by an unrelated genome.

ADD REPLY
0
Entering edit mode

Thank you so much

ADD REPLY

Login before adding your answer.

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