本文整理匯總了Python中CGATPipelines.PipelineMappingQC.buildPicardCoverageStats方法的典型用法代碼示例。如果您正苦於以下問題：Python PipelineMappingQC.buildPicardCoverageStats方法的具體用法？Python PipelineMappingQC.buildPicardCoverageStats怎麽用？Python PipelineMappingQC.buildPicardCoverageStats使用的例子？那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在類CGATPipelines.PipelineMappingQC的用法示例。
在下文中一共展示了PipelineMappingQC.buildPicardCoverageStats方法的1個代碼示例，這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚，您的評價將有助於係統推薦出更棒的Python代碼示例。
示例1: buildCoverageStats
# 需要導入模塊: from CGATPipelines import PipelineMappingQC [as 別名]
# 或者: from CGATPipelines.PipelineMappingQC import buildPicardCoverageStats [as 別名]
def buildCoverageStats(infile, outfile):
'''Generate coverage statistics for regions of interest from a
bed file using Picard'''
# TS check whether this is always required or specific to current baits file
# baits file requires modification to make picard accept it
# this is performed before CalculateHsMetrics
to_cluster = USECLUSTER
baits = PARAMS["roi_baits"]
modified_baits = infile + "_temp_baits_final.bed"
regions = PARAMS["roi_regions"]
statement = '''samtools view -H %(infile)s > %(infile)s_temp_header.txt;
awk 'NR>2' %(baits)s |
awk -F '\\t' 'BEGIN { OFS="\\t" } {print $1,$2,$3,"+",$4;}'
> %(infile)s_temp_baits.bed;
cat %(infile)s_temp_header.txt %(infile)s_temp_baits.bed
> %(modified_baits)s; checkpoint ;
rm -rf %(infile)s_temp_baits.bed %(infile)s_temp_header.txt
'''
P.run()
PipelineMappingQC.buildPicardCoverageStats(
infile, outfile, modified_baits, modified_baits)
IOTools.zapFile(modified_baits)