Conversation
There was a problem hiding this comment.
Pull request overview
Updates the trimming step so that UMI information stored in the BAM RX tag can be carried through BAM→FASTQ conversion (and then through fastp) by explicitly copying RX into FASTQ output.
Changes:
- Bump
fastpconda dependency from1.0.1to1.1.0. - Remove the
fastp_pathoverride and invokefastpdirectly from the environment. - Update
samtools fastqinvocation to attempt to copyRXtags into FASTQ (-T RX).
Reviewed changes
Copilot reviewed 2 out of 2 changed files in this pull request and generated 2 comments.
| File | Description |
|---|---|
| modules/fastp.nf | Update fastp/samtools invocation to preserve RX tag during BAM→FASTQ→fastp trimming. |
| conf/test.config | Remove unused fastp_path param from the test profile. |
💡 Add Copilot custom instructions for smarter, more guided reviews. Learn how to get started.
|
|
||
| samtools fastq -n ${bam} | \\ | ||
| ${fastp_path}fastp --stdin \\ | ||
| samtools fastq -T RX -n ${bam} | \\ |
| samtools fastq -T RX -n ${bam} | \\ | ||
| fastp --stdin \\ |
There was a problem hiding this comment.
Probably the test is failing because the snapshot needs to be updated for the fastp version change (which gets us back to the bioconda fastp rather than our local install since the adapter dimer fix was merged). But I'm waiting on the bwameth change -- until then this UMI fix won't work and it breaks.
|
i re-ran the tests... i think they are failing due to the fastp version difference: "\t\t"fastp_version": "1.0.1",", | "\t\t"fastp_version": "1.1.0",", |
@bwlang Yes, I agree, thanks! UMI fix won't work without the bwameth change. I'm still hoping that goes into conda soon, but if we need it before then a local version of bwameth could be used. |
relies on this change to bwameth, so will need to wait for bioconda release or use another installation.
brentp/bwa-meth@2bf5ae3