Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
R
rnaseq
  • Project overview
    • Project overview
    • Details
    • Activity
    • Releases
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 9
    • Issues 9
    • List
    • Boards
    • Labels
    • Service Desk
    • Milestones
    • Iterations
  • Merge Requests 0
    • Merge Requests 0
  • Requirements
    • Requirements
    • List
  • CI / CD
    • CI / CD
    • Pipelines
    • Jobs
    • Schedules
    • Test Cases
  • Security & Compliance
    • Security & Compliance
    • Dependency List
    • License Compliance
  • Operations
    • Operations
    • Incidents
    • Environments
  • Packages & Registries
    • Packages & Registries
    • Container Registry
  • Analytics
    • Analytics
    • CI / CD
    • Code Review
    • Insights
    • Issue
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Members
    • Members
  • Collapse sidebar
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
  • BICF
  • Astrocyte
  • rnaseq
  • Issues
  • #13

Closed
Open
Opened May 02, 2019 by Jeon Lee@s167719Developer

An error in 'trimgalore.sh'

Hi, An user, skim17, reported an error generated during the 'trim' step such as [Process trim (3) terminated with an error exit status (1)]. After looking at the relevant scripts, we suspect that the error came from some lines in 'trimgalore.sh' below:

r1base="${fq1%.fastq*}" r2base="${fq2%.fastq*}" source /etc/profile.d/modules.sh module load trimgalore/0.4.1 cutadapt/1.9.1

if [ -s $fq2 ] then trim_galore --paired -q 25 --illumina --gzip --length 35 ${fq1} ${fq2} mv ${r1base}_val_1.fq.gz ${pair_id}.trim.R1.fastq.gz mv ${r2base}_val_2.fq.gz ${pair_id}.trim.R2.fastq.gz else trim_galore -q 25 --illumina --gzip --length 35 ${fq1} mv ${r1base}_trimmed.fq.gz ${pair_id}.trim.R1.fastq.gz cp ${pair_id}.trim.R1.fastq.gz ${pair_id}.trim.R2.fastq.gz fi

because the user's fastq file names are in a format of 'SampleID_ReplicateNumber_ReadNumber.fq.gz', e.g. 'ff_Cre_1_1.fq.gz'. In this case 'r1base' and 'r2base' are invalid.

Edited May 02, 2019 by Jeon Lee
Assignee
Assign to
None
Milestone
None
Assign milestone
Time tracking
None
Due date
None
Reference: BICF/Astrocyte/rnaseq#13