From 33dbfce0bbf294afaec94e51ec35dc927b814da0 Mon Sep 17 00:00:00 2001 From: sophie22 Date: Sat, 7 May 2022 19:18:19 +0100 Subject: [PATCH] add context to output file --- genes_coverage.py | 2 +- genes_suboptimal_coverage30x.txt | 1 + 2 files changed, 2 insertions(+), 1 deletion(-) diff --git a/genes_coverage.py b/genes_coverage.py index 78edc49..6d45521 100644 --- a/genes_coverage.py +++ b/genes_coverage.py @@ -2,7 +2,6 @@ # Import libraries/packages for use in the code import sys -from numpy import outer import pandas as pd # v1.3.4 ### Read inputs from the command line @@ -30,5 +29,6 @@ ### Write gene symbols with suboptimal coverage to file outfile = f"genes_suboptimal_coverage{coverage_threshold}x.txt" with open(outfile, 'w') as fh: + fh.write(f"The following genes had regions not fully covered at 30x: \n") for gene in below_threshold_genes: fh.write(gene + "\n") diff --git a/genes_suboptimal_coverage30x.txt b/genes_suboptimal_coverage30x.txt index b7deb6d..f2f8704 100644 --- a/genes_suboptimal_coverage30x.txt +++ b/genes_suboptimal_coverage30x.txt @@ -1,3 +1,4 @@ +The following genes had regions not fully covered at 30x: SELENON MSTO1 NEB