本文整理匯總了Python中CGATPipelines.PipelineMappingQC.buildPicardAlignmentStats方法的典型用法代碼示例。如果您正苦於以下問題：Python PipelineMappingQC.buildPicardAlignmentStats方法的具體用法？Python PipelineMappingQC.buildPicardAlignmentStats怎麽用？Python PipelineMappingQC.buildPicardAlignmentStats使用的例子？那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在類CGATPipelines.PipelineMappingQC的用法示例。
在下文中一共展示了PipelineMappingQC.buildPicardAlignmentStats方法的4個代碼示例，這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚，您的評價將有助於係統推薦出更棒的Python代碼示例。
示例1: buildPicardStats
# 需要導入模塊: from CGATPipelines import PipelineMappingQC [as 別名]
# 或者: from CGATPipelines.PipelineMappingQC import buildPicardAlignmentStats [as 別名]
def buildPicardStats(infile, outfile):
'''build alignment stats using picard.
Note that picards counts reads but they are in fact alignments.
'''
if PARAMS["pool_reads"]:
reffile = os.path.join(
os.path.dirname(infile), "agg-agg-agg.filtered.contigs.fa")
else:
reffile = P.snip(infile, ".bam") + ".fa"
PipelineMappingQC.buildPicardAlignmentStats(infile,
outfile,
reffile)示例2: runPicardOnRealigned
# 需要導入模塊: from CGATPipelines import PipelineMappingQC [as 別名]
# 或者: from CGATPipelines.PipelineMappingQC import buildPicardAlignmentStats [as 別名]
def runPicardOnRealigned(infile, outfile):
to_cluster = USECLUSTER
job_options = getGATKOptions()
# TS no multithreading so why 6 threads?
# job_threads = 6
tmpdir_gatk = P.getTempDir('/ifs/scratch')
# threads = PARAMS["gatk_threads"]
outfile_tumor = outfile.replace("Control", PARAMS["mutect_tumour"])
infile_tumor = infile.replace("Control", PARAMS["mutect_tumour"])
track = P.snip(os.path.basename(infile), ".bam")
track_tumor = track.replace("Control", PARAMS["mutect_tumour"])
genome = "%s/%s.fa" % (PARAMS["bwa_index_dir"],
PARAMS["genome"])
PipelineMappingQC.buildPicardAlignmentStats(infile, outfile, genome)
PipelineMappingQC.buildPicardAlignmentStats(infile_tumor,
outfile_tumor, genome)
# check above functions then remove statement
statement = '''
cat %(infile)s
| python %%(scriptsdir)s/bam2bam.py -v 0 --method=set-sequence
| CollectMultipleMetrics
INPUT=/dev/stdin
REFERENCE_SEQUENCE=%%(bwa_index_dir)s/%%(genome)s.fa
ASSUME_SORTED=true
OUTPUT=%(outfile)s
VALIDATION_STRINGENCY=SILENT
>& %(outfile)s;
cat %(infile_tumor)s
| python %%(scriptsdir)s/bam2bam.py -v 0
--method=set-sequence --output-sam
| CollectMultipleMetrics
INPUT=/dev/stdin
REFERENCE_SEQUENCE=%%(bwa_index_dir)s/%%(genome)s.fa
ASSUME_SORTED=true
OUTPUT=%(outfile_tumor)s
VALIDATION_STRINGENCY=SILENT
>& %(outfile_tumor)s;''' % locals()示例3: runPicardOnRealigned
# 需要導入模塊: from CGATPipelines import PipelineMappingQC [as 別名]
# 或者: from CGATPipelines.PipelineMappingQC import buildPicardAlignmentStats [as 別名]
def runPicardOnRealigned(infile, outfile):
to_cluster = USECLUSTER
job_memory = PARAMS["gatk_memory"]
tmpdir_gatk = P.getTempDir()
outfile_tumor = outfile.replace(
PARAMS["sample_control"], PARAMS["sample_tumour"])
infile_tumor = infile.replace(
PARAMS["sample_control"], PARAMS["sample_tumour"])
track = P.snip(os.path.basename(infile), ".bam")
track_tumor = track.replace(
PARAMS["sample_control"], PARAMS["sample_tumour"])
genome = "%s/%s.fa" % (PARAMS["bwa_index_dir"],
PARAMS["genome"])
PipelineMappingQC.buildPicardAlignmentStats(infile, outfile, genome)
PipelineMappingQC.buildPicardAlignmentStats(infile_tumor,
outfile_tumor, genome)示例4: buildPicardAlignmentStats
# 需要導入模塊: from CGATPipelines import PipelineMappingQC [as 別名]
# 或者: from CGATPipelines.PipelineMappingQC import buildPicardAlignmentStats [as 別名]
def buildPicardAlignmentStats(infile, outfile):
'''Gather BAM file alignment statistics using Picard '''
PipelineMappingQC.buildPicardAlignmentStats(infile, outfile,
os.path.join(PARAMS["bwa_index_dir"],
PARAMS["genome"] + ".fa"))