tophat reporting output problem
0
0
Entering edit mode
6.4 years ago
pavanraaz • 0

Hi

I new to rna seq analysis, i have reference genome and tried to map reads. After running Tophat i got below error. what might be the reason and how to solve it. I need your valuable suggestions.

 Joining segment hits
 Reporting output tracks
    [FAILED]
Error running /usr/local/bin/tophat_reports --min-anchor 8 --splice-mismatches 0 --min-report-intron 50 --max-report-intron 500000 --min-isoform-fraction 0.15 --output-dir ./tophat_out/ --max-multihits 20 --max-seg-multihits 40 --segment-length 25 --segment-mismatches 2 --min-closure-exon 100 --min-closure-intron 50 --max-closure-intron 5000 --min-coverage-intron 50 --max-coverage-intron 20000 --min-segment-intron 50 --max-segment-intron 500000 --read-mismatches 2 --read-gap-length 2 --read-edit-dist 2 --read-realign-edit-dist 3 --max-insertion-length 3 --max-deletion-length 3 -z gzip --gtf-annotations /home/pavan/cajan_bowtie2-index/cajan.gff --gtf-juncs ./tophat_out/tmp/cajan.juncs --no-closure-search --no-coverage-search --no-microexon-search --sam-header ./tophat_out/tmp/cajan_genome.bwt.samheader.sam --report-discordant-pair-alignments --report-mixed-alignments --samtools=/usr/local/bin/samtools_0.1.18 --bowtie2-max-penalty 6 --bowtie2-min-penalty 2 --bowtie2-penalty-for-N 1 --bowtie2-read-gap-open 5 --bowtie2-read-gap-cont 3 --bowtie2-ref-gap-open 5 --bowtie2-ref-gap-cont 3 /home/pavan/cajan_bowtie2-index/cajan.fa ./tophat_out/junctions.bed ./tophat_out/insertions.bed ./tophat_out/deletions.bed ./tophat_out/fusions.out ./tophat_out/tmp/accepted_hits ./tophat_out/tmp/left_kept_reads.m2g.bam,./tophat_out/tmp/left_kept_reads.m2g_um.mapped.bam,./tophat_out/tmp/left_kept_reads.m2g_um.candidates.bam ./tophat_out/tmp/left_kept_reads.bam
Loaded 125959 GFF junctions from ./tophat_out/tmp/cajan.juncs.

 

RNA-Seq tophat software error • 2.8k views
ADD COMMENT
0
Entering edit mode

Try running that command manually so you can see the actual error message. For what it's worth, this seems to be a common place for tophat2 to die.

ADD REPLY

Login before adding your answer.

Traffic: 2201 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