def main():
args = docopt(__doc__)
bam_f = args['--bam']
include_f = args['--include']
exclude_f = args['--exclude']
out_prefix = args['--out']
gzip = args['--gzip']
do_sort = args['--sort']
keep_sorted = args['--keep']
sort_threads = int(args['--threads'])
print BtLog.status_d['22'] % bam_f
out_f = BtIO.getOutFile(bam_f, out_prefix, None)
if include_f and exclude_f:
print BtLog.error('43')
elif include_f:
sequence_list = BtIO.parseList(include_f)
BtIO.parseBamForFilter(bam_f, out_f, sequence_list, None, gzip,
do_sort, keep_sorted, sort_threads)
elif exclude_f:
sequence_list = BtIO.parseList(exclude_f)
BtIO.parseBamForFilter(bam_f, out_f, None, sequence_list, gzip,
do_sort, keep_sorted, sort_threads)
else:
BtIO.parseBamForFilter(bam_f, out_f, None, None, gzip, do_sort,
keep_sorted, sort_threads)
def main():
args = docopt(__doc__)
bam_f = args['--bam']
include_f = args['--include']
exclude_f = args['--exclude']
out_prefix = args['--out']
gzip = args['--gzip']
do_sort = args['--sort']
keep_sorted = args['--keep']
sort_threads = int(args['--threads'])
print BtLog.status_d['22'] % bam_f
out_f = BtIO.getOutFile(bam_f, out_prefix, None)
if include_f and exclude_f:
print BtLog.error('43')
elif include_f:
sequence_list = BtIO.parseList(include_f)
BtIO.parseBamForFilter(bam_f, out_f, sequence_list, None, gzip, do_sort, keep_sorted, sort_threads)
elif exclude_f:
sequence_list = BtIO.parseList(exclude_f)
BtIO.parseBamForFilter(bam_f, out_f, None, sequence_list, gzip, do_sort, keep_sorted, sort_threads)
else:
BtIO.parseBamForFilter(bam_f, out_f, None, None, gzip, do_sort, keep_sorted, sort_threads)
def main():
#print data_dir
args = docopt(__doc__)
blobdb_f = args['--input']
prefix = args['--out']
ranks = args['--rank']
taxrule = args['--taxrule']
hits_flag = args['--hits']
seq_list_f = args['--list']
concoct = args['--concoct']
cov = args['--cov']
notable = args['--notable']
experimental = args['--experimental']
# Does blobdb_f exist ?
if not isfile(blobdb_f):
BtLog.error('0', blobdb_f)
out_f = BtIO.getOutFile(blobdb_f, prefix, None)
# Are ranks sane ?
if 'all' in ranks:
temp_ranks = RANKS[0:-1]
ranks = temp_ranks[::-1]
else:
for rank in ranks:
if rank not in RANKS:
BtLog.error('9', rank)
# Does seq_list file exist?
seqs = []
if (seq_list_f):
if isfile(seq_list_f):
seqs = BtIO.parseList(seq_list_f)
else:
BtLog.error('0', seq_list_f)
# Load BlobDb
blobDb = BtCore.BlobDb('new')
print BtLog.status_d['9'] % (blobdb_f)
blobDb.load(blobdb_f)
blobDb.version = blobtools.__version__
# Is taxrule sane and was it computed?
if (blobDb.hitLibs) and taxrule not in blobDb.taxrules:
BtLog.error('11', taxrule, blobDb.taxrules)
# view(s)
viewObjs = []
print BtLog.status_d['14']
if not (notable):
tableView = BtCore.ViewObj(name="table", out_f=out_f, suffix="table.txt", body=[])
viewObjs.append(tableView)
if (experimental):
experimentalView = BtCore.ExperimentalViewObj(name = "experimental", view_dir=out_f)
viewObjs.append(experimentalView)
if (concoct):
concoctTaxView = BtCore.ViewObj(name="concoct_tax", out_f=out_f, suffix="concoct_taxonomy_info.csv", body=dict())
viewObjs.append(concoctTaxView)
concoctCovView = BtCore.ViewObj(name="concoct_cov", out_f=out_f, suffix="concoct_coverage_info.tsv", body=[])
viewObjs.append(concoctCovView)
if (cov):
for cov_lib_name, covLibDict in blobDb.covLibs.items():
out_f = BtIO.getOutFile(covLibDict['f'], prefix, None)
covView = BtCore.ViewObj(name="covlib", out_f=out_f, suffix="cov", body=[])
blobDb.view(viewObjs=[covView], ranks=None, taxrule=None, hits_flag=None, seqs=None, cov_libs=[cov_lib_name], progressbar=True)
if (viewObjs):
blobDb.view(viewObjs=viewObjs, ranks=ranks, taxrule=taxrule, hits_flag=hits_flag, seqs=seqs, cov_libs=[], progressbar=True)
print BtLog.status_d['19']
def main():
#print data_dir
args = docopt(__doc__)
blobdb_f = args['--input']
prefix = args['--out']
ranks = args['--rank']
taxrule = args['--taxrule']
hits_flag = args['--hits']
seq_list_f = args['--list']
concoct = args['--concoct']
cov = args['--cov']
notable = args['--notable']
experimental = args['--experimental']
# Does blobdb_f exist ?
if not isfile(blobdb_f):
BtLog.error('0', blobdb_f)
out_f = BtIO.getOutFile(blobdb_f, prefix, None)
# Are ranks sane ?
if 'all' in ranks:
temp_ranks = RANKS[0:-1]
ranks = temp_ranks[::-1]
else:
for rank in ranks:
if rank not in RANKS:
BtLog.error('9', rank)
# Does seq_list file exist?
seqs = []
if (seq_list_f):
if isfile(seq_list_f):
seqs = BtIO.parseList(seq_list_f)
else:
BtLog.error('0', seq_list_f)
# Load BlobDb
blobDb = BtCore.BlobDb('new')
print BtLog.status_d['9'] % (blobdb_f)
blobDb.load(blobdb_f)
blobDb.version = blobtools.__version__
# Is taxrule sane and was it computed?
if (blobDb.hitLibs) and taxrule not in blobDb.taxrules:
BtLog.error('11', taxrule, blobDb.taxrules)
# view(s)
viewObjs = []
print BtLog.status_d['14']
if not (notable):
tableView = BtCore.ViewObj(name="table",
out_f=out_f,
suffix="table.txt",
body=[])
viewObjs.append(tableView)
if (experimental):
experimentalView = BtCore.ExperimentalViewObj(name="experimental",
view_dir=out_f)
viewObjs.append(experimentalView)
if (concoct):
concoctTaxView = BtCore.ViewObj(name="concoct_tax",
out_f=out_f,
suffix="concoct_taxonomy_info.csv",
body=dict())
viewObjs.append(concoctTaxView)
concoctCovView = BtCore.ViewObj(name="concoct_cov",
out_f=out_f,
suffix="concoct_coverage_info.tsv",
body=[])
viewObjs.append(concoctCovView)
if (cov):
for cov_lib_name, covLibDict in blobDb.covLibs.items():
out_f = BtIO.getOutFile(covLibDict['f'], prefix, None)
covView = BtCore.ViewObj(name="covlib",
out_f=out_f,
suffix="cov",
body=[])
blobDb.view(viewObjs=[covView],
ranks=None,
taxrule=None,
hits_flag=None,
seqs=None,
cov_libs=[cov_lib_name],
progressbar=True)
if (viewObjs):
blobDb.view(viewObjs=viewObjs,
ranks=ranks,
taxrule=taxrule,
hits_flag=hits_flag,
seqs=seqs,
cov_libs=[],
progressbar=True)
print BtLog.status_d['19']
