Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by gconcepcion View PostIt seems as though pulse features were not added into your *.cmp.h5
What was the process you used to generate your cmp.h5. Did you use pbalign with the --forQuiver option?Code:tofu_wrap.py --nfl_fa isoseq_nfl.fasta --ccs_fofn reads_of_insert.fofn --bas_fofn input.fofn -d clusterOut --quiver --gmap_db /zs32/data-analysis/liucy_group/llhuang/Reflib/gmapdb --gmap_name gmapdb_h19 --output_seqid_prefix human isoseq_flnc.fasta final.consensus.fa
Last edited by GenoMax; 07-28-2016, 09:18 AM.
Comment
-
Originally posted by bowhan View Post
reads_of_insert.fofn
Comment
-
Originally posted by bowhan View Post
Comment
-
Originally posted by lingling huang View PostI am sure that the fofn files are correct. I try it again, Just as expected, error occurred I can't understand the cause of the error displayed on screen, as usual, I attach them at below. As you can see, the code has been running for half.Could you explain why this is the case? Thank you!
[ATTACH]4447[/ATTACH]
[ATTACH]4448[/ATTACH]
[ATTACH]4449[/ATTACH]
This one is on DALIGNER itself. The version shipped within ToFU is pretty outdated. Several bugs have been found and fixed since then.
Please obtain the more recently version from https://github.com/PacificBiosciences/DALIGNER/, compile it and replace the ones under your virtualenv/bin directory.
Comment
-
Originally posted by lingling huang View Postinput.fofn
[ATTACH]4445[/ATTACH]
reads_of_insert.fofn
[ATTACH]4446[/ATTACH]
Try removing the blank lines at the end and re-running.Last edited by gconcepcion; 07-29-2016, 12:36 PM.
Comment
-
Originally posted by bowhan View PostTo me this is a different error from the one described originally.
This one is on DALIGNER itself. The version shipped within ToFU is pretty outdated. Several bugs have been found and fixed since then.
Please obtain the more recently version from https://github.com/PacificBiosciences/DALIGNER/, compile it and replace the ones under your virtualenv/bin directory.
error:
Comment
-
Originally posted by lingling huang View PostHi, bowhan . I have updated DALIGNER files, and ran again. Finally, I obtain the 'combined' folder in the directory 'clusterOut' . It seems that I basically got the result, but I found it was lack of "all_sizes.quivered_hq.fastq.5merge.collapsed.XXX " files. Disappointed, I failed again. I guess there is a problem in the step of 'collapse_isoforms_by_sam.py‘ command or in the process of mapping to the reference genome, right? How can I fix it up?
Please try the method described in this thread. If the error remains, re-install ToFU under a new directory. Make sure you follow the instructions (https://github.com/PacificBioscience...ranscript-tofu) strictly.
Additionally, I remember that you had issues with outdated Cython. Consider upgrading that BEFORE installing ToFU.
Originally posted by lingling huang View PostAnother question is why I can't find 'final.consensus.fa' in current pathway since I wrote it in my command. As usual, I attach there pictures at below. I am looking forward to your help, thank you very much!
https://github.com/PacificBioscience...ICE-and-Quiver)
BTW, I wouldn't use the default bin size, which is 1kb. Reads in different bins are not seeing each other during the clustering process. Use something like `--bin_manual "(0,2,4,6,8,99)"`.
Comment
-
Originally posted by bowhan View PostThe error message indicates a (solved) issue described here:
Please try the method described in this thread. If the error remains, re-install ToFU under a new directory. Make sure you follow the instructions (https://github.com/PacificBioscience...ranscript-tofu) strictly.
Additionally, I remember that you had issues with outdated Cython. Consider upgrading that BEFORE installing ToFU.
That's fine. `final.consensus.fa` is the output of `pbtranscript.py cluster` but not ToFU. You should be expecting a file named `all_sizes.quivered_hq.fastq` under the `combined` folder. For more information, please read this:
https://github.com/PacificBioscience...ICE-and-Quiver)
BTW, I wouldn't use the default bin size, which is 1kb. Reads in different bins are not seeing each other during the clustering process. Use something like `--bin_manual "(0,2,4,6,8,99)"`.
Comment
Latest Articles
Collapse
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
22 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment