Monday, November 13, 2023

Notes on snapatac2

When I ran:

data = snap.pp.import_data(

    fragment_file,

    chrom_sizes=snap.genome.hg38,

    #file="2_filtered_data/ATAC/atac_test.h5ad",  # Optional

    sorted_by_barcode=False,)

------------------------

fragment_file being fragments.tsv.gz and fragments.tsv generated the error

PanicException: called `Result::unwrap()` on an `Err` value: Custom { kind: Other, error: "Unknown frame descriptor" }

instead i used fragments.srt.bed.gz generated from :

time sort -k4,4 fragments.tsv | gzip - > fragments.srt.bed.gz,

resolved the issue

No comments:

Post a Comment