You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
New optional param motifs for rule sc2rf_recombinants.
New param weeks for new rule plot.
Removed prev_linelist param.
Output
Switch from a pdf report to powerpoint slides for better automation.
Create summary plots.
Split report rule into linelist and report.
Output svg plots.
Workflow
New rule plot.
Changed growth calculation from a comparison to the previous week to a score of sequences per day.
Assign a cluster_id according to the first sequence observed in the recombinant lineage.
Define a recombinant lineage as a group of sequences that share the same:
Lineage assignment
Parents
Breakpoints or phylogenetic placement (subtree)
For some sequences, the breakpoints are inaccurate and shifted slightly due to ambiguous bases. These sequences can be assigned to their corresponding cluster because they belong to the same subtree.
For some lineages, global prevalence has exceeded 500 sequences (which is the subtree size used). Sequences of these lineages are split into different subtrees. However, they can be assigned to the correct cluster/lineage, because they have the same breakpoints.
Confirmed not to use deletions define recombinants and breakpoints (differs from published)?
Programs
Move sc2rf_recombinants.py to postprocess.py in ktmeaton fork of sc2rf.
Add false positives filtering to sc2rf_recombinants based on parents and breakpoints.