Analyzing data produced by the UMI based SCRB-seq (Soumillon et al. 2014) protocol.
Hello All,
I am working on analyzing single cell RNA-seq data generated using the SCRB-seq protocol. The data that I have has already gone through data pre-processing including demultiplexing the BCL files using Picard, alignment with STAR, QC with RNASeqQC and marking duplicates with UmiAwareMarkDuplicatesWithMateCigar also with Picard. The only sequence files that I have to access to are the BAM files with duplicates marked using Picard. I would like to run RSEM or HT-seq for transcript quantification. Is it a good idea to just remove the Picard marked duplicates and then run RSEM on the data? I appreciate all help and suggestions.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Biobambam vs Picard in RNA Sequencing Analysis for BAM to FASTQ
written by kcarey •
-
Genecount-difference between HT-seq count, RSEM, and Kallisto
written by tvibhaps •
-
How to check for Over or Under-marking Duplicates in BAM files
written by noob@bioinfo •
-
Removing umi split off as a separate fastq (RNA-seq)
written by YJ •
-
RSEM, interpreting quantification results
written by argonvibio •
-
Which genome index should be used for "rsem-calculate-expression" command after RNA-seq reads mapp…
written by taoi2 •
-
Generating a gene expression matrix from Picard *duplicates_marked.bam file
written by halo22
-
Unable to create star refernce file with rsem
written by Debbie •
-
picard tools HsMetrics not using duplicate reads to calculate coverage
written by Floydian_slip
-
Somatic Sniper Producing Slightly Different Call Numbers Depending On Whether Duplicates Are Marked…
written by Travis
