Closed jennieli421 closed 3 weeks ago
max_chunks
has to be greater than min_chunks
, and should not be set to 1.
You also don't seem to have set any targets?
We have been using max_chunks=1
for a very long time and the results were fine. Here I am only providing targeted region's reference, so I don't need to specify targets
Here I am only providing targeted region's reference, so I don't need to specify targets
I'm not quite sure what you are trying to do here. If you do not give any targets then you will not get any unblocking activity. An empty
targets = [] will mean that no targets are checked and so the later logic of on or off target will not be evaluated.
An example toml enriching for two chromosomes from a reference is available here: https://github.com/LooseLab/readfish/blob/main/docs/_static/example_tomls/human_chr_selection.toml
In short, you must have something in the targets = [] for any unblocking activity to occur.
I have been using this setting (with split region) for a long time and it is working as expected. The design for not having targets can be traced back to this conversation: https://github.com/LooseLab/readfish/issues/242#issuecomment-1610971660
Basically I have provided references of only the desired sequences, so everything mapped to the reference is accepted. Thus, the problem I want to fix is why this setting does not work when I only have one region.
Hmmm
I indeed stand corrected.
Let me check this further and come back to you.
Hi, Could you pleaase confrim your readfish version and also share your readfish log with us - the log should be called the
Thanks
Matt
tested_one_or_more_regions_ef170d0f.zip
Unfortunately, I was testing different region settings during the run, so this run had multiple settings. I hope this is fine.
Hi @jennieli421 - do these logs contain an example of it working, with a region of control=true
, where you see unblocking?
And can I just check the version of readfish
you were using?
readfish --version
Final question, what is the mean read length of the library you are testing?
readfish --version
readfish 2023.1.1
Both should contain multiple trials where I see unblocking, as long as there is a region of control=True
.
This issue is stale because it has been open 30 days with no activity. Remove stale label or comment or this will be closed in 5 days.
This issue was closed because there has been no response for 5 days after becoming stale.
I wanted to perform a run with only enrichment, so I had one region specified in TOML like the following. But I did not see any unblock actions from the live state panel. The purple peak indicating voltage reversals did not appear at all in the histogram. However, everything works fine as long as I add a region with
control=true
. I also tried with two regions withcontrol=false
and that also didn't have any rejection.