Hi all! I'm trying to remove adapter sequences from my RNA-seq dataset, and i am only interested in reads that are 20-30bp in length. So I used cutadapt with the following specifications:
cutadapt -a $adapter -m 20 -M 30
In the cutadapt user guide, it states that
--maximum-length LENGTH or -M LENGTH Discard processed reads that are longer than LENGTH. Reads that are too long even before adapter removal are also discarded.
my raw reads are all 50bp, so they are already longer than my specified M value (30). Does this mean that cutadapt will get rid of all my reads?
When I run cutadapt though, I still get a good # of reads left over (e.g. 80-90%), so i'm not sure if I understand this M parameter incorrectly.
Thanks for the input guys,
•
written
by
max_19
1 answer
"processed reads" are what's left after trimming but before filtering for length. Since you have a sRNA prep most of the reads will get trimmed and likely be below 30 bases thereafter (thus being under the maximum length).
The documentation is wrong there.
Log in to answer this question.
-
Trimmomatic or Cutadapt not Trimming my Adapters even though I see them clearly in the Fastq-File
written by obsto123 •
-
cutadapt not trimming all adapters in PE Illumina sequencing
written by noodle
-
Error in sequence length distribution
written by bioinformatics.queries
-
Variable read length distribution after cutadapt running for my ATAC-seq datasets
written by yaoyao20152031 •
-
STAR or Bowtie for small RNA seq?
written by max_19
-
Cutadapt wrong behaviour
written by dllopezr •
-
Content of T was larger than others in NSR-RNAseq data, why?
written by wm
-
cutadapt for trimming reads to specific length
written by Varun Gupta
-
Many 0 length reads after trimming
written by marcon •
-
How To Delete The All Fastq Reads Which Includes A Potential 50Bp Illumina Single End Pcr Primer 1
written by Tonyzeng •
