Guest User

Untitled

a guest
Jun 21st, 2016
197
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. import subprocess, os
  2. from ftplib import FTP
  3.  
  4. ftp_site = 'ftp.ncbi.nlm.nih.gov'
  5. ftp = FTP(ftp_site)
  6. ftp.login()
  7. ftp.cwd('genomes/genbank/bacteria')
  8. dirs = ftp.nlst()
  9. for organism in dirs:
  10.     latest = os.path.join(organism, "latest_assembly_versions")
  11.     for path in ftp.nlst(latest):
  12.         accession = path.split("/")[-1]
  13.         fasta = accession+"_genomic.fna.gz"
  14.         subprocess.call(['rsync',
  15.                          '--recursive',
  16.                          '-vv',
  17.                          '--copy-links',
  18.                          '--dry-run',
  19.                          '-f=+ '+fasta,
  20.                          #'-f=- *',
  21.                          'ftp.ncbi.nlm.nih.gov::genomes/genbank/bacteria/'+latest,
  22.                          '--log-file=scratch/test_dir/log.txt',
  23.                          'scratch/' + organism])
Add Comment
Please, Sign In to add comment