Question: Explanation for kmer enrichment after skewer trimming of the RNA-seq reads
0
gravatar for tobytaogf
9 months ago by
tobytaogf0
tobytaogf0 wrote:

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

fastqc kmer skewer • 254 views
ADD COMMENTlink modified 3 months ago by hongshan.jiang0 • written 9 months ago by tobytaogf0

How to add images to a Biostars post

ADD REPLYlink written 9 months ago by genomax74k
0
gravatar for hongshan.jiang
3 months ago by
hongshan.jiang0 wrote:

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 COMMENTlink written 3 months ago by hongshan.jiang0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1258 users visited in the last hour