This is a test version of Biostars. For the public version, visit https://www.biostars.org.
program repair in featureCounts

Hello,

I try to use featureCounts to summarize the human pair-end read.

The manual of Subread said "We also provide an utility program repair to allow users to pair up the reads before feeding them to featureCounts".

Where is the program "repair"?

Thanks in advance!

rna-seq
• 1,828 views • link
updated by Biostar Biostar • written by yueli7 yueli7

How did you map your reads? Most likely you don't need to repair the bam file.

How did you install Subread? The repair program is under $SUBREAD_FOLDER/bin/utilities/.

• 0 views • link
written by h.mon h.mon

Thanks!

I ran STAR to map my reads.

The output files are: 2-pass_Aligned.out.bam, 2-pass_Aligned.sortedByCoord.out.bam, 2-pass_Aligned.toTranscriptome.out.bam.

I just want to double check: I will use 2-pass_Aligned.sotredByCoord.out.bam to do the following work, it that correct?

Thanks!

• 0 views • link
written by yueli7 yueli7

1 answer

Indeed, you don't need to repair your reads, and yes, use the "2-pass_Aligned.sotredByCoord.out.bam" with featureCounts. By the way, you could have used STAR --quantMode GeneCounts to get a table of counts similar to that of featureCounts.

• 0 views • link
written by h.mon h.mon

Thanks for your great help!

• 0 views • link
written by yueli7 yueli7

The reaction from h.mon got moved to an answer, so can you accept it if it solves your question?

• 0 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this