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

BBmap_Sendsketch not finding SRA #85

Open
buchanri opened this issue Jul 18, 2024 · 0 comments
Open

BBmap_Sendsketch not finding SRA #85

buchanri opened this issue Jul 18, 2024 · 0 comments

Comments

@buchanri
Copy link

After running and getting this error, sendsketch is not getting a hit for a SRA I chose. It then moves on to this step but errors because it didn't get any hits for the SRA to begin with.

ERROR ~ Error executing process > 'PATHOGENSURVEILLANCE:PREPARE_INPUT:INITIAL_CLASSIFICATION (SRA_NCBISRR8356409)'

Caused by:
  Process `PATHOGENSURVEILLANCE:PREPARE_INPUT:INITIAL_CLASSIFICATION (SRA_NCBISRR8356409)` terminated with an error exit status (1)


Command executed:

  sendsketch_filter.R SRA_NCBISRR8356409.txt
  
  KINGDOM="$(cat kingdom.txt)"
  CLASS="$(cat classification.txt)"
  
  mv species.txt SRA_NCBISRR8356409_species.txt
  mv genera.txt SRA_NCBISRR8356409_genera.txt
  mv families.txt SRA_NCBISRR8356409_families.txt
  
  cat <<-END_VERSIONS > versions.yml
  "PATHOGENSURVEILLANCE:PREPARE_INPUT:INITIAL_CLASSIFICATION":
      r-base: $(echo $(R --version 2>&1) | sed 's/^.*R version //; s/ .*$//')
  END_VERSIONS

Command exit status:
  1

Command output:
  (empty)

Command error:
  Error in writeLines(data$taxonomy[1], "classification.txt") : 
    can only write character objects
  Execution halted

Work dir:
  /media/changlab/changlab4/paradarc/paper1_bra/scripts/pathogensurveillance/work/62/831f3690c47314df4733c9d3915a94

Tip: you can replicate the issue by changing to the process work dir and entering the command `bash .command.run`

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

Here is the BBmap_sendsketch work output:

Query: SRR8356409.fastq.gz      DB: RefSeq      SketchLen: 8    Seqs: 118       Bases: 17818    gSize: 61       GC: 0.545       Quality: 0.3483 AvgCount: 5.375 Depth: 5.354    File: SRR8356409.fastq.gz
No hits.

work/b6/7cc4098376958ceb73c5573facf2f6/SRA_NCBISRR8356409.txt 

But it should be able to find it because on the the ncbi website it shows that it is a good option. https://www.ncbi.nlm.nih.gov/sra?linkname=bioproject_sra_all&from_uid=492229

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

No branches or pull requests

1 participant