search

huangyh09 / brie

BRIE: Bayesian Regression for Isoform Estimate in Single Cells
https://brie.readthedocs.io
Apache License 2.0
41 stars 15 forks source link

Freeze during h5ad converting #65

Closed shiquan closed 1 month ago

shiquan commented 1 month ago

Hi,

I am processing data generated with Smartseq2. For some reason, I merged all cells into one sorted and indexed BAM file, with each cell distinguished by the RG tag. I used the following command to count the events, but the process always freezes at the last step. Below is the snapshot and the command I used:

brie-count -b cells.csv --cellTAG=RG --UMItag=RG -p 40 -o ./ -s merged.bam -a SE.filtered.gff3.gz

image

Thanks and looking forward to hearing from you soon.

BR, Quan

huangyh09 commented 1 month ago

Hi Quan, Thanks for sharing the issue. I'm unsure how to fix it quickly, but using UMItag as RG may cause another big problem. Since all reads in one cell have the same RG, it will treat them as one UMI if you use it UMItag. I guess you may have to split them into individual bam files, then input them into the brie-count.

Yuanhua

shiquan commented 1 month ago

Dear Yuanhua,

Thank you for your quick reply. After setting up different threads, I believe the issue may be caused by multithreading. The process freezes when I set -p to 40, but it works fine with other configurations. Additionally, I've noticed an issue with the UMItag, and I am currently working on splitting the cells. I appreciate your advice.

Quan