This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filter reads from fastq files (10x) with a specific cell barcode

Hello everyone,

I am analyzing perturb-seq data (10x). Unless cells were FACS sorted in the first hand, many of them (~40%) have not been associated to any sgRNA by cellranger. I observed that the BFP expression of the cells with no guide associated is comparable to other cells, showing that a CRISPR guide might be present but not detected / filtered out by cellranger ?

I would like to dig into this sgRNA attribution from the fastq files. To start, I would like to: 1) Select reads of a specific cell barcode (I chose a cell that has no guide according to cellranger, and a high depth) from the fastq file 2) Align these reads to a reference of all my sgRNAs (using bowtie2 ?) allowing 1 mismatch/indel

Do you have any recommandations for step 1 or 2 ? Would you suggest doing something else (I am very new to bioinformatics) ?

Best, Paul

fastq perturb-seq scrnaseq

1 answer

Log in to answer this question.