Trouble with SOLiD reads

[BioFalcon]

Hi Mike,

I have been trying to run ShortStack using SOLiD reads, but it seems to be having trouble reading the reads on the get_read_s_aln subroutine, here is the last part of the log.

Adapter trimming input reads file Sample1.csfasta.gz with adapter CTGTAGGC Complete: 9711 out of 9994211 reads successfully trimmed
FATAL in sub-routine get_read_s_aln : Could not determine format of reads file ShortStack_Sample1/Sample1.c_trimmedsfasta.gz

Any idea of what might be the problem?

Thanks

[MikeAxtell]

Wow, SOLiD data. Haven't thought about SOLiD in many years and that part of the code base has not been maintained at all. One thing I noticed from your output is a very small percentage of input reads that were successfully trimmed. Is that what you expected? It seems like an issue (although not one related to this bug). I don't have any example SOLiD data to work with here. I think for me to properly de-bug this, I'd need some SOLiD data. Do you want to share yours (you could post on dropbox or whatever), and let me know what genome assembly you're working with? That way I could reproduce the bug here and try to fix it. Thanks, Mike

…

On Wed, Sep 26, 2018 at 5:23 PM Francisco Falcon ***@***.***> wrote: Hi Mike, I have been trying to run ShortStack using SOLiD reads, but it seems to be having trouble reading the reads on the get_read_s_aln subroutine, here is the last part of the log. Adapter trimming input reads file Sample1.csfasta.gz with adapter CTGTAGGC Complete: 9711 out of 9994211 reads successfully trimmed FATAL in sub-routine get_read_s_aln : Could not determine format of reads file ShortStack_Sample1/Sample1.c_trimmedsfasta.gz Any idea of what might be the problem? Thanks — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#82>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AGiXial1DzyLOz4p3RR1NfLMNKHof6Zfks5ue_BTgaJpZM4W7g38> .

-- Michael J. Axtell, Ph.D. Professor of Biology Penn State University http://sites.psu.edu/axtell

[MikeAxtell]

Will be dropping support for colorspace data altogether in the next major release.