Hello everyone,
I am new to single cell RNA sequencing analysis and have been stuck on how to convert from .bam files to a gene-cell expression matrix. I am currently looking at data from the site: https://trace.ncbi.nlm.nih.gov/Traces/sra/?run=SRR10349606. I have downloaded the .bam file, and have been following the tutorial https://scrnaseq-course.cog.sanger.ac.uk/website/construction-of-expression-matrix.html. The tutorial talks about how to go from a fastq file to the gene-cell expression matrix, however it is my understanding that the .bam file I downloaded has already been aligned to a genome. I also looked into using featureCounts (bioconductor), but this required a cellbarcodeList, which I am unsure how to obtain.
Can someone point me in the right direction? I think I have a conceptual misunderstanding about this process.
Thank you!
0 answers
No answers yet.
Log in to answer this question.
-
Only one read per run - Trying to use CellRangerv7
written by Sky •
-
cell ranger low correct barcodes and failed
written by shersky •
-
Integrate CITESeq reads to Seurat object
written by Vinothkumar •
-
Ask for help! Need a tutorial for processing Drop-seq data and converting sequencing data into expr…
written by lyxwyc777 •
-
Filter out mouse reads from single cell RNA-seq of PDX tumor
written by fifty_fifty •
-
SRR data assembly
written by shamsurrahaman132 •
-
PacBio file type from NCBI run download
written by noodle
-
How to figure out which adapter has been used
written by juan.crescente •
-
Problem: "Pair-end" reads scRNA seq data (Drop-seq)
written by jayeshkumarsundaram •
-
How can I find SRA MD5 checksums for FASTQ files?
written by rivanov •
You can use featurecounts directly on the BAM file which is what the tutorial is doing. Then you can use umi-tools to count the output to generate the gene-cell expression matrix -if your data has UMIs-.