I love this forum very much,everyone here is very kind.
I have run the MACS on my bacterial whose genome is about 26mb
and there are about 10000000 reads mapped to the genome.
However, the out file of MACS contains only 53 peak regions,much less than the 3000 predicted gene.
Maybe there are some problem about my outcome.
my command is
macs14 -t testsorted.bam -f BAM -g 3.0e6 -n test -w
Could you give me some advice?
(my data is rna-seq, I don't have the control file, so I did't use the -c option,is it ok?)
ps.
I also countered this odd problem.
when I run
macs14 -t testsorted.bam -f BAM -g 3.0e6 -n test -w
in the bash,it is ok.
but,when I write the same command in the .sh file,and qsub it.
I will get this error.
/usr/local/bin/python: bad interpreter: No such file or directory
I am not familiar with python,could anyone tell me how to deal with this problem? I have used qsub successfully for many times before today.
I have run the MACS on my bacterial whose genome is about 26mb
and there are about 10000000 reads mapped to the genome.
However, the out file of MACS contains only 53 peak regions,much less than the 3000 predicted gene.
Maybe there are some problem about my outcome.
my command is
macs14 -t testsorted.bam -f BAM -g 3.0e6 -n test -w
Could you give me some advice?
(my data is rna-seq, I don't have the control file, so I did't use the -c option,is it ok?)
ps.
I also countered this odd problem.
when I run
macs14 -t testsorted.bam -f BAM -g 3.0e6 -n test -w
in the bash,it is ok.
but,when I write the same command in the .sh file,and qsub it.
I will get this error.
/usr/local/bin/python: bad interpreter: No such file or directory
I am not familiar with python,could anyone tell me how to deal with this problem? I have used qsub successfully for many times before today.
Comment