Skip to content

Segmentation fault errors #16

Description

@jingzhejiang

Hi,

I run virmap on c5d.24xlarge (96 vCPU, 192 GiB, 4* 800Gb) instance with an 838Gb SWAP setting

# swap setting
$ sudo mkswap -f /dev/nvme3n1; \
$ sudo swapon /dev/nvme3n1
$ free -h
              total        used        free      shared  buff/cache   available
Mem:           188G         77G         33G         74M         76G        109G
Swap:          838G        2.2G        835G

It works for many libraries, but didn't for another similar size libraries (maybe with too many viral contigs ? I'm not sure.)
For the error libraries, they threw out many empty files after 12S-WGA.filtered.fa step.

$ ll -th 12S-WGA/
-rw-r--r-- 1 ec2-user ec2-user 2.0K 11月  5 13:28 12S-WGA.err
-rw-r--r-- 1 ec2-user ec2-user 1.3K 11月  5 13:28 12S-WGA.log
-rw-r--r-- 1 ec2-user ec2-user  607 11月  5 13:28 12S-WGA.taxonomy.err.bz2
-rw-r--r-- 1 ec2-user ec2-user    0 11月  5 13:27 12S-WGA.final.fa
-rw-r--r-- 1 ec2-user ec2-user    0 11月  5 13:27 12S-WGA.diamondBlastx.out.bz2
-rw-r--r-- 1 ec2-user ec2-user 1.3K 11月  5 13:27 12S-WGA.diamondBlastx.err
-rw-r--r-- 1 ec2-user ec2-user   14 11月  5 13:27 12S-WGA.blastn.out.bz2
-rw-r--r-- 1 ec2-user ec2-user   34 11月  5 13:27 12S-WGA.blastn.err
-rw-r--r-- 1 ec2-user ec2-user  915 11月  5 13:27 12S-WGA.selfAlign.err.bz2
-rw-r--r-- 1 ec2-user ec2-user    0 11月  5 13:27 12S-WGA.selfAlign.fa
-rw-r--r-- 1 ec2-user ec2-user 1.4K 11月  5 13:27 12S-WGA.iterateImprove.err.bz2
-rw-r--r-- 1 ec2-user ec2-user    0 11月  5 13:27 12S-WGA.improved.fa
-rw-r--r-- 1 ec2-user ec2-user 6.6K 11月  5 13:27 12S-WGA.filter.err.bz2
-rw-r--r-- 1 ec2-user ec2-user    0 11月  5 13:21 12S-WGA.filtered.fa
-rw-r--r-- 1 ec2-user ec2-user 283M 11月  5 13:21 12S-WGA.blastnFilter.out.bz2
-rw-r--r-- 1 ec2-user ec2-user  58M 11月  5 11:56 12S-WGA.diamondFilter.out.bz2
-rw-r--r-- 1 ec2-user ec2-user  70M 11月  5 11:44 12S-WGA.combined.fa
-rw-r--r-- 1 ec2-user ec2-user 1.9M 11月  5 11:44 12S-WGA.combine.err
-rw-r--r-- 1 ec2-user ec2-user 1.6K 11月  5 11:34 12S-WGA.assembly.err.bz2
-rw-r--r-- 1 ec2-user ec2-user  62M 11月  5 11:34 12S-WGA.contigs.fa
-rw-r--r-- 1 ec2-user ec2-user 6.1M 11月  5 09:04 12S-WGA.superScaffolds.err.bz2
-rw-r--r-- 1 ec2-user ec2-user 9.7M 11月  5 09:04 12S-WGA.superScaffolds.fa
-rw-r--r-- 1 ec2-user ec2-user  11M 11月  5 08:58 12S-WGA.superScaffolds.init.fa
-rw-r--r-- 1 ec2-user ec2-user 2.2M 11月  5 08:58 12S-WGA.pseudoScaffolds.fa
-rw-r--r-- 1 ec2-user ec2-user 200K 11月  5 08:52 12S-WGA.both.centroids.txt
-rw-r--r-- 1 ec2-user ec2-user 221M 11月  5 08:49 12S-WGA.aa.sam.bz2
-rw-r--r-- 1 ec2-user ec2-user 3.0K 11月  5 08:26 12S-WGA.bbmap.err
-rw-r--r-- 1 ec2-user ec2-user 8.9M 11月  5 08:26 12S-WGA.nuc.sam.bz2
-rw-r--r-- 1 ec2-user ec2-user 341M 11月  5 08:12 12S-WGA.normalized.fa.bz2
-rw-r--r-- 1 ec2-user ec2-user 2.3K 11月  5 08:12 12S-WGA.normalize.err
-rw-r--r-- 1 ec2-user ec2-user 1.7G 11月  5 08:07 12S-WGA.derep.fa.bz2
-rw-r--r-- 1 ec2-user ec2-user  397 11月  5 08:07 12S-WGA.derep.err
-rw-r--r-- 1 ec2-user ec2-user 2.5G 11月  5 07:39 12S-WGA.fa.bz2

Here is the record of error log file

TIME 12S-WGA decompress: 166.73 seconds
TIME 12S-WGA dereplicate: 1692.31 seconds
TIME 12S-WGA normalize: 281.11 seconds
TIME 12S-WGA bbmap to virus: 822.88 seconds
TIME 12S-WGA diamond to virus: 1364.85 seconds
TIME 12S-WGA construct superscaffolds: 913.65 seconds
TIME 12S-WGA megahit assembly: 9022.49 seconds
TIME 12S-WGA dedupe assembly: 17.12 seconds
TIME 12S-WGA merge assembly: 580.39 seconds
TIME 12S-WGA diamond filter map: 693.05 seconds
TIME 12S-WGA blastn filter map: 5120.28 seconds
sh: line 1: 72606 Segmentation fault      determineTaxonomy.pl filter /scratch/VirmapDb/Taxonomy.virmap /scratch/tmp/MEUE7kr6ri/12S-WGA.entropy.fa /scratch/Virmap_result/12S-WGA/12S-WGA.diamondFilter.out.bz2 /scratch/Virmap_result/12S-WGA/12S-WGA.blastnFilter.out.bz2 96 > /scratch/Virmap_result/12S-WGA/12S-WGA.filtered.fa 2>> /scratch/tmp/MEUE7kr6ri/12S-WGA.filter.err
TIME 12S-WGA filter contigs: 377.65 seconds
FAlite: Empty
TIME 12S-WGA iterative improvement: 6.71 seconds
TIME 12S-WGA self align and quantify: 1.19 seconds
FAlite: Empty
TIME 12S-WGA blastn full: 0.02 seconds
FAlite: Empty
cat: /scratch/tmp/MEUE7kr6ri/12S-WGA.diamondBlastx.out: No such file or directory
FAlite: Empty
FAlite: Empty
FAlite: Empty
lbzip2: skipping "/scratch/tmp/MEUE7kr6ri/12S-WGA.diamondBlastx.out": lstat(): No such file or directory
lbzip2: skipping "/scratch/tmp/MEUE7kr6ri/12S-WGA.diamondBlastx.remain.out": lstat(): No such file or directory
TIME 12S-WGA diamond full: 0.04 seconds
TIME 12S-WGA determine taxonomy: 26.89 seconds
TIME 12S-WGA Overall Virmap time: 21087.35 seconds

I also checked the end of file 12S-WGA.filter.err.bz2, and found error messages at the last line:

...................
THREAD 77 - GOT contig_2641;src=megahit;length=2873;cov=53.3;gc=0.646;taxId=10239
THREAD 90 - GOT gi|NoNucHits|gb|MF418016.1|Bacillus.phage.vB_BceM.HSE3;taxId=2170705;length=168394
THREAD 16 - GOT contig_9275;src=megahit;length=2881;cov=51;merged=1;taxId=10239
THREAD 58 - GOT contig_19539;src=megahit;length=1666;cov=54.1;gc=0.627;taxId=10239
THREAD 67 - GOT gi|NEWDB|gb|KY935474.1|Influenza.A.virus.A.Porto.Alegre.LACENRS.1155.2015.H3N2.segment.7.matrix.protein.2.M2.and.matrix.protein.1.M1.genes.complete.cds.;taxId=1979338;length=1047
Thread 94 terminated abnormally: IO error: While opening a file for sequentially reading: /dev/shm/vBuhJeZAAe/TAX/CURRENT: Too many open files at /usr/local/bin/determineTaxonomy.pl line 458.

Is this a memory error, even after an 838Gb SWAP is set? Can I overcome the problem of memory error by adjusting SWAP setting, or I can only change to some expensive, big memory instance?

ANOTHER IMPORTANT QUESTION is: Can virmap continue to execute from the interrupted steps, thus saving some time and money?

Thank you for your reply!

12S-WGA.err.gz
12S-WGA.filter.err.gz

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