Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
Chris Boursnell
@cboursnell
@beig2048 No it doesn't, sorry. Reads need to be accurately assigned to contigs and there isn't enough information in single end reads to be able to do that.
jcornel3
@jcornel3
OK, I just renamed that one library just in case, but that seems to have solved the problem and it's working great! Thanks. Also, I'm using Ubuntu 16.04
Aaron
@aaron_mckenna_twitter
@cboursnell thanks; I had fixed it by bumping the SNAP version used, so that makes sense.
jcornel3
@jcornel3
Hey all, a couple quick questions: does Transrate allow for interleaved reads as input? How about globbing read fastas for --left and --right? For example, could the inputs be --left left.fq --right right.fq? Or, do they actually need to be comma separated lists? Thanks!
jcornel3
@jcornel3
lol, just realized that my attempts at writing an asterisk failed miserably in my question.
Chris Boursnell
@cboursnell
@jcornel3 Transrate doesn't allow for interleaved reads or globbing I'm afraid.
@jcornel3 you can type asterisks by doing \* btw and you can also edit you message by clicking the three dot menu on the right of each message
jcornel3
@jcornel3
Ok, thanks! For some reason "edit" is grayed out.
Chris Boursnell
@cboursnell
You can only edit in the first few minutes after posting a message
Chenming Cui
@cui_cm_twitter
@blahah I filtered my input reads into 1million, 2million...5 million to assemble my transcriptom, then I tried to evaluate these 5 transcriptom by transRate, transRate can generate output CSV files, but just didn't give "score", "optimal score", cutoff and weight, showing "NA"... can you point out what is the problem? Actually transRate works with my non-filtered reads ones... I am wondering if these 5 transcriptom were assembly with too less reads and in very poor quality then transRate can not rate them in score.
Chenming Cui
@cui_cm_twitter
@cboursnell I filtered my input reads into 1million, 2million...5 million to assemble my transcriptom, then I tried to evaluate these 5 transcriptom by transRate, transRate can generate output CSV files, but just didn't give "score", "optimal score", cutoff and weight, showing "NA"... can you point out what is the problem? Actually transRate works with my non-filtered reads ones... I am wondering if these 5 transcriptom were assembly with too less reads and in very poor quality then transRate can not rate them in score.
countloki
@countloki
So I deleted librt.so.1 and TransRate seemed to be running fine. However I noticed the salmon.log was 4.6GB of error messages like this:
[2017-03-16 00:14:21.704] [jointLog] [warning] [Sampler.hpp]: Failed to sample an alignment for this read; this shouldn't happen
currentMass = -nan, r = 0.677047
I'm afraid maybe some will miss this and use erroneous stats
I find my own results suspicious anyway as they are paradoxical to my Detonate RSEM scores.
Still hoping for a solution to the salmon not installed message. Working from Ubuntu 16.04 LTS.
MGCBrown
@MGCBrown
Hello. I'm looking to evaluate my assembly with RNAseq reads. If I have paired end reads from multiple individuals, do I have concatenate the reads into a single --left.fq and --right.fq file, or can I specify multiple --left and --right fastq files separated by a comma? Thanks.
Rik
@blahah
@MGCBrown you can pass in multiple files, like --left file1_left.fq,file2_left.fq --right file1_right.fq,file2_right.fq
the list of files to each argument should be comma separated, no spaces
you can also provide multiple assemblies and/or reference the same way
Rik
@blahah
@countloki thanks for reporting - we'll make a new release soon that will solve the salmon issues
Stuart Willis
@stuartwillis
@blahah and @cboursnell have you guys had any luck finding why salmon gives the "numeric overflow" error on some assemblies (issue #7 on the Transfuse github)? I have found it is reproducible running salmon on directly on the bam file.
Stuart Willis
@stuartwillis
at least one assembly also repeatably gives the following error, also from salmon: "terminate called after throwing an instance of 'std::bad_alloc'"
Rob Patro
@rob-p
@stuartwillis, any idea if that is still true with a newer salmon?
The version in transrate is very old at this point. Lots of improvements and bug fixes since.
@cboursnell @blahah, I'm booked this week, but perhaps we could work on integrating a newer salmon (and fixing other bugs) the following week, since some ppl seem to have recurring issues with 1.0.3
Stuart Willis
@stuartwillis
@rob-p, good question, I didn't try; I can check that if would justify upgrading transrate to a newer salmon
Rob Patro
@rob-p
Great. If you still see the error, please file a bug report for salmon so I can see if I can repro the issue upstream.
Stuart Willis
@stuartwillis
Indeed a more recent version of salmon (0.8.2) was successful on a bam/assembly that repeatedly resulted in a "numeric overflow" with 0.4.2.
However, I also found in the same case that deleting the transrate results (bam and all) and starting over circumvented the problem, while remaking the assembly (transabyss, in this case) circumvented this error in another case. Not sure what to make of that
It does speak to the prescience of @cboursnell to make transfuse able to pick up mid-way, though
Rik
@blahah
@rob-p definitely we should do that - I've been swamped for like 10 months, but this week I've finished up a project so I should have free cycles for it
wow we are super out of sync with salmon
@rob-p how close are we to being able to skip SNAP altogether? in-memory streaming alignment -> assignment -> assessment is the dream we talked about a year ago
Rob Patro
@rob-p
We are ~90% of the way done I'd say. I have 2 students working on this (which also includes some algorithmic advances to rapmap).
Their stuff adds the ability to "selectively" align quasi-mappings. Given that, it'll be a few days of work on the salmon end to enable the posterior sampling with the new code.
Rik
@blahah
@rob-p sweet - that will be such a great update
Stuart Willis
@stuartwillis
if you guys are tweaking the implementation of salmon into transrate, would it be worth it to add a specification for strand specific data? Not sure how salmon's behavior changes w or w/o on stranded data
jdeyar
@jdeyar

Hello @blahah,

Referring to other posts regarding an issue with 1.0.3, I am having a similar SNAP error while running transrate. My error code is noted below.

Welcome to SNAP version 1.0dev.96.

mmap failed.
MemMapDataReader: fail to map /pylon2/bi4s8fp/jsdsmail/Transrate/Trim_TS6_T24.left.pe.qc.fq at 0,3765724304
SNAP exited with exit code 1 from line 2283 of file SNAPLib/DataReader.cpp

Any help/advice on this error would be much appreciated. I got the same error when trying to run it on the previous version of transrate as well.

Thanks!

Rik
@blahah
@jdeyar hmm sorry about that, I'm going to try to work through the transrate issues next week (I'm no longer paid to work on it, and have been swamped with other stuff for a while)
@jdeyar hopefully this is one we can get dealt with quickly by just updating SNAP and Salmon and cleaning up a few packaging issues
jdeyar
@jdeyar
@blahah thank you for your prompt response. Hopefully it is as easy as an update for SNAP and other linked software. Please comment in this forum whenever it is updated! Much appreciated!
XinpengQi
@XinpengQi

@blahah any updates so for? I am running Transrate v1.0.3. The algorithms for contig/assembly score are so attractive! Eager to try them out in my data! SNAP seems fine, but Salmon gives warning "[jointLog] [warning] [Sampler.hpp]: Failed to sample an alignment for this read; this shouldn't happen
currentMass = -nan, r = 0.694954"

snap.log:

Loading index from directory... 2s. 384271831 bases, seed size 23
Aligning.
Total Reads Aligned, MAPQ >= 10 Aligned, MAPQ < 10 Unaligned Too Short/Too Many Ns Extra Alignments %Pairs Reads/s Time in Aligner (s)
693,340,846 280,486,471 (40.45%) 330,208,110 (47.63%) 78,660,098 (11.35%) 3,986,167 (0.57%) 620,752,208 78.36% 124,387 5,574
Welcome to SNAP version 1.0dev.96.

salmon.log

[2017-05-17 13:40:55.996] [fileLog] [info] quantification processed 354000000 fragments so far
[2017-05-17 13:41:01.787] [jointLog] [info]
Completed first pass through the alignment file.
Total # of mapped reads : 354762625
of uniquely mapped reads : 144015998
ambiguously mapped reads : 210746627
[2017-05-17 13:41:01.761] [fileLog] [info] quantification processed 354762625 fragments so far
[2017-05-17 13:41:33.214] [jointLog] [info] Computed 9628374 rich equivalence classes for further processing
[2017-05-17 13:41:33.214] [jointLog] [info] Counted 354762625 total reads in the equivalence classes
[2017-05-17 13:41:33.227] [jointLog] [info] starting optimizer
[2017-05-17 13:41:39.341] [jointLog] [info] Marked 0 weighted equivalence classes as degenerate
[2017-05-17 13:41:39.618] [jointLog] [info] iteration = 0 | max rel diff. = 800.763
[2017-05-17 13:41:53.311] [jointLog] [info] iteration 50, recomputing effective lengths
[2017-05-17 13:42:09.716] [jointLog] [info] iteration = 100 | max rel diff. = 0.358082
[2017-05-17 13:42:40.040] [jointLog] [info] iteration = 200 | max rel diff. = 0.305365
[2017-05-17 13:43:10.754] [jointLog] [info] iteration = 300 | max rel diff. = 0.0754343
[2017-05-17 13:43:41.135] [jointLog] [info] iteration = 400 | max rel diff. = 0.0490107
[2017-05-17 13:44:09.308] [jointLog] [info] iteration 500, recomputing effective lengths
[2017-05-17 13:44:10.037] [jointLog] [info] iteration = 500 | max rel diff. = 0.0219555
[2017-05-17 13:44:39.741] [jointLog] [info] iteration = 600 | max rel diff. = 0.127602
[2017-05-17 13:45:09.829] [jointLog] [info] iteration = 699 | max rel diff. = 0.0077808
[2017-05-17 13:45:09.864] [jointLog] [info] finished optimizer
[2017-05-17 13:45:10.734] [jointLog] [info] Sampling alignments; outputting results to ./postSample.bam
[2017-05-17 13:45:11.509] [jointLog] [warning] [Sampler.hpp]: Failed to sample an alignment for this read; this shouldn't happen
currentMass = -nan, r = 0.694954

Rik
@blahah
cc @rob-p \o any ideas?
Rob Patro
@rob-p
Yes (cc @blahah), this is a bug that can occur in VB inference that was fixed in a more recent version of Salmon. I fear the only immediate solution here may be to generate the relevant file manually with a more recent Salmon release --- though we should probably upgrade the one included in Transrate soon
btw, the preprint that describes the relevant algorithm for having salmon do alignment internally went up yesterday --- http://biorxiv.org/content/early/2017/05/17/138800
I think that we can do a patch release before that, but this (rather big new feature) could land soon.
Rik
@blahah
so @rob-p if I cut a new transrate with updated salmon it will fix this?
Rob Patro
@rob-p
This should be the case, yes --- however, before you do that, we should really test the new pipeline on one or two data sets first. A lot has been updated, but the module to sample alignments hasn't been touched in a while.
Rik
@blahah
ok