Ka Ming Nip
Ka Ming Nip
Thanks for reporting this! Yes, this happens when there are too few reads.
I was able to replicate this, but this is not a bug. The assembled sequences are too short and they all end up in `rnabloom.transcripts.short.fa` (instead of `rnabloom.transcripts.fa`). I have...
Not at all. The length threshold is the only determining factor for assigning sequences to these two files.
There is already an option for that (i.e. `-length`) and its default value is 200, which is what separates the sequences in the two files. All RNA-seq assemblers I can...
Sorry, I thought you were asking whether RNA-Bloom use any evidence to determine that threshold.
The sum isn't the same because each base of a given read can be assigned to more than one contig, especially when it aligns to the overlap between two contigs....
In the above diagram, contigs 1 and 2 can still be join together into a transcript if there is sufficient overlap at the contig edges. Double-counting can also happen when...
Can you please print the first few lines of `rnabloom_k17.hist`? e.g. ``` head rnabloom_out/rnabloom_k17.hist ```
This number `9223372036854775808` is just larger than the largest 64-bit long integer in Java, i.e. `9,223,372,036,854,775,807`. Can you please report your ntCard version? ``` ntcard --version ```
This looks like a bug in ntCard. Can you share your read sequences with me?