Originally posted by roryk
View Post
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by fongchun View PostDid you end up writing your own script that merges the junction.bed files from different sample while preserving the score? If yes, would you be willing to share that script since this is exactly what I need to do.
cat junctions_file1.bed junction_file_2.bed junction_file_3.bed| tbed2juncs | combineJuncs > combined.juncs
combine.juncs will be a BED file and contains all of the junctions from the junction files with the score the sum of the junction scores. The names of the junctions are the locations of the two nucleotides at the edges of the exons that are joined together.
Leave a comment:
-
Did you end up writing your own script that merges the junction.bed files from different sample while preserving the score? If yes, would you be willing to share that script since this is exactly what I need to do.
Leave a comment:
-
I'm new just thought I'd try helping. Apparently it isn't the bed file format I am familiar with.
Leave a comment:
-
Originally posted by DrD2009 View PostYou could use Galaxy (http://main.g2.bx.psu.edu/). Upload the files (probably as file format tabular) you want to combine and then use the 'Concatenate queries' tool found under Text Manipulation. There's probably a short script out there to combine them, but this should work for what you want to do too.
-Brandon
Leave a comment:
-
Originally posted by Cole Trapnell View PostTopHat comes with an (undocumented, I realize) script called bed_to_juncs. Running it on a TopHat BED file will produce a .juncs-format file. To merge these, you can simply cat them together, sort, and pipe to uniq.
I am aware of the bed_to_juncs script but when converting BED files to .junc files, it doesn't preserve the score and overhang information that I would like to combine. BTW, could you confirm if my understanding of the score field in BED files is correct? Thanks a lot!
-- LeoLast edited by HTS; 12-14-2009, 02:27 PM.
Leave a comment:
-
TopHat comes with an (undocumented, I realize) script called bed_to_juncs. Running it on a TopHat BED file will produce a .juncs-format file. To merge these, you can simply cat them together, sort, and pipe to uniq.
Leave a comment:
-
You could use Galaxy (http://main.g2.bx.psu.edu/). Upload the files (probably as file format tabular) you want to combine and then use the 'Concatenate queries' tool found under Text Manipulation. There's probably a short script out there to combine them, but this should work for what you want to do too.
-Brandon
Leave a comment:
-
How to combine junctions.bed files produced by TopHat
Hi,
Basically I would like to double check before writing my own script to do so. Since I need to pool together samples with different read lengths, I have to run TopHat separately for them (using the -j option wherever appropriate). I already know how to merge the resulting .sam files with samtools and generate a combined coverage.wig file. Before putting effort to combine the junctions.bed files, I would like to know:
1. If there are tools/scripts to do this already, either within TopHat or outside.
2. If none, I would like to confirm if the score field of the junctions.bed file is simply the number of uniquely mapping reads that are aligned to the junction, or if muitlreads are also counted/weighted.
Please feel free to share your knowledge/experience/comments. Thanks a lot!
-- LeoTags: None
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
25 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
29 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
24 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
|
Leave a comment: