-
Notifications
You must be signed in to change notification settings - Fork 10
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
Chipseq updates #111
base: main
Are you sure you want to change the base?
Chipseq updates #111
Conversation
@adthrasher can you split the |
Moved to #119 |
* main: feat: Style updates (#122) fix: calc_gene_length algorithm. And better docs for gene length and TPM (#120) Iterative samtools merge (#118) chore: move from pr #111 (#119) feat: experimental dnaseq-standard WF (#116) docs: new stuff for style-guide, best-practices, and template (#117) docs: More docs (#114) fix(qc): finish deleting gc_bias_metrics feat(qc): drop requirement of reference_fasta feat: put strandedness value into HTSeq file name (#90) chore: bump all kraken images to latest (#113)
max_retries=max_retries | ||
} | ||
String rg_id_field = sub(sub(pair.right, ".*ID:", "ID:"), "\t.*", "") | ||
String rg_id_field = sub(sub(tuple.right, ".*ID:", "ID:"), "\t.*", "") | ||
String rg_id = sub(rg_id_field, "ID:", "") | ||
call samtools.addreplacerg as single_end { input: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Rename single_end
to something else
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
poke
"MISSING_PLATFORM_VALUE", | ||
"INVALID_PLATFORM_VALUE", | ||
"INVALID_MAPPING_QUALITY", | ||
"MATES_ARE_SAME_END", | ||
"MISMATCH_FLAG_MATE_NEG_STRAND", | ||
"MISMATCH_MATE_ALIGNMENT_START" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The first two seem "ok" to ignore. Less comfortable ignoring the last four. Are all these still being produced by seaseq?
max_retries=max_retries | ||
} | ||
String rg_id_field = sub(sub(pair.right, ".*ID:", "ID:"), "\t.*", "") | ||
String rg_id_field = sub(sub(tuple.right, ".*ID:", "ID:"), "\t.*", "") | ||
String rg_id = sub(rg_id_field, "ID:", "") | ||
call samtools.addreplacerg as single_end { input: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
poke
Updates to reflect changes to SEAseq and add support for SEAseq's paired-end mode.