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?
•
written
by
tobytaogf •
1 answer
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
Log in to answer this question.
-
Mimicing the Blastn "lineage report" from online Blast service
written by mschmid
-
Parsing RNAhybrid output
written by baurumon •
-
Bam QC, GC content
written by Martina •
-
More informative GSEA enrichment ggplot
written by stephenwongscience •
-
RNA-seq analysis using RUVSeq: am I doing it right?
written by Ankur B. Sharma •
-
Blastn .XML to FASTA
written by baurumon •
-
How to create PSSM matrix
written by sepehr.salsabili •
-
Is this common to see two types of adapters in the same sequencing data?
written by lishen0709 •
-
Forward and Reverse strand on gene annotation
written by DanielC
-
DESeq2 log2FoldChange vs Salmon log2 TPM
written by liartom2 •
How to add images to a Biostars post