Hello all, I’m currently running trinity on 3 reps of PE leaf tissue RNA from a species with no genome sequence in order to create a transcriptome for sequence analysis.
This is my first try at running Trinity on our universities computing cluster.
Here is my Trinity input code, and I’m running on a compute node with 1.5 TB memory:
Trinity --seqType fq --SS_lib_type RF --max_memory 1500G --left $raw/NM1_PG1_S64_R1.fastq,$raw/NM2_PG2_S65_R1.fastq,$raw/NM3_PG3_S66_R1.fastq --right $raw/NM1_PG1_S64_R2.fastq,$raw/NM2_PG2_S65_R2.fastq,$raw/NM3_PG3_S66_R2.fastq --CPU 6 --no_bowtie --output $trin
I’ve gotten through read normalization (~10% of reads kept, which for a large genome plant I think I expect because there are probably a lot of very highly expressed similar genes), and inchworm (I think?), and the output file said it’s started working through chrysalis. However, for many hours the output keeps giving me something like this:
If it indicates bad_alloc(), then Inchworm ran out of memory. You’ll need to either reduce the size of your data set or run Trinity on a server with more memory available. The inchworm process failed.warning, cmd: /software/sl-7.x86_64/modules/trinity/2.5.1/util/support_scripts/../../Trinity –single “//Trinity_PG/read_partitions/Fb_1/CBin_1253/c125330.trinity.reads.fa” –output “/Trinity_PG/read_partitions/Fb_1/CBin_1253/c125330.trinity.reads.fa.out” –CPU 1 –max_memory 1G –run_as_paired –SS_lib_type F –seqType fa –trinity_complete –full_cleanup –no_bowtie failed with ret: 256, going to retry. succeeded(10076), failed(3324) 10.3695% completed.
Files are being produced in the “read_partitions” directory. Is this fine or should I be concerned about these warnings?
Read more here: Source link