This is a test version of Biostars. For the public version, visit https://www.biostars.org.
feature counts script output help please help

Where are my feature count files going when I use this script:

featureCounts -a ncbi_dataset/data/GCA_002738495.1/genomic.gff \
-o Feature_counts.txt \
-p -s 2 -T 10 \
*_sorted.bam

It's running but can't find the file...no feature_counts.txt folder or cd?

RNAseq_CD_1_CCGCGGTT_Arabidopsis_I1187_L1_sorted.bam

Are they in these folders? I thought these folders were just the sam to bam conversion ones.

feature counts output script
• 2,354 views • link
updated by Mensur Dlakic Mensur Dlakic • written by smithkthedale • smithkthedale

what should the new files end in? .featurecounts?

• 0 views • link
written by smithkthedale • smithkthedale

I think featureCounts prints a log to standard output, which if you provide to use can help us figure out the problem.

• 0 views • link
written by rpolicastro rpolicastro

how do I use verbose in my script. when I type in verbose it says command not found

• 0 views • link
written by smithkthedale • smithkthedale

It says so already right below your question:

Maybe add the --verbose switch to your command and tell us what was the screen output.

That means adding --verbose to your featureCounts command anywhere before you type *_sorted.bam.

• 0 views • link
written by Mensur Dlakic Mensur Dlakic

1 answer

If you type featureCounts without any arguments, it will give help:

  -o <string>         Name of output file including read counts. A separate file
                      including summary statistics of counting results is also
                      included in the output ('<string>.summary'). Both files
                      are in tab delimited format.

Based on your commands, it should be saved in Feature_counts.txt. Maybe add the --verbose switch to your command and tell us what was the screen output.

• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Log in to answer this question.

More posts like this