-
Notifications
You must be signed in to change notification settings - Fork 5
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Consistent use of primers throughout pipeline #102
Comments
Identifying which primers belonging to an amplicon are observed among all fragments belonging to that amplicon is a subtly different task than identifying which primer appears in a fragment. I agree that it would have been nice to retain the original consistent logic but the solution that we arrived at ultimately requires that these operations follow different logic. |
IIRC the plan was always to remove primers in the way I just described. And be consistent throughout. |
There has been a large amount of discussion and planning about how we identify amplicons and which primers are observed from them and how we identify which bases belong to primers inside a given fragment. These are quite different problems. If you would like to modify the behaviour of the amplicon coordinates that are passed to cylon during the consensus step, then let's decouple that request from the adapter issues. |
This issue is about
|
Yes, that is what we do |
yeah the point is, we excliude it when passing to cylon, but it comes back by the time we get to selfQC |
Expected behaviour: the amplicon scheme and the primers used are inferred early on in the pipeline. Some primers can be identified as not present in the reads. These primers should be completely excluded from everything onwards.
Actual behaviour: inconsistent use of whether primers are present or not.
See comment "The logic for identifying primers before cylon is run is different than the later primer id logic that we use to mask primer positions in the consensus." #99 (comment)
The text was updated successfully, but these errors were encountered: