vertebrata[Organism]) NOT insects[Organism]) AND ("200"[Modification Date] : "3000"[Modification Date])''' # Maximum number of returned results. retmax = 5000 # Instantiate search object. sra_search = SRASearch(query=query, retmax=retmax, email=email_bruno) # Execute search itself. sra_search.esearch() # Fetch metadata from packages. packages = [SRAPackage(sra_id) for sra_id in sra_search.idlist] # Store packages in data frame for filtering. package_filter = FilterPackages(packages) # copy working data frame. df = package_filter.data_frame # Filter booleans. filtered_df = df[df.library_layout == 'PAIRED'][df.nreads > 1][df.read_average >= 70] # Sort data buy lineage. sorted_df = filtered_df.sort('lineage') # Write CSV out. package_filter.filtered_data_frame = sorted_df package_filter.write_csv(basename + '.csv') # Write unique list of taxa.