trimming 151bp to shorter k-mer in ATAC-seq analysis
0
0
Entering edit mode
8 weeks ago
QX ▴ 60

Hi all,

I have dataset ATAC-seq for 151bp and process to the analysis. However, I read in some tutorial that "fragment length distribution of ATAC-seq libraries usually peaks at around 45–50 bp , and as a result, a majority of fragments are already fully covered by 2x36mer reads." For this reason, shall I cut down the read from 151 to a shorter length (like 65 or 35) before proceed the mapping?

ATAC-seq • 367 views
ADD COMMENT
0
Entering edit mode

Yes, read trimming is a normal step in ATAC-seq workflow, although I don't know that a majority of fragments are covered by 36mer reads.

You can try cutadapt, fastp, or another tool to trim reads with read-through.

ADD REPLY
0
Entering edit mode

yes I do trimming, but usually for adapters only. Do we need to trim to a specific number, for e.g. 36?

ADD REPLY
1
Entering edit mode

No, trim for adapters and call it a day.

ADD REPLY
0
Entering edit mode

so only trim for adapters and the length remain does not matter? do you know why they mention the term '2x36mer reads', does it make any biological sense of ATAC-seq?

ADD REPLY

Login before adding your answer.

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