How to identify the miRNA counts column?

0

Entering edit mode

7.6 years ago

Shahzad ▴ 30

I have miRNA data of a plant. I aligned the reads using Bowtie (short read aligner) against miRbase database miRNAs. Then I used samtools -> bamtools -> bedtools to get the counts. Now I have counts imported into R but I am confused which column contains the counts. Please help to name all the variables in count file. enter image description here

Thank You.

counts miRNA NGS biostat bedtools • 2.3k views

ADD COMMENT • link 7.6 years ago by Shahzad ▴ 30

1

Entering edit mode

What program from bedtools did you use?

In general you would want to use featureCounts or HTSeq-count to get your read counts from aligned BAM files.

ADD REPLY • link 7.6 years ago by GenoMax 141k

0

Entering edit mode

I used "multicov" from bedtools.

ADD REPLY • link 7.6 years ago by Shahzad ▴ 30

0

Entering edit mode

I suggest that you use featureCounts instead. You can provide multiple bam files to featureCount. This will generate a matrix of counts that will have genes as rows and each sample as a new column.

ADD REPLY • link 7.6 years ago by GenoMax 141k

Login before adding your answer.