Preprocessing#
Import Pycallingcards as:
import pycallingcards as cc
Clean qbed with some extra chromosomes. |
|
Call peaks from qbed data. |
|
Combine two peaks. |
|
Separate two peaks. |
|
Annotate the peak data using bedtools [Quinlan and Hall, 2010]. |
|
Make cell(sample) by using peak anndata for calling cards. |
|
Filter peaks based on the number of cells or counts. |
|
Calculate sum of peaks per cluster or average peaks per cell in different cluster and give it to anndata object. |
|
Down sample insertion data. |