Skip to content
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

update guppy #204

Open
replikation opened this issue Jan 27, 2022 · 4 comments
Open

update guppy #204

replikation opened this issue Jan 27, 2022 · 4 comments
Assignees
Labels
enhancement New feature or request

Comments

@replikation
Copy link
Owner

new guppy version 6 does not require arrangement files for demultiplexing and has some changes

demulti like this

process collect_fastq {
        label 'demultiplex'
    input:
        tuple val(name), val(technology), path(dir)
    output:
        tuple val(name), path("*.fastq.gz"), emit: reads
    script:
        if (params.single)
        """
        find -L ${dir} -name '*.fastq' -exec cat {} +  | gzip > ${name}.fastq.gz
        find -L ${dir} -name '*.fastq.gz' -exec zcat {} + | gzip >> ${name}.fastq.gz
        """
        else if (!params.single)
        """
        BARCODE_DIRS=\$(find -L ${dir} -name "barcode??" -type d)
        
        if [ -z "\${BARCODE_DIRS}" ]; then 
            guppy_barcoder -t ${task.cpus} -r -i ${dir} -s fastq_outbreak  \
                    --detect_mid_strand_barcodes \
                    --min_score_barcode_mid 50 \
                    --trim_adapters \
                    --trim_barcodes \
                    --disable_pings  

            for barcodes in fastq_outbreak/barcode??; do
                find -L \${barcodes} -name '*.fastq' -exec cat {} + | gzip >> \${barcodes##*/}.fastq.gz
            done
        else
            for barcodes in \${BARCODE_DIRS}; do
                find -L \${barcodes} -name '*.fastq' -exec cat {} + | gzip >> \${barcodes##*/}.fastq.gz
                find -L \${barcodes} -name '*.fastq.gz' -exec zcat {} + | gzip >> \${barcodes##*/}.fastq.gz
            done
        fi
        """
        stub:
        """
        touch ${name}.fastq.gz
        """
}


@replikation replikation added the enhancement New feature or request label Jan 27, 2022
@replikation
Copy link
Owner Author

also needs to be updated with the most recent medaka - artic version and model

@hoelzer
Copy link
Collaborator

hoelzer commented Feb 24, 2023

I want to bump this issue because we now have a Medaka update in the ARTIC process. It would be good to have updated Guppy so people can use the R1041 models for basecalling and the matching Medaka models.

@replikation should be relatively straight-forward, or? We should just do proper testing where we can also assist.

@replikation
Copy link
Owner Author

they changed some flags with the new guppy for adapter trimming so it not be a drop in replacment

@hoelzer
Copy link
Collaborator

hoelzer commented Feb 24, 2023

I see, so would need proper adjustment and testing for all the different supported settings. But

docker pull nanozoo/guppy_gpu:6.4.6-1--2c17584

seems quite recent for testing.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants