Skip to content

ERROR ~ Error executing process > 'DENSE_WF:DENSE:EXTRACT_CDS (18)' #7

@Xillanne

Description

@Xillanne

I'm using fna and gff3 dowload via dataset from NCBI. Three of my genomes invoque the same issue, I don't know what to do. Can you please help me.

Caused by:
Process DENSE_WF:DENSE:EXTRACT_CDS (18) terminated with an error exit status (134)

Command executed:

chmod -R +x /home/desmarais/.nextflow/assets/i2bc/dense/bin

Use GffRead.

Need a fai index

samtools faidx Mayamaea_pseudoterrestris.fna

Keep GFF sequences that are also in the FASTA file.

awk 'BEGIN{FS=OFS=" "} {print "^"$1,""} END {print "^#"}' Mayamaea_pseudoterrestris.fna.fai | grep -f - Mayamaea_pseudoterrestris.gff3 > gff_filterA

Also remove mRNA with undefined strand and features with abnormal end

awk -F" '.next" '
FNR==NR {max[$1]=$2}
FNR!= NR && ( /^#/ || ($7 !~ /?/ && $5 <= max[$1]) )
' Mayamaea_pseudoterrestris.fna.fai gff_filterA > gff_filterB

Extract the genomic CDS.

-V discard any mRNAs with CDS having in-frame stop codons

-x : write a fasta file with spliced CDS for each GFF transcript

gffread -V -g Mayamaea_pseudoterrestris.fna -x Mayamaea_pseudoterrestris_CDS.fna gff_filterB

Remove CDS missing a terminal stop codon and get a translated version of the FASTA.

discard_CDS_missing_terminal_stop_codon.sh Mayamaea_pseudoterrestris_CDS.fna

Command exit status:
134

Command output:
(empty)

Command error:
free(): invalid next size (normal)
.command.sh: line 21: 72 Aborted gffread -V -g Mayamaea_pseudoterrestris.fna -x Mayamaea_pseudoterrestris_CDS.fna gff_filterB

Work dir:
/data/desmarais/PlastiDia/DENSE/work/88/e57bcaa319205ab5843415784d4dac

Container:
/home/desmarais/.nextflow/assets/i2bc/dense/containers/lopesi2bc-denovogenes.img

Tip: when you have fixed the problem you can continue the execution adding the option -resume to the run command line

-- Check '.nextflow.log' file for details

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions