# HG changeset patch
# User public-health-bioinformatics
# Date 1573613085 18000
# Node ID 4f963b583186353bc6e784469a77ffa208738276
# Parent b2d56a44a872eb253d099ec67dea7d85bc9934aa
"planemo upload for repository https://github.com/public-health-bioinformatics/galaxy_tools/blob/master/tools/screen_abricate_report commit c27095ceaf23aa2c910e0ea866db9ec8d3e7816b"
diff -r b2d56a44a872 -r 4f963b583186 screen_abricate_report.py
--- a/screen_abricate_report.py Thu Oct 31 15:43:15 2019 -0400
+++ b/screen_abricate_report.py Tue Nov 12 21:44:45 2019 -0500
@@ -3,11 +3,9 @@
from __future__ import print_function
import argparse
-import os
+import csv
import re
-import sys
-import csv
-from pprint import pprint
+
def parse_screen_file(screen_file):
screen = []
@@ -17,21 +15,26 @@
screen.append(row)
return screen
+
def get_fieldnames(input_file):
with open(input_file) as f:
reader = csv.DictReader(f, delimiter="\t", quotechar='"')
row = next(reader)
fieldnames = row.keys()
return fieldnames
-
+
+
def main(args):
screen = parse_screen_file(args.screening_file)
- abricate_report_fieldnames = get_fieldnames(args.abricate_report)
gene_detection_status_fieldnames = ['gene_name', 'detected']
- with open(args.abricate_report, 'r') as f1, open(args.screened_report, 'w') as f2, open(args.gene_detection_status, 'w') as f3:
+ with open(args.abricate_report, 'r') as f1, \
+ open(args.screened_report, 'w') as f2, \
+ open(args.gene_detection_status, 'w') as f3:
abricate_report_reader = csv.DictReader(f1, delimiter="\t", quotechar='"')
- screened_report_writer = csv.DictWriter(f2, delimiter="\t", quotechar='"', fieldnames=abricate_report_reader.fieldnames)
- gene_detection_status_writer = csv.DictWriter(f3, delimiter="\t", quotechar='"', fieldnames=gene_detection_status_fieldnames)
+ screened_report_writer = csv.DictWriter(f2, delimiter="\t", quotechar='"',
+ fieldnames=abricate_report_reader.fieldnames)
+ gene_detection_status_writer = csv.DictWriter(f3, delimiter="\t", quotechar='"',
+ fieldnames=gene_detection_status_fieldnames)
screened_report_writer.writeheader()
gene_detection_status_writer.writeheader()
@@ -45,14 +48,16 @@
gene_detection_status['detected'] = True
screened_report_writer.writerow(abricate_report_row)
gene_detection_status_writer.writerow(gene_detection_status)
- f1.seek(0) # return file pointer to start of abricate report
+ f1.seek(0) # return file pointer to start of abricate report
if __name__ == '__main__':
parser = argparse.ArgumentParser()
parser.add_argument("abricate_report", help="Input: Abricate report to screen (tsv)")
parser.add_argument("--screening_file", help="Input: List of genes to screen for (tsv)")
- parser.add_argument("--screened_report", help="Output: Screened abricate report including only genes of interest (tsv)")
- parser.add_argument("--gene_detection_status", help="Output: detection status for all genes listed in the screening file (tsv)")
+ parser.add_argument("--screened_report", help=("Output: Screened abricate report ",
+ "including only genes of interest (tsv)"))
+ parser.add_argument("--gene_detection_status", help=("Output: detection status for all genes ",
+ "listed in the screening file (tsv)"))
args = parser.parse_args()
main(args)
diff -r b2d56a44a872 -r 4f963b583186 screen_abricate_report.xml
--- a/screen_abricate_report.xml Thu Oct 31 15:43:15 2019 -0400
+++ b/screen_abricate_report.xml Tue Nov 12 21:44:45 2019 -0500
@@ -18,8 +18,8 @@
-
-
+
+