I knew about (1), but then how can I use the parameters to run it? and I had the trimmed data already. For (2) it worked, but then I got a memory error message; for (3) how shall I feed the files to Trinity? -- left pairedA,pairedB --right unpairedA,unpairedB? Or simply --left pairedA --right pairedB?
Trinity assembly from Trimmomatic output
Hello,
I am running Trinity from paired sequences which have been trimmed from low quality and adaptors using Trimmomatic. The output of Trimmomatic is formed by 4 files -- two paired and two unpaired -- and I understand Trinity can get multiple files as comma separated. Here are the commands I wrote:
java -jar /usr/share/java/trimmomatic-0.33.jar PE -threads 16 -phred33 source.1.fq source.2.fq pairedA pairedB unpairedA unpairedB ILLUMINACLIP:NexteraPE-PE.fa:2:30:10:1:true LEADING:5 TRAILING:5 MINLEN:36
followed by
Trinity --seqType fq --max_memory 10G --no_bowtie --no_version_check --left pairedA, pairedB --right unpairedA, unpairedB --CPU 10
but I got this error:
Sunday, December 6, 2015: 09:30:45 CMD: java -Xmx64m -XX:ParallelGCThreads=2 -jar /usr/local/lib/Trinity/util/support_scripts/ExitTester.jar 0
Sunday, December 6, 2015: 09:30:46 CMD: java -Xmx64m -XX:ParallelGCThreads=2 -jar /usr/local/lib/Trinity/util/support_scripts/ExitTester.jar 1
Error, cannot locate file: at /usr/bin/Trinity line 2150.
main::create_full_path(ARRAY(0x1102cf0), 1) called at /usr/bin/Trinity line 1106
which I think is due to the wrong input of for the Trinity command; any tips on how to run this process properly?
Thank you
• 4,457 views
•
link
1 answer
- Trinity can run Trimmomatic for you with
--trimmomatic - The cause for your error is you have a space between file names on the Trinity call.
- If if you fix the space, your Trinity command still will be wrong, you can not pair
pairedA.fqwithunpairedA.fq- this is what you are doing with your--leftand--rightparameters
• 0 views
•
link
Log in to answer this question.