Is there a tool for calculating the number of uniquely mapped reads, ambigiously mapped reads, and unmapped reads from sam/bam?
I found a tool called MMQuant but this is suppose to be drop in replacement (though, it's not) for featureCounts. I'm trying to calculate the number of uniquely mapped reads, the number of ambigiously mapped reads, and the number of unmapped reads.
Is there a tool that can calculate these basic measures for sam/bam files?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
featureCounts percentage of successful assigned alignments lower than expected
written by gt •
-
trouble with 'too short' reads using STAR. each mate maps fine, but the paired end mapping does not
written by james.alan.gregory •
-
FeatureCounts results: does it look correct? Why so many unassigned unmapped reads?
written by tanya_fiskur •
-
Difference between read counts STAR vs featureCounts
written by sg197 •
-
Is there a preferred RNA-seq aligner for (very) short reads and potential homologous regions?
written by O.rka
-
Realign Aligned.sortedByCoord.out.bam to Aligned.toTranscriptome.out.bam for RSEM input
written by mkvalluru •
-
multi-mapped reads counting methods with RNA-seq
written by guillaume.rbt
-
Questions about STAR and Bowtie (not Bowtie2) mapping
written by debitboro
-
Problems with STAR mapping output
written by manekineko
-
Should Samtools Pileup Be Performed On Uniquely Mapped Reads Or All The Reads?
written by Doctoroots
why not just use featureCounts in the first place? (== might be useful to mention why do you want an alternative for it when asking for alternatives)
htseq-count ?
I’m getting a large number of ambiguously mapped reads and a very inconsistent gene expression between replicates. I suspect it’s a result of these ambiguous reads. I want to know exactly what proportion of the read pairs are ambiguously mapped.
I hadn't realized this information is in the summary output of featureCounts.
use combination of samtools and broad institute's website for flags
BBTools (
reformat.shhas a number of stats options) for aligned data files.