Thanks for helping Felix. A simple typo on my end unfortunately. Problem solved!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
For those with weaker heart (those that cannot use complex scripts in Linux and need a graphic interface) here is another (free) program for trimming qualities:
An efficient SFF/FastQ viewer and editor (GUI)
The gray/green curves in the second graphic shows the average quality before and after trimming the low quality ends.
Only works on Fasta, FastQ, SFF for the moment.
Sorry.Last edited by create.share; 05-30-2015, 12:19 AM.
Comment
-
Hi everyone,
not sure if this is the right thread but I'll give it a go:
I'm trying to get rid of adapters in my sequences. There aren't many with adapters but I'd like to save what I can.
When I run my command I get the following error:
==================================================================
Your job looked like:
------------------------------------------------------------
# LSBATCH: User input
/gpfs/scratch/cbh12wsu/trim/trim_galore -phred33 --illumina --paired -phred33 --path_to_cutadapt /gpfs/scratch/cbh12wsu/trim/cutadapt.py -o adaptrim.fastq -length 50 1_L.fastq 2_R.fastq
------------------------------------------------------------
Exited with exit code 2.
The output (if any) follows:
Path to Cutadapt set as: '/gpfs/scratch/cbh12wsu/trim/cutadapt.py' (user defined)
File "/gpfs/scratch/cbh12wsu/trim/cutadapt.py", line 99
print(rest, match.read.name, file=self.file)
^
SyntaxError: invalid syntax
Cutadapt seems to be working fine (tested command '/gpfs/scratch/cbh12wsu/trim/cutadapt.py --version')
File "/gpfs/scratch/cbh12wsu/trim/cutadapt.py", line 99
print(rest, match.read.name, file=self.file)
^
SyntaxError: invalid syntax
Failed to write to file '1_L.fastq_trimming_report.txt': No such file or directory
==================================================================
Does indicate a problem with the cutadapt.py?
Thanks!
Comment
-
Hmm, when you just type "/gpfs/scratch/cbh12wsu/trim/cutadapt.py" on the command line do you see something like:
Code:cutadapt version 1.8.1 Copyright (C) 2010-2015 Marcel Martin <[email protected]> cutadapt removes adapter sequences from high-throughput sequencing reads. Usage: cutadapt -a ADAPTER [options] [-o output.fastq] input.fastq
Comment
-
Hi Felix,
thanks for replying so quick!
I get the following error when I try your command:
File "/gpfs/scratch/cbh12wsu/trim/cutadapt.py", line 99
print(rest, match.read.name, file=self.file)
^
SyntaxError: invalid syntax
I think I should mention that I made cutadapt.py an executable. And copied it from it's original folder into the trim folder.
Comment
-
I put it back where it belongs and tried again. Same error as before.
With normal installation you mean the pip command? I tried that before but I'm not allowed to do that on the cluster. I can load a cutadapt module but I was hoping making it an executable would solve the problem.
Comment
-
I've installed it and when I run felix's command I get:
cutadapt version 1.8.1
...
But I still get this error:
Path to Cutadapt set as: '/gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt' (user defined)
1.8.1
Cutadapt seems to be working fine (tested command '/gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt --version')
Failed to write to file '1_L.fastq_trimming_report.txt': No such file or directory
Used this command:
bsub -N -q short -o cut3.txt "/gpfs/scratch/cbh12wsu/trim/trim_galore -phred33 --illumina --paired -phred33 --path_to_cutadapt /gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt -o adaptrim.fastq -length 50 1_L.fastq 2_R.fastq"
More ideas?
Comment
-
no that shouldn't matter... Could it be that your queuing system does not take you to the same directory (option like -cwd)? What happens if you just run the command without using the queue?
By the way you can shorten the command like this:
Code:/gpfs/scratch/cbh12wsu/trim/trim_galore --paired --path_to_cutadapt /gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt [COLOR="Red"]-o adaptrim.fastq[/COLOR] -length 50 1_L.fastq 2_R.fastq
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-25-2024, 11:49 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
04-25-2024, 11:49 AM
|
||
Started by seqadmin, 04-24-2024, 08:47 AM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
04-24-2024, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
62 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
61 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
Comment