5 bowtie跑colorspace出现的错误

你好,我用bowtie跑colorspace时出现了如下错误:

[2017-08-09 21:35:34] Beginning TopHat run (v2.1.1)
-----------------------------------------------
[2017-08-09 21:35:34] Checking for Bowtie
          Bowtie version:     1.2.1.0
[2017-08-09 21:35:34] Checking for Bowtie index files (genome)..
[2017-08-09 21:35:34] Checking for reference FASTA file
[2017-08-09 21:35:34] Generating SAM header for /home/chenyijun/bowtie_index/hg19/hg19_c
[2017-08-09 21:35:34] Reading known junctions from GTF file
[2017-08-09 21:35:36] Preparing reads
     left reads: min. length=50, max. length=50, 11562869 kept reads (5895 discarded)
[2017-08-09 21:37:53] Building transcriptome data files P5_1/tmp/hg19
[2017-08-09 21:37:58] Building Bowtie index from hg19.fa
[2017-08-09 21:43:18] Mapping left_kept_reads to transcriptome hg19 with Bowtie
    [FAILED]
Error running bowtie:
Reads file contained a pattern with more than 1024 quality values.
Please truncate reads and quality values and and re-run Bowtie
terminate called after throwing an instance of 'int'

运行程序的命令如下:

tophat -G ${gtf} --bowtie1 -C -Q --output-dir P5_1 --no-novel-juncs --no-coverage-search --num-threads 10 ${BowtieIndex} SRR1020087_1.csfasta,SRR1020088_1.csfasta,SRR1020089_1.csfasta,SRR1020090_1.csfasta,SRR1020091_1.csfasta,SRR1020092_1.csfasta SRR1020087_1.qual,SRR1020088_1.qual,SRR1020089_1.qual,SRR1020090_1.qual,SRR1020091_1.qual,SRR1020092_1.qual > P5_1.log 2>&1

而且tophat和bowtie都是没有问题的。

希望您解答。


请先 登录 后评论

1 个回答

张海伦 - 生物信息分析员

是不是read 质量值编码有问题啊。你瞄下。

请先 登录 后评论
  • 1 关注
  • 0 收藏,3935 浏览
  • 陈东君 提出于 2017-08-09 21:51

相似问题