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?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
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.
-
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?
Comment
-
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.
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 08:47 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Today, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
57 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
53 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Comment