hyunhwan-jeong / CB2

CB2 is an R package which provides functions for hit gene identification and quantification of sgRNA (single-guided RNA) abundances for CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) pooled screen data analysis. Details are in Jeong et al. (2019) <doi:10.1101/gr.245571.118> and Baggerly et al. (2003) <doi:10.1093/bioinformatics/btg173>.
https://cran.r-project.org/web/packages/CB2/index.html
Other
7 stars 1 forks source link

maximum guide length #15

Open LdeDiaz opened 1 year ago

LdeDiaz commented 1 year ago

I have had great success using CB2 for CRISPR studies and greatly appreciate your work on this package!

I have a project where I am looking to analyze enrichment of guide pairs, rather than individual guides. In another workflow, I simply merge the two guide sequences in my fasta reference files and also merge the two reads together from the fastq files, making a 42mer in each file that is quantified by a simple kmer match. I have so far been unsuccessful in using CB2 to quantify the guides when they are merged. I have have been able to confirm that the individual guides can be quantified properly by CB2 and that the kmer matching workflow I have used in the past is able to quantify the merged guides properly.

Can you please let me know what is the limit in guide length allowable to search with CB2?

Many thanks!

hyunhwan-jeong commented 1 year ago

Hi @LdeDiaz,

Thanks for reaching to us! The problem is caused due to CB2 only supports guides <= 30mer due to an algorithm limitation. I am interested in supporting the feature in CB2, it would be great if you can provide the following information:

Best regards,

Hyun-Hwan Jeong

hyunhwan-jeong commented 1 year ago

@LdeDiaz, I wonder if you are still interested in supporting the function. If so, could you provide the details regarding my questions?

Thank you,

Hyun-Hwan Jeong