Question: How to extract a specific gene sequence for a specific cell from scRNA-seq fastq file
gravatar for cook.675
6 months ago by
cook.67530 wrote:

I have a scRNA seq experiment where I would like to look at the raw sequencing data for one specific gene.

We have 3' ends, and the sequence may or may not be different depending on the cell type we are looking at.

I have the list (barcodes) from the cells we are interested in looking at, to see what the sequence is, and if its different for the different cell types.

How would I find and extract this information from the fastq files? We have paired end sequencing

Thanks in advance!

sequencing scrna-seq • 245 views
ADD COMMENTlink modified 6 months ago by mahmoud.s.fahmy0 • written 6 months ago by cook.67530
gravatar for mahmoud.s.fahmy
6 months ago by
mahmoud.s.fahmy0 wrote:

A similar post, here.

You would have to first align the raw reads to the reference genome.

ADD COMMENTlink written 6 months ago by mahmoud.s.fahmy0

Thanks for that link, one quick question:

All the reads we have are transcript reads, I could then map to the transcriptome instead of the genome?

Im going through this now, I created the SAM file from paired end fastq's following this tutorial:

And am now converting them to BAM, and I will have to search for the next step

ADD REPLYlink written 6 months ago by cook.67530
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1226 users visited in the last hour