Question: Peak call? or Just convert to bedfile by bedtools bamtobed? to generate the accessibility matrix of single-end single cell ATAC-seq
0
gravatar for soya
5 days ago by
soya0
soya0 wrote:

I'm wondering about the "single-end" single cell ATAC-seq (scATAC-seq) analysis. The reads obtained from each cell is at most 1 read at each region, by nature.

Some paired-end scATAC-seq pipelines execute the peak-calling by MACS2 etc., but it seems to be hard to detect peaks from single-end single cell data. Actually, the sum of obtained peak region was about only 3% of the read coverage simply calculated by bedtools bamtobed when I tried with my sample data.

Then, how should I treat the "open chromatin region" of each single cell?

Calculate the peak region with the aggregate of single cell, and then consider whether there is a read in each cell by bedtools intersect is the best?

single cell atac-seq • 56 views
ADD COMMENTlink written 5 days ago by soya0
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: 1415 users visited in the last hour