Hello!
I am trying to run some bam files though a pipeline that requires the bam files to contain the assignments of entire reads not read mates. I am using the Slide-seq2 data (from https://singlecell.broadinstitute.org/single_cell/study/SCP815/highly-sensitive-spatial-transcriptomics-at-near-cellular-resolution-with-slide-seqv2#study-summary ) which is not a single cell protocol and I am not sure what the bam files look like.
I would greatly appreciate some help to determine whether the bam files contain the assignments of the entire reads /entire alignment/ or read mates. And if they contain the read mates, how can I generate bam files containing the entire read alignments?
Thanks a lot!
•
written
by
Arman •
0 answers
No answers yet.
Log in to answer this question.
-
Tool: Loom - brings fragmented spatial transcriptomics workflows into one platform - need feedbacks
written by Olivia •
-
Forum: Cataloguing available spatial transcriptomics tools, then linking them together into a workf…
written by LauferVA
-
Decimal numbers in the third column of Gene/cell count matrix (mtx format).
written by sarapascu98 •
-
Single Cell Portal: Download troubles using curl
written by needHelpWkeggGSEA •
-
RNA velocity for smart-seq2 dataset
written by Dhanusha •
-
Extracting reads from BAM file based on partial read name
written by julianneradford •
-
combining / merging overlapping reads in bam file to one in silico long read
written by Lilly •
-
Extract all BAM reads that intersect a given region using the BAI index
written by bugmaster •
-
featureCounts: 0% successfully assigned fragments on multiple PE .BAM files
written by imesi •
-
Problem using rsem function for validate sam file
written by pier-luc.dudemaine •
It could have misunderstood your question, but if you really just want to see the number of reads mapped with their mates and the number of singletons, you can simply use
samtools flagstat.