Question: using featureCounts in command line
0
gravatar for A
5.1 years ago by
A3.9k
A3.9k wrote:

hi dudes,

i want to count the number of my reads that have been mapped to each gene, then i downloaded featureCounts from http://sourceforge.net/projects/subread/files/subread-1.4.6-p5/ , untared that then using my accepted_hits.sam and genes.gtf as input, i replaced the command from documentation like below:

Summarize a single-end read dataset using 5 threads:
featureCounts -T 5 -t exon -g gene_id -a annotation.gtf -o counts.txt mapping_results_SE.sam

[izadi@lbox161 Downloads]$ ls
subread-1.4.6-p5-source  subread-1.4.6-p5-source.tar.gz
[izadi@lbox161 Downloads]$ cd subread-1.4.6-p5-source
[izadi@lbox161 subread-1.4.6-p5-source]$ ls
LICENSE  README.txt  annotation  doc  ouput  src  test
[izadi@lbox161 subread-1.4.6-p5-source]$ cd src
[izadi@lbox161 src]$ ln -s /usr/data/nfs6/izadi/angel/ouput/
[izadi@lbox161 src]$ featureCounts -T 5 -t exon -g gene_id -a genes.gtf -o counts.txt accepted_hits.sam 
bash: featureCounts: command not found...
[izadi@lbox161 src]$ 

after untaring there no bin file there? do you know why?

[izadi@lbox161 izadi]$ tar -zxvf subread-1.4.6-p5-source.tar.gz

 

may you please read my command and detect my fault here?

thank you

 

myposts featurecounts • 9.6k views
ADD COMMENTlink modified 5.1 years ago • written 5.1 years ago by A3.9k
5

please, again, learn from your previous experience :

* no happen when runing tophat in fedora

* the number of reads that have mapped to the each gene

I'm going to close that one.

ADD REPLYlink modified 5.1 years ago • written 5.1 years ago by Pierre Lindenbaum131k

ohhhh Pierre please,

im really helpless and this website is the only where i can get help because im totally alone with my problem

ADD REPLYlink written 5.1 years ago by A3.9k
6

"im totally alone with my problem"

Just you, man command + Google and manuals. Start with googling:

  1. compiling a program from a source
  2. mastering ls and find command. Because you know, the secret incantation featureCounts can show its magic power only if it actually exist on your system and you invoke it with another short prayer which we call /path/to/your/executable or (but this is a real Black Magic) include some extra info in a book containing good spells, some like to call $PATH. Something like this (bash):
export PATH=/path/to/directory/where/your/executable/is:$PATH

Trust me, these are good spells to know. Without them not even basic magic tricks are possible. And nobody will love you in Hogwarts.

Edit: removed space after "PATH= " + formatted as code

ADD REPLYlink modified 5.1 years ago • written 5.1 years ago by Darked894.2k

thank you so much for your caring...yes you all right

ADD REPLYlink written 5.1 years ago by A3.9k
3
gravatar for A
5.1 years ago by
A3.9k
A3.9k wrote:


Enter the src subdirectory under the home directory of the package and then issue the following command to build it on a Linux/unix computer:
make -f Makefile.Linux

then the bin will be created and Enter that subdirectory under the home directory, now run your command

as simple as eating a piece of cake but not cheese one!

ADD COMMENTlink written 5.1 years ago by A3.9k
0
gravatar for EagleEye
5.1 years ago by
EagleEye6.7k
Sweden
EagleEye6.7k wrote:
When you are in 'src' folder ofsubread program. Try refering your current directory. Example run like this ./featureCounts -T 5 -t exon -g gene_id -a genes.gtf -o counts.txt accepted_hits.sam
ADD COMMENTlink written 5.1 years ago by EagleEye6.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1098 users visited in the last hour