Dear all,
I have been using tophat for my RNA-seq data. The command line I'm using is the standard: ./tophat -p 8 -G /path_to_genes.gtf/ -o <output_folder> /path_to_genome_index/ XXX.fastq
Tophat runs happily, finds bowtie and samtools, checks for bowtie index, Fasta reference, generates SAM header... all goes well and I have 4 outputs: accepted_hits.bam deletions.bam insertions.bam and junctions.bam... However, when the sam version of accepted_hits.bam is going to be written, samtools fails and I get a message saying sh:samtools: command not found...
Having gone through with the bam files till the very end (cufflinks then cuffdiff and R), I seem to get the analysis running but I don't understand why samtools fails during the mapping.... Tophat clearly recognises it at first.. and ./samtools <options> work nicely for other applications.
Just to add I have used samtools-0.1.17.tar.bz2 (I'm on a iMac) untar ed it and cp ed it into my working path (I haven't got permission for /usr/bin), I have also cp ed libbam.a into my working directory...
Can anyone help me
Cheers
lj907
I have been using tophat for my RNA-seq data. The command line I'm using is the standard: ./tophat -p 8 -G /path_to_genes.gtf/ -o <output_folder> /path_to_genome_index/ XXX.fastq
Tophat runs happily, finds bowtie and samtools, checks for bowtie index, Fasta reference, generates SAM header... all goes well and I have 4 outputs: accepted_hits.bam deletions.bam insertions.bam and junctions.bam... However, when the sam version of accepted_hits.bam is going to be written, samtools fails and I get a message saying sh:samtools: command not found...
Having gone through with the bam files till the very end (cufflinks then cuffdiff and R), I seem to get the analysis running but I don't understand why samtools fails during the mapping.... Tophat clearly recognises it at first.. and ./samtools <options> work nicely for other applications.
Just to add I have used samtools-0.1.17.tar.bz2 (I'm on a iMac) untar ed it and cp ed it into my working path (I haven't got permission for /usr/bin), I have also cp ed libbam.a into my working directory...
Can anyone help me
Cheers
lj907
Comment