cortes-ciriano-lab/SComatic

Why is the default value of the parameter "min_mq" 255

Topigs opened this issue · 1 comments

Topigs commented

Thanks so much for such a wonderful tool!
In view the script "scripts/BaseCellCounter/BaseCellCounter py", I am very confused with the default value of the command line option "-- min_mq". Why set it to 255 instead of something more reasonable, I mean like 30.

Dear user,
Thanks for using SComatic. This default mapping quality is ideal for detecting uniquely mapping reads in scRNA-seq data (usually aligned with STAR).

-- min_mq 30 would be a better parameter for DNA-based approaches like scATAC-seq.

Thanks,
Fran