feature counts script output help please help
1
0
Entering edit mode
6 weeks ago

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 • 467 views
ADD COMMENT
0
Entering edit mode

what should the new files end in? .featurecounts?

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode
6 weeks ago
Mensur Dlakic ★ 20k

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.

ADD COMMENT

Login before adding your answer.

Traffic: 1208 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6