是不是read 质量值编码有问题啊。你瞄下。
你好,我用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都是没有问题的。
希望您解答。