tombo
tombo copied to clipboard
TOMBO MAPPING ERROR
Hello, I am running this tombo command with etect_modifications
[00:54:29] Loading minimap2 reference.
[00:54:47] Getting file list.
[00:55:52] Loading default canonical ***** DNA ***** model.
[00:55:58] Re-squiggling reads (raw signal to genomic sequence alignment).
100%|██████████| 5910216/5910216 [182:29:30<00:00, 9.00it/s]
******************** WARNING ********************
Unexpected errors occured. See full error stack traces for first (up to) 50 errors in "unexpected_tombo_errors.4056.err"
[15:25:29] Final unsuccessful reads summary (37.6% reads unsuccessfully processed; 2222411 total reads):
17.2% (1016283 reads) : Read event to sequence alignment extends beyond bandwidth
11.0% ( 650175 reads) : Alignment not produced
9.4% ( 555418 reads) : Poor raw to expected signal matching (revert with tombo filter clear_filters
)
0.0% ( 223 reads) : Too much raw signal for mapped sequence
0.0% ( 171 reads) : Read failed sequence-based signal re-scaling parameter estimation.
0.0% ( 44 reads) : Unexpected error
0.0% ( 41 reads) : Not enough raw signal around potential genomic deletion(s)
0.0% ( 37 reads) : Fewer changepoints found than requested
0.0% ( 19 reads) : Read contains too many potential genomic deletions
[15:25:29] Saving Tombo reads index to file.
/public/home/lishigui/anaconda3/envs/deepsignalpenv/lib/python3.6/site-packages/numpy/core/fromnumeric.py:3373: RuntimeWarning: Mean of empty slice.
out=out, **kwargs)
/public/home/lishigui/anaconda3/envs/deepsignalpenv/lib/python3.6/site-packages/numpy/core/_methods.py:170: RuntimeWarning: invalid value encountered in double_scalars
ret = ret.dtype.type(ret / rcount)
Any ideas why this happened?
Thanks, Lee