Skip to content

Commit

Permalink
chages to exome.py
Browse files Browse the repository at this point in the history
  • Loading branch information
jscaber committed Sep 20, 2023
1 parent af5661a commit a0524b9
Showing 1 changed file with 0 additions and 15 deletions.
15 changes: 0 additions & 15 deletions cgatpipelines/tasks/exome.py
Original file line number Diff line number Diff line change
Expand Up @@ -337,21 +337,6 @@ def vcfToTable(infile, outfile, genome, columns, gatkmem):
##############################################################################


def snpSift(infile, outfile, vcf, bgzip=True, memory="6G"):
'''Annotate VCF with other VCF'''
job_memory=memory
statement = '''SnpSift annotate
-Xmx%(job_memory)s %(vcf)s
%(infile)s > %(outfile)s 2> %(outfile)s.log;'''
if bgzip:
statement += '''bgzip %(outfile)s;
tabix -p vcf %(outfile)s.gz'''
P.run(statement)


##############################################################################


def selectVariants(infile, outfile, genome, select):
'''Filter de novo variants based on provided expression'''
statement = '''gatk SelectVariants
Expand Down

0 comments on commit a0524b9

Please sign in to comment.