How does featurecounts deal with the output from Bowtie?
Hi there,
If the following command line is used to generate the SAM file.
bowtie -q -v 0 -k 10 -S -t <index.name> <input.fastq> <out.sam>
featureCounts program will be used to count the reads by enabling the multi-mapping reads counting:
featureCounts -t miRNA -g Name -O -M -a <mirbasefile> -o <outfile> <samfiles>
I wonder how featureCounts with the tags such as -O and -M will hand this situation as we know that this SAM file doesn't have "NH" tag? Will the command line above work properly for the microRNA quantification?
Many thanks,
Tom
1 answer
Log in to answer this question.
More posts like this
-
low Successfully assigned alignments
written by Srujana •
-
Exon counts do not match gene and transcript counts
written by firestar
-
featureCounts
written by Hien •
-
FeatureCounts for miRNAseq
written by Muc09 •
-
The number of peaks in microRNA-Seq trimmed data
written by wangdp123
-
FeatureCounts count reads per gene (Including Introns)
written by paulranum11
-
FeatureCounts summarizing at exon level
written by yamar •
-
Extremely low assigned reads in featureCounts
written by xqyjxau •
-
why the featurecounts result is empty
written by zizigolu
-
htseq-count problem with --samout option
written by Ian
