Explanation for kmer enrichment after skewer trimming of the RNA-seq reads
1
0
Entering edit mode
2.7 years ago
tobytaogf • 0

Hi. I have some Phragmites RNA seq data. I tried to use Skewer to trim the adapter of these reads, as the code " skewer-master/skewer -m pe -l 50 Phragmites_RNA/raw_data/PL_R1.fastq Phragmites_RNA/raw_data/PL_R2.fastq -o Phragmites_RNA/raw_data/skewer/PL" After trimming, I run fastqc of these trimmed data to check the qualities. I found all these data have the kmer flag and position either on the front or the end position. However, the adaptor section did not give me the flag. I am wondering what cause these kmer issues and do I need to worry these for my downstream analysis. Or parameter setting of the skewer is wrong?

PL_R1 PM_R1 PL_R2

skewer fastqc kmer • 621 views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode
2.3 years ago

In your case, there is NO need to specify -m or -l option!

Adapter: -m, --mode <str> trimming mode; for paired-end files -- pe: paired-end; mp: mate-pair; ap: amplicon (DEFAULT: pe) Filtering: -l, --min <int> The minimum read length allowed after trimming; (DEFAULT: 18)

Please try the following command instead: skewer-master/skewer Phragmites_RNA/raw_data/PL_R1.fastq Phragmites_RNA/raw_data/PL_R2.fastq -o Phragmites_RNA/raw_data/skewer/PL

ADD COMMENT

Login before adding your answer.

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