max length in cutadapt
1
0
Entering edit mode
4.0 years ago
max_19 ▴ 170

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,

RNA-Seq sequencing adapters • 1.7k views
ADD COMMENT
2
Entering edit mode
4.0 years ago

"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).

ADD COMMENT
0
Entering edit mode

That makes sense, thank you. I think the part that confuses me is the "Reads that are too long even before adapter removal are also discarded."

ADD REPLY
1
Entering edit mode

The documentation is wrong there.

ADD REPLY

Login before adding your answer.

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