Hi all
I tried Tophat on 35 bp Illumina reads and 50bp Solid reads. Failed by giving warning Junction database is empty. I posted this error message with logs but no response. As I see few others also posted the same error but the responses given were not so useful.
Could any one please assist me regarding finding new splice junctions by using 35 bp or 50 bp RNASeq libraries (Illumina or Solid)
TopHat-30bp RNASeq read library Error
TopHat-50bp RNASeq read library Error
Then I shifted to Splicemap. It is also failed by giving error message
SpliceMap error: same for both 35 bp and 50bp Libraries
I tried Tophat on 35 bp Illumina reads and 50bp Solid reads. Failed by giving warning Junction database is empty. I posted this error message with logs but no response. As I see few others also posted the same error but the responses given were not so useful.
Could any one please assist me regarding finding new splice junctions by using 35 bp or 50 bp RNASeq libraries (Illumina or Solid)
TopHat-30bp RNASeq read library Error
Code:
Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Warning: found a read < 20bp in SL001_R00089_RHE011_01pgx2_F3.fasta Traceback (most recent call last): File "/home/bogugk/XXX/tophat-1.0.13/bin/tophat", line 1635, in ? sys.exit(main()) File "/home/bogugk/XXX/tophat-1.0.13/bin/tophat", line 1562, in main params.read_params = check_reads(params.read_params, left_reads_list) File "/home/bogugk/XXX/tophat-1.0.13/bin/tophat", line 719, in check_reads if line_num % 2 == 1: KeyboardInterrupt
Code:
]$ tophat -o tophat_out_SRR017234_1.fas tq /home/bogugk/rnaseq/software/bowtie-0.12.5/indexes/hg18_c SRR017234_1.fastq [Wed Aug 4 20:40:52 2010] Beginning TopHat run (v1.0.13) ----------------------------------------------- [Wed Aug 4 20:40:52 2010] Preparing output location tophat_out_SRR017234_1.fastq/ [Wed Aug 4 20:40:52 2010] Checking for Bowtie index files [Wed Aug 4 20:40:52 2010] Checking for reference FASTA file [Wed Aug 4 20:40:52 2010] Checking for Bowtie Bowtie version: 0.12.5.0 [Wed Aug 4 20:40:52 2010] Checking reads seed length: 35bp format: fastq quality scale: phred33 (default) [Wed Aug 4 20:42:21 2010] Mapping reads against hg18_c with Bowtie [Wed Aug 4 20:42:21 2010] Joining segment hits [Wed Aug 4 20:42:21 2010] Searching for junctions via segment mapping Warning: junction database is empty! [Wed Aug 4 20:44:06 2010] Joining segment hits [Wed Aug 4 20:44:06 2010] Reporting output tracks [FAILED] Error: Report generation failed with err = 1 Traceback (most recent call last): File "/home/bogugk/XXX/tophat-1.0.13/bin/tophat", line 1635, in ? sys.exit(main()) File "/home/bogugk/XXX/tophat-1.0.13/bin/tophat", line 1607, in main params.gff_annotation) File "/home/bogugk/XXX/tophat-1.0.13/bin/tophat", line 1044, in compile_reports exit(1) TypeError: 'str' object is not callable
SpliceMap error: same for both 35 bp and 50bp Libraries
Code:
k@Trivia SpliceMap3313_linux-64]$ ./bin/runSpliceMap run.cfg ---== Welcome to SpliceMap 3.3.1.3 ==--- Developed by Kin Fai Au and John C. Mu http://www.stanford.edu/group/wonglab/SpliceMap/ __________ Loading configuration file... Genome Directory: ../SpliceMap3313_linux-64/genome/ Number of mismatch allowed: 2 Mapping option: B Minimum (25th-percentile) intron size: 20000 Maximum (99th-percentile) intron size: 400000 Ref path: name: all.gene.refFlat.txt package path: ./bin/ name: runSpliceMap Read format: FASTA Number of threads: 2 Separating non-unique coverage Will print Cufflinks SAM file COMMAND: mkdir temp mkdir: cannot create directory `temp': File exists List 1: /home/bogugk/rnaseq/shiyan_data/csfasta_Folder/SL001_R00089_RHE011_01pgx2_F3.fasta List 2: Preparing the reads!... Using bases 1 to 1 (inclusive) of each read for mapping. 1 bases in total. I'm sorry, SpliceMap 3.3.1.3 only supports read lengths >= 50 Please contact: http://www.stanford.edu/group/wonglab/SpliceMap/ if you need this requirement
Comment