Hello, I have a data of snatacSeq of 15 batch of patients (5 patients per batch that are multiplexed ), I need to demultplexe them using genotypes I am running demuxlet on my single-cell ATAC-seq (snATAC-seq) data, to demultiplexe my batchs
but the output file only contains the header line and no results.
In the log file , I see: NOTICE [2025/06/28 11:14:06] - Finished processing 0 droplets total NOTICE [2025/06/28 11:14:06] - Finished writing output files
and in the output.best i have only the header : BARCODE RD.TOTL RD.PASS RD.UNIQ N.SNP BEST SNG.1ST SNG.LLK1 SNG.2ND SNG.LLK2 SNG.LLK0 DBL.1ST DBL.2ND ALPHA LLK12 LLK1 LLK2 LLK10 LLK20 LLK00 PRB.DBL PRB.SNG1
My input BAM and VCF files seem to be correctly formatted.
There is some flags differences between RNA and ATAC bams but le CB is here, the GT is also present in my vcfs (I have one vcf bu batch)
Are there any specific requirements or known issues when using demuxlet with snATAC-seq data?
Could you advise on possible reasons why no droplets are being processed or written to the output?
Thank you for your help!
Looks like you posted this as an issue to the project GitHub site (which is probably your best option): https://github.com/statgen/demuxlet/issues/115
If you get an answer from the developer, please come back and post it here.