Date:
01/09/09 04:20
|
Submitted by:
glyn_dk
|
Assigned to:
glyn_dk
|
Category:
ALG file handling
|
Priority:
9
|
Ticket group:
reproducible
|
Resolution:
Resolved
|
Summary:
import Barski wigfiles, only one track is loaded to graphs
|
Original submission:
- Load multiple wigfiles (varstep) from the Barski ChIP-seq data
- Load a bedfile
- 'Convert and add to CATCH'
BUG:
Three pwig files are generated, but only one track is shown in the graphs view.
No exceptions or error messages are shown.
Update Jan 9, 2009 by glyn_dk:
I have located the error. In the the Barski wigfiles the tracks do not have signal values for all positions. Since profiles of less than two signal values are ignored, the pwig files may contain different numbers of profiles.
When adding pwigs to dataset files that do not have the correct number of profiles are ignored. I have added an error message to this special case: "ERROR: mismatch in number of profiles in selected tracks"
In general it should be possible for the user to use these tracks together, so we have to construct a workaround.
Suggestion for workaround:
- the generation of the pwig files shall remain the same
- the user can choose either of two options when loading data: 1) discard profiles with no signal values defined 2) pad with zeros in tracks where no signal values are defined.
|
|
Please log in to add comments and receive followups via email.
|
|
Followups
|
Comment
|
Date
|
By
|
Implemented workaround:
- pwig files now include a profile name ( e.g. # chr2:220006943-220008943) for all selected positions, regardless of whether the selected region has signal values defined within it.
- when readSequences reads a sequence (a profile region) with less than two signals defined, two dummy signals are inserted with value 0.
- the interpolation of the zero values gives a track with zero signal at all reference positions.
Now users can define any crazy region or load tracks with lots of missing data (e.g. ChIP-seq tracks) without encountering an error message or runtime exception.
|
01/11/09 07:24
|
glyn_dk
|
|
HIGH PRIORITY: Must be fixed before next release!
|
01/09/09 04:22
|
glyn_dk
|
|
|
No results for "Dependent on ticket" |
No results for "Dependent on Task" |
|
No other tickets are dependent on this ticket
|
|
Ticket change history
|
Field
|
Old value
|
Date
|
By
|
|
status_id
|
Unset |
01/11/09 07:24
|
glyn_dk
|
|
assigned_to
|
unset |
01/11/09 07:24
|
glyn_dk
|
|
resolution_id
|
Unset |
01/11/09 07:24
|
glyn_dk
|
|
close_date
|
12/31/69 19:00 |
01/11/09 07:24
|
glyn_dk
|
|
priority
|
5 |
01/09/09 04:22
|
glyn_dk
|
|