SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
FastQC: 3' bias (weird peaks) in Kmer content graphics Ann7 RNA Sequencing 2 06-19-2013 06:31 AM
FastQC: odd kmer content zshuhua Introductions 3 05-13-2013 08:36 PM
kmer content warning in FastQC vallejov RNA Sequencing 0 04-05-2013 11:10 AM
weird kmer-content peak in RNA-seq data kareldegendt Bioinformatics 2 08-22-2012 04:45 PM
Kmer Distribution Problem cyyuan Illumina/Solexa 3 05-05-2012 11:13 PM

Reply
 
Thread Tools
Old 07-12-2013, 11:12 AM   #1
feralBiologist
Member
 
Location: UK

Join Date: Jun 2011
Posts: 61
Default Weird kmer distribution (using fastqc)

I ran FatsQC on a set of 4 Illumina PE exomes and got really weird graphs for the kmer distribution. I've never seen this pattern before - does anyone know what is going on and what to do about it?







feralBiologist is offline   Reply With Quote
Old 07-12-2013, 03:25 PM   #2
fkrueger
Senior Member
 
Location: Cambridge, UK

Join Date: Sep 2009
Posts: 597
Default

I am not entirely sure about what is going on since I don't have the full length Illumina adapter sequences handy, but my guess is that the plot shows an adapter or primer dimer of some sort. Seeing that GGAAG is part of the start of the Illumina adapter (GATCGGAAGAGC..., wihtout the A-tail at the start) and the spacing is always the same this might simply be full-length adapter dimers. Such sequences should not align anywhere, so you probably don't need to worry about them. This doesn't mean that you don't have to trim your sample at all, this is difficult to say without seeing the rest of the plots.
fkrueger is offline   Reply With Quote
Old 07-14-2013, 06:07 AM   #3
feralBiologist
Member
 
Location: UK

Join Date: Jun 2011
Posts: 61
Default removing over-represented sequences does not help

@fkrueger: I found that two TruSeq adaptors were listed among the over-represented sequences in the fastqc report so I removed them for one of the samples and re-ran fastqc. The graph changed somewhat but there are still non-random peaks (below). Any ideas?


feralBiologist is offline   Reply With Quote
Old 07-14-2013, 06:20 AM   #4
fkrueger
Senior Member
 
Location: Cambridge, UK

Join Date: Sep 2009
Posts: 597
Default

I don't know how you removed the adapters, but the k-mers could still result from adapter dimers. Trim Galore for example only removes the entire adapter sequence AGATCGGAAGAGC but only if it starts with the A from the A-tailing process. In your case the sequence seems to start with GATC(GGAAG) and thus wouldn't be affected by adapter trimming.

The CACAC and ACACA is often a highly occurring sequence in MeDIP-Seq or the like. I wouldn't worry too much about the plot for the moment, and only come back to it if you find that you've got an unusually low mapping efficiency or similar.
fkrueger is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:27 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO