I am using a slurm script to submit Trinity assembly of my dataset. Can you please comment on what might be wrong? Is it still Java?
#!/bin/sh
#SBATCH --job-name=Trinity_32 # Job name
#SBATCH -N 1 # Maximum CPU cores for job
#SBATCH -n 12 #all cores from same node
#SBATCH --mem=20gb # Memory limit
#SBATCH --time=2-00:00:00 # Time limit days-hrs:min:sec
#SBATCH -e errorlog_%j.err # standard error
#SBATCH --output=loggeneraltrinity32 # Standard output and error log
for file in *R1.fq
do
. ~/.profile
withpath="${file}"
filename=${withpath##*/}
base="${filename%*_*.fq}"
echo "${base}"
time Trinity --seqType fq --bypass_java_version_check --max_memory 18G --left "${base}"_R1.fq --right "${base}"_R2.fq --CPU 12 --no_normalize_reads --KMER_SIZE 32 --output general_trinity_"${base}"_32
done;
--------------------------------------------------------
-------------------- Chrysalis -------------------------
-- (Contig Clustering & de Bruijn Graph Construction) --
--------------------------------------------------------
inchworm_target: /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/both.fa
bowite_reads_fa: /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/both.fa
chrysalis_reads_fa: /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/both.fa
* Running CMD: /home/hasnahana/softwares/trinityrnaseq-Trinity-v2.4.0/util/support_scripts/filter_iworm_by_min_length_or_cov.pl /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/inchworm.K32.L32.DS.fa 100 10 > /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/inchworm.K32.L32.DS.fa.min100 * Running CMD: bowtie2-build -o 3 /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/inchworm.K32.L32.DS.fa.min100 /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/inchworm.K32.L32.DS.fa.min100 1>/dev/null
* Running CMD: bash -c " set -o pipefail;bowtie2 --local -k 2 --threads 12 -f --score-min G,46,0 -x /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/inchworm.K32.L32.DS.fa.min100 /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/both.fa | samtools view -@ 12 -F4 -Sb - | samtools sort -m 805306368 -@ 12 -no - - > /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/iworm.bowtie.nameSorted.bam"
139641624 reads; of these:
139641624 (100.00%) were unpaired; of these:
281744 (0.20%) aligned 0 times
57569891 (41.23%) aligned exactly 1 time
81789989 (58.57%) aligned >1 times
99.80% overall alignment rate
[bam_sort_core] merging from 96 files...
[W::bam_hdr_read] EOF marker is absent. The input is probably truncated.
[E::bgzf_read] bgzf_read_block error -1 after 0 of 4 bytes
[bam_merge_core] error: './samtools.63291.393.tmp.0054.bam' is truncated.
Error, cmd: bash -c " set -o pipefail;bowtie2 --local -k 2 --threads 12 -f --score-min G,46,0 -x /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/inchworm.K32.L32.DS.fa.min100 /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/both.fa | samtools view -@ 12 -F4 -Sb - | samtools sort -m 805306368 -@ 12 -no - - > /scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/headcropped_assembly/som_assembly/general_trinity_15TRIM_SOM_32/chrysalis/iworm.bowtie.nameSorted.bam" 2>tmp.62365.stderr died with ret 256 at /home/hasnahana/softwares/trinityrnaseq-Trinity-v2.4.0/PerlLib/Pipeliner.pm line 166
Pipeliner::run('Pipeliner=HASH(0xa5c030)') called at /home/hasnahana/softwares/trinityrnaseq-Trinity-v2.4.0/Trinity line 1750
main::run_chrysalis('/scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/...', '/scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/...', 200, 500, undef, '/scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/...', '/scratch/hasnahana/transcriptome/som/combinedOldNew/assembly/...') called at /home/hasnahana/softwares/trinityrnaseq-Trinity-v2.4.0/Trinity line 1595
main::run_Trinity() called at /home/hasnahana/softwares/trinityrnaseq-Trinity-v2.4.0/Trinity line 1262
eval {...} called at /home/hasnahana/softwares/trinityrnaseq-Trinity-v2.4.0/Trinity line 1261
Trinity run failed. Must investigate error above.