Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
how many reads per sample would you recommend in single cell RNA-Seq (e.g. Cel-Seq) super0925 RNA Sequencing 0 09-06-2016 02:02 PM
Problem with HTSeq in RIPseq pipeline Rivalyn Bioinformatics 7 04-17-2014 06:22 AM
Problem with Resuming TopHat pipeline with unmapped reads bjoernoest Bioinformatics 13 08-13-2013 06:35 AM
Corona Lite 4.0 Pairing Pipeline problem nedoluzhko Bioinformatics 7 09-30-2009 11:19 AM

Thread Tools
Old 01-17-2017, 10:27 AM   #1
Junior Member
Location: Long Island

Join Date: Nov 2014
Posts: 4
Default Problem with CEL-Seq-pipeline

I wonder anyone is using the CEL-Seq pipeline. I am new to CEL-Seq as well as to bioinformatics. I got stuck at the very first step. Can anyone help?

Many thanks!

Below is the error message:

2017-01-17 14:04:09,541 - pijp - INFO - ===== Started pijpleiding with config file : ./config_file.txt =====
2017-01-17 14:04:09,557 - pijp - INFO - ===== pipeline at /sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq-pipeline-master/ =====
2017-01-17 14:04:09,557 - pijp - INFO - ========== writing log to file at /sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq_t1/pijp.log ==========
2017-01-17 14:04:09,558 - pijp - INFO - Section : bc_demultiplex
2017-01-17 14:04:09,558 - pijp - INFO - parameters : {"strand_field": "2", "sample_sheet": "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/sample_sheet", "bc_index_file": "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/bc_index_file", "umi_length": "0", "bc_length": "8", "min_bc_quality": "10", "output_dir": "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq_t1", "lane_field": "1", "index_field": "0", "input_files": ["/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI1_L001_R1_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI1_L001_R2_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI1_L002_R1_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI1_L002_R2_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI3_L001_R1_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI3_L001_R2_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI3_L002_R1_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI3_L002_R2_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI7_L001_R1_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI7_L001_R2_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI7_L002_R1_001.fastq", "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/T1/RPI7_L002_R2_001.fastq"], "stats_file": "", "fname_delimiter": "_", "cut_length": "35"}
Traceback (most recent call last):
File "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq-pipeline-master/", line 129, in <module>
File "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq-pipeline-master/", line 93, in main
File "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq-pipeline-master/", line 37, in main
fastq_dict = create_fastq_dict(input_files, fname_delimiter, int(lane_field), int(index_field), int(strand_field))
File "/sonas-hs/dubnau/hpc_norepl/data/mshih/CEL-Seq/CEL-Seq-pipeline-master/", line 115, in create_fastq_dict
assert ("R1" in fastq_dict[fastq_file].strand), "File name does not contain R1 in field %d. Aborting (%r not R1)" %(strand_field, fastq_dict[fastq_file].strand)
AssertionError: File name does not contain R1 in field 2. Aborting ('R2' not R1)
MaxwellShih is offline   Reply With Quote

cel-seq, python, rna-seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:15 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO