Dear all,
I ran an alignment of Sequel data using blasr. Now I want to get some metrics about the sam file using samstat 1.5.1 (running on Ubuntu 16.04 LTS; 16 Gb of RAM).
I get a message that reads:
[2017-08-01 13:26:25] SAMstat 1.5.1, Copyright (C) 2014 Timo Lassmann <[email protected]>
[2017-08-01 13:26:25] cmd: samstat Ecoli_2K.sam
[2017-08-01 13:26:25] Ecoli_2K.sam
--------------------------------------------------
[2017-08-01 13:26:25] Starting to collect data.
[2017-08-01 13:26:25] Training a HMM on mapq > 20 reads.
Segmentation fault (core dumped)
I suspect the problem is the file size (>150,000 reads, file size: 13.5 Gb), so I trimmed the sam file using "head -n". I am only successful at getting a report when I have 20 reads...
{I have no problem getting a report when I work with Illumina-based alignments}
Is this a PacBio-specific problem?
Is the problem the memory? is there a way to increase the RAM usage of samstat (similar to a -vm function in some Java applications)?
is there another solution? Are there other programs out there that can do the same as samstat 1.5 (the reports are very cool!)?
Thanks,
Gustavo
PS: there is the same kind of problem if the alignment is done with BWA MEM...
I ran an alignment of Sequel data using blasr. Now I want to get some metrics about the sam file using samstat 1.5.1 (running on Ubuntu 16.04 LTS; 16 Gb of RAM).
I get a message that reads:
[2017-08-01 13:26:25] SAMstat 1.5.1, Copyright (C) 2014 Timo Lassmann <[email protected]>
[2017-08-01 13:26:25] cmd: samstat Ecoli_2K.sam
[2017-08-01 13:26:25] Ecoli_2K.sam
--------------------------------------------------
[2017-08-01 13:26:25] Starting to collect data.
[2017-08-01 13:26:25] Training a HMM on mapq > 20 reads.
Segmentation fault (core dumped)
I suspect the problem is the file size (>150,000 reads, file size: 13.5 Gb), so I trimmed the sam file using "head -n". I am only successful at getting a report when I have 20 reads...
{I have no problem getting a report when I work with Illumina-based alignments}
Is this a PacBio-specific problem?
Is the problem the memory? is there a way to increase the RAM usage of samstat (similar to a -vm function in some Java applications)?
is there another solution? Are there other programs out there that can do the same as samstat 1.5 (the reports are very cool!)?
Thanks,
Gustavo
PS: there is the same kind of problem if the alignment is done with BWA MEM...
Comment