Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Two peaks on FastQC plot "Per sequence GC content"

    Hi,
    I just got illumina DNA genome re-sequencing data. All the items in FastQC reports passed but "Per sequence GC content". There are two peaks on the plot of "Per sequence GC content". The major peak centers around 40% GC content, while the minor peak centers around 70% GC content.

    I would appreciate it if you can explain to me how this happened and what I should do to correct it or discard the minor peak.

    Thanks in advance!

  • #2
    It suggests that maybe you have some kind of contamination.

    What %GC content are you expecting for the species you are sequencing?

    I would do adapter trimming/quality trimming and rerun FastQC afterwards to see whether that gets rid of the problem or not.

    Comment


    • #3
      Its probably the adapters. Do some trimming and it will go away.

      Comment


      • #4
        If the secondary peak is very sharp it's probably a specific contaminant - often something which is found by the overrepresented sequences module.

        If the peak is fairly sharp and not too far from your main distribution it could be long read through into adapters as suggested above.

        If the secondary peak is quite broad then it might be that you have contamination with a different species. You could use something like fastq_screen to check for other species you work with regularly, but this won't pick up other odd sources of contamination.

        Comment


        • #5
          The “Per base sequence content” and “Per base GC content” graphs should not show contamination of the adapters?

          Comment


          • #6
            Originally posted by MichalGordon View Post
            The “Per base sequence content” and “Per base GC content” graphs should not show contamination of the adapters?
            They might show some effects. If you have adapter dimers then you'll see the adapter sequence superimposed on the sequence content graphs. If your adapters have markedly different GC content than your library in general then you might also see an overall effect on the GC level.

            In the latest fastqc release there is a graph specifically to measure adapter content which will show exactly what proportion of the library is composed of read-through adapter which will illustrate this much better than trying to use sequence content plots.

            Comment


            • #7
              Thank you!

              Comment


              • #8
                Originally posted by simonandrews View Post
                They might show some effects. If you have adapter dimers then you'll see the adapter sequence superimposed on the sequence content graphs. If your adapters have markedly different GC content than your library in general then you might also see an overall effect on the GC level.

                In the latest fastqc release there is a graph specifically to measure adapter content which will show exactly what proportion of the library is composed of read-through adapter which will illustrate this much better than trying to use sequence content plots.
                I am having a similar problem with my run (2x150), As you can see there are two peaks in my run. I expect to have a 40% of GC content (bacterial genome) but I don know why did I obtain these two peaks.

                [PASS] Basic Statistics
                [PASS] Per base sequence quality
                [PASS] Per sequence quality scores
                [FAIL] Per base sequence content
                [FAIL] Per base GC content
                [WARNING] Per sequence GC content
                [PASS] Per base N content
                [WARNING] Sequence Length Distribution
                [WARNING] Sequence Duplication Levels
                [WARNING] Overrepresented sequences
                [WARNING] Kmer Content

                Oversequencing is probably not the problem because in fact I obtained less reads as expected. Could it be due to a adaptor problem? Any clue would be really appreciated
                Attached Files

                Comment


                • #9
                  I think it will be helpful if you could provide more information such as library type, input material, kit used for library prep and graphs from new version of FastQC.

                  Comment


                  • #10
                    Originally posted by nucacidhunter View Post
                    I think it will be helpful if you could provide more information such as library type, input material, kit used for library prep and graphs from new version of FastQC.
                    I updated FastQC to the 11.2 version and my error disappeared. I wonder it was an old version problem...

                    Comment


                    • #11
                      Originally posted by chariko View Post
                      I updated FastQC to the 11.2 version and my error disappeared. I wonder it was an old version problem...
                      The per base GC plot was removed in the latest version since it mostly replicated information which was in the per base composition plot. You should still be able to see the biased positions as a deviation in the composition of C or G content at the same positions, but it's possible it's not enough of a deviation to trigger a warning.

                      Comment


                      • #12
                        Originally posted by simonandrews View Post
                        The per base GC plot was removed in the latest version since it mostly replicated information which was in the per base composition plot. You should still be able to see the biased positions as a deviation in the composition of C or G content at the same positions, but it's possible it's not enough of a deviation to trigger a warning.
                        As you can see in the per base composition plot the C content goes down on position 5 (as seen in the per base GC plot before and goes up on position 9. I assume as the manual tells, the first 12 positions could be a selection bias.
                        I assume everything is OK then since the GC content in the specie s around 40%,


                        It was an Nextera MiSeq bacterial genome sequencing experiment.

                        Thank you very much for your help
                        Attached Files
                        Last edited by chariko; 08-19-2014, 01:40 AM.

                        Comment


                        • #13
                          Everything is ok

                          Comment


                          • #14
                            Hi!

                            I have two problems: one is two peaks in the per sequence GC-content and another is a weird profile which I'm attaching here.

                            We're trying out Agilent's SureSelect enrichment protocol for Exome-Seq and have just concluded our first run on samples that were already done before using Illumina's Nextera kit (so we have another run with which to compare our results). The first run was sequenced on the Illumina HiSeq while this run was done on a MiSeq. Also, the first run was a 100bp paired end run while this was 150bp paired end run. Anyway, upon running a QC on the Fastq files I got this weird profile for the per-sequence GC content. I had already removed the low-quality reads and trimmed the adaptors but that didn't change anything. The only thing that helped was trimming 25 nucleotides from each end of the reads. Since we lose a lot of information that way, I'd prefer not to do this and want to ask if anyone has seen anything like this. I have no idea what might cause this.
                            Attached Files
                            Last edited by Khillo81; 10-14-2014, 04:55 AM.

                            Comment


                            • #15
                              This is sometimes a sign of contamination, though if trimming the reads reduces it, that's a bit odd. Is this supposed to be human data? Human should peak around 50%, which does not correspond to either of your peaks. The most important question is what organism this is supposed to be, and what it's average GC% is.

                              Also, please post an insert-size histogram, which will help determine if the problem is caused by short inserts. You can get one quickly using BBMerge:

                              bbmerge.sh in1=read1.fq in2=read2.fq ihist=ihist.txt

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                New Genomics Tools and Methods Shared at AGBT 2025
                                by seqadmin


                                This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                                The Headliner
                                The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                                03-03-2025, 01:39 PM
                              • seqadmin
                                Investigating the Gut Microbiome Through Diet and Spatial Biology
                                by seqadmin




                                The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
                                02-24-2025, 06:31 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, 03-03-2025, 01:15 PM
                              0 responses
                              178 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 02-28-2025, 12:58 PM
                              0 responses
                              272 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 02-24-2025, 02:48 PM
                              0 responses
                              656 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 02-21-2025, 02:46 PM
                              0 responses
                              267 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X