comparison bdss_client_sra_pe.Rmd @ 16:c2150a0f9c7c draft

planemo upload for repository https://github.com/statonlab/docker-GRReport/tree/master/my_tools/rmarkdown_bdss_client commit ea6df97960569d247d64be1549aa90768c9786f4-dirty

15:fe5f30584308

    # build URL from SRA accession
    url = paste0('ftp://ftp.ncbi.nlm.nih.gov/sra/sra-instant/reads/ByRun/sra/',
                 substr(id, 1, 3), '/',
                 substr(id, 1, 6), '/', id, '/', id, '.sra')
    # download sra file with bdss
    bdss_command = paste0(' /main/sites/galaxy/galaxy/tools/_conda/bin/bdss transfer -u ', url)
    system(bdss_command, intern = TRUE)
    # convert .sra to .fastq/.fasta
    if('FORMAT' == 'fasta') {
      command = paste0('fastq-dump --fasta --split-files ', '-O read_files_directory ', id)
    } else {

    if (length(all_fastq_files) == 1) {
      # this is a single end SRA file. build url to download fastq
      print(paste0('only one fastq file found for this SRA accession: ', id))
      print(paste0('Downloading skipped for ', id))
      # url = paste0(url_base, '/', all_fastq_files)
      # bdss_command = paste0(' /main/sites/galaxy/galaxy/tools/_conda/bin/bdss transfer -u ', url)
      # run_bdss = system(bdss_command, intern = TRUE)
      # print(run_bdss)
    } else {
      print(paste0('Two fastq files found for this SRA accession: ', id))
      url_1 = paste0(url_base, '/', all_fastq_files[1])
      bdss_command_1 = paste0(' /main/sites/galaxy/galaxy/tools/_conda/bin/bdss transfer -u ', url_1)
      run_bdss_1 = system(bdss_command_1, intern = TRUE)
      print(run_bdss_1)
      url_2 = paste0(url_base, '/', all_fastq_files[1])
      bdss_command_2 = paste0(' /main/sites/galaxy/galaxy/tools/_conda/bin/bdss transfer -u ', url_2)
      run_bdss_2 = system(bdss_command_2, intern = TRUE)
      print(run_bdss_2)
    }
  }

16:c2150a0f9c7c

    # build URL from SRA accession
    url = paste0('ftp://ftp.ncbi.nlm.nih.gov/sra/sra-instant/reads/ByRun/sra/',
                 substr(id, 1, 3), '/',
                 substr(id, 1, 6), '/', id, '/', id, '.sra')
    # download sra file with bdss
    bdss_command = paste0('/tool_deps/_conda/bin/bdss transfer -u ', url)
    system(bdss_command, intern = TRUE)
    # convert .sra to .fastq/.fasta
    if('FORMAT' == 'fasta') {
      command = paste0('fastq-dump --fasta --split-files ', '-O read_files_directory ', id)
    } else {

    if (length(all_fastq_files) == 1) {
      # this is a single end SRA file. build url to download fastq
      print(paste0('only one fastq file found for this SRA accession: ', id))
      print(paste0('Downloading skipped for ', id))
      # url = paste0(url_base, '/', all_fastq_files)
      # bdss_command = paste0('/tool_deps/_conda/bin/bdss transfer -u ', url)
      # run_bdss = system(bdss_command, intern = TRUE)
      # print(run_bdss)
    } else {
      print(paste0('Two fastq files found for this SRA accession: ', id))
      url_1 = paste0(url_base, '/', all_fastq_files[1])
      bdss_command_1 = paste0('/tool_deps/_conda/bin/bdss transfer -u ', url_1)
      run_bdss_1 = system(bdss_command_1, intern = TRUE)
      print(run_bdss_1)
      url_2 = paste0(url_base, '/', all_fastq_files[1])
      bdss_command_2 = paste0('/tool_deps/_conda/bin/bdss transfer -u ', url_2)
      run_bdss_2 = system(bdss_command_2, intern = TRUE)
      print(run_bdss_2)
    }
  }