ShortRead::FastqStreamer() causes segfault
0
0
Entering edit mode
@watal-m-iwasaki-6651
Last seen 5.8 years ago
Japan

Hi all,

I am trying to filter FASTQs with library(ShortRead). One of my FASTQs contains an entry with width=0 and it causes segfault in the following code:

it = FastqStreamer('reads_1.fasta.gz', 10)
srq = yield(it)
as(quality(srq), 'matrix')

I think I could remove those reads in advance with srq[width(srq) > 0], but it seems to spoil the advantage of srFilter() and compose(). Interestingly, I found that an ShortReadQuality object read by readFastq() did not reproduce this error. So splitting FASTQ to small pieces may be another workaround. Do you have any suggestion?

 

Thanks,
Watal

shortread segfault • 773 views
ADD COMMENT
0
Entering edit mode

I wasn't able to reproduce this with a fastq file I created artificially; can you please update your question to include the output of sessionInfo(), as well as providing a way to reproduce your data (e.g., by creating a zero-width record after running example(readFastq), writing out using writeFastq(), or by making your data available via dropbox or referencing it in a public data base or communicating with me martin.morgan at roswellpark dot org directly).

One you have updated your question, please use the 'ADD COMMENT' link to indicate that the question has been updated.

ADD REPLY

Login before adding your answer.

Traffic: 223 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6