Mercurial > repos > drosofff > fetch_fasta_from_ncbi
comparison retrieve_fasta_from_NCBI.py @ 10:2c5375809c03 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/fetch_fasta_from_ncbi commit e0a1114b735bf1af257456174f64e5ef8d205754-dirty
| author | drosofff |
|---|---|
| date | Wed, 28 Oct 2015 11:24:45 -0400 |
| parents | cc43a7a11324 |
| children | 0bec3cba5c56 |
comparison
equal
deleted
inserted
replaced
| 9:f57ebdc7ad56 | 10:2c5375809c03 |
|---|---|
| 164 counter += 1 | 164 counter += 1 |
| 165 self.logger.info("Server Transaction Trial: %s" % ( counter ) ) | 165 self.logger.info("Server Transaction Trial: %s" % ( counter ) ) |
| 166 try: | 166 try: |
| 167 response = urllib2.urlopen(req) | 167 response = urllib2.urlopen(req) |
| 168 fasta = response.read() | 168 fasta = response.read() |
| 169 if "Resource temporarily unavailable" in fasta: | 169 if ("Resource temporarily unavailable" in fasta) or (not fasta.startswith(">") ): |
| 170 serverTransaction = False | 170 serverTransaction = False |
| 171 else: | 171 else: |
| 172 serverTransaction = True | 172 serverTransaction = True |
| 173 except urllib2.HTTPError as e: | 173 except urllib2.HTTPError as e: |
| 174 serverTransaction = False | 174 serverTransaction = False |
| 175 self.logger.info("urlopen error:%s, %s" % (e.code, e.read() ) ) | 175 self.logger.info("urlopen error:%s, %s" % (e.code, e.read() ) ) |
| 176 except httplib.IncompleteRead as e: | 176 except httplib.IncompleteRead as e: |
| 177 serverTransaction = False | 177 serverTransaction = False |
| 178 self.logger.info("IncompleteRead error: %s" % ( e.partial ) ) | 178 self.logger.info("IncompleteRead error: %s" % ( e.partial ) ) |
| 179 if self.dbname != "pubmed": | |
| 180 assert fasta.startswith(">"), fasta | |
| 181 fasta = self.sanitiser(self.dbname, fasta) # | 179 fasta = self.sanitiser(self.dbname, fasta) # |
| 182 time.sleep(1) | 180 time.sleep(1) |
| 183 return fasta | 181 return fasta |
| 184 | 182 |
| 185 def sanitiser(self, db, fastaseq): | 183 def sanitiser(self, db, fastaseq): |
