diff --git a/workflow/scripts/map_reads.py b/workflow/scripts/map_reads.py index a3a6abf85410682627095b02c995f34241387a31..9c09b0a9254ad0749cd36356212dd494b3162bf7 100644 --- a/workflow/scripts/map_reads.py +++ b/workflow/scripts/map_reads.py @@ -218,6 +218,13 @@ def main(): shlex.split("samtools flagstat %s" % (bam_filename)), stdout=temp_file) + #Genome/Bad fastq File Check + file_check = open(bam_mapstats_filename).readlines() + percent = file_check[4].split('(')[1] + percent = percent.split('%')[0] + if float(percent) < 10: + raise Exception ('Mapped Genes too low: Check for correct Genotype') + # Remove sai files for sai_file in sai: os.remove(sai_file)