Skip to content

Error during index built #456

Description

@desmodus1984

Hi,
I am interested in doing genome annotation, and Braker3 mentions Hisat2 to map RNA reads to a genome assembly.
I installed Hisat2 in conda, and then I used the following script to make the index for the genome:
hisat2-build -p 64 -f HydLep.fx.fasta HydLep
Using 64 cores is too much?
I was expecting it to run fast, and when I checked the error/log of the job in PBS,
The log looks okay-ish:
Sorting block time: 00:00:12
Returning block of 6318003 for bucket 457
Sorting block time: 00:00:05
Returning block of 2123596 for bucket 460
Sorting block time: 00:00:08
Returning block of 4302114 for bucket 463
Sorting block time: 00:00:09
Returning block of 4797812 for bucket 465
Sorting block time: 00:00:08
Returning block of 4785495 for bucket 466
Sorting block time: 00:00:12
Returning block of 6871697 for bucket 459
Sorting block time: 00:00:09
Returning block of 4793335 for bucket 467
Sorting block time: 00:00:11
Returning block of 6138098 for bucket 462
Sorting block time: 00:00:11
Returning block of 6187927 for bucket 461
Sorting block time: 00:00:13
Returning block of 7085826 for bucket 464

But the error mentions an issue:
offMask: 0xfffffff0
ftabChars: 10
eftabLen: 0
eftabSz: 0
ftabLen: 1048577
ftabSz: 4194308
offsLen: 153398404
offsSz: 613593616
lineSz: 64
sideSz: 64
sideGbwtSz: 48
sideGbwtLen: 192
numSides: 12783201
numLines: 12783201
gbwtTotLen: 818124864
gbwtTotSz: 818124864
reverse: 0
linearFM: Yes
Total time for call to driver() for forward index: 00:42:16
/var/spool/torque/mom_priv/jobs/1022481.scyld.localdomain.SC: line 43: syntax error: unexpected end of file

Since I was expecting it to finish fast, I added a loop to start mapping the RNA data to the genome,

for i in /data/common/juanpablo.aguilar/analysis/RNA/*_1_trim.fastq;
        do
        base=$(basename $i "_1_trim.fastq")

hisat2 --phred33 --dta -x HydLep -S ${base}.sam \
        -1 /data/common/juanpablo.aguilar/analysis/RNA/${base}_1_trim.fastq -2 /data/common/juanpablo.aguilar/analysis/RNA/${base}_2_trim.fastq \

samtools view -b -h -q 20 -F 4 ${base}.sam > ${base}.bam

but I am not sure that the index creating worked well. I was expecting a clear, explicit, message that the index was done, and/or that the loop failed.

Could you tell me how to find if the index was properly created?

Thank you.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions