Mercurial > repos > iuc > nanocompore_sampcomp

annotate NanocomporeDB_process.py @ 4:89b4ed4145b3 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 51078f4c0eab9cc41a0cb2190000e428322631a4
author iuc
date Thu, 23 Jun 2022 17:27:25 +0000
parents e7d7d87dc1a3
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
1 #!/usr/bin/env python
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
2
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
3 import argparse
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
4 import os
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
5
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
6 from nanocompore.SampCompDB import SampCompDB
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
7
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
8
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
9 def is_valid_file(file_name):
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
10 if os.path.isfile(file_name):
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
11 return os.path.abspath(file_name)
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
12 else:
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
13 raise FileNotFoundError(os.path.abspath(file_name))
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
14
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
15
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
16 def is_valid_directory(dir_name):
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
17 if os.path.isdir(dir_name):
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
18 return os.path.abspath(dir_name)
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
19 else:
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
20 raise NotADirectoryError(os.path.abspath(dir_name))
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
21
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
22
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
23 if __name__ == '__main__':
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
24
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
25 parser = argparse.ArgumentParser(
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
26 description='save nanocompre sampcomp \
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
27 results as interval outputs \
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
28 \nSample call: \"python Nannocompore-plot.py --db-path \
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
29 ./out_SampComp.db --ref-fasta ref.fa --annotation-bed annot.bed \
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
30 --out-dir ./plots/')
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
31
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
32 parser.add_argument('--ref-fasta', required=True, type=is_valid_file,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
33 help='The reference genome used for read alignment.')
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
34 parser.add_argument('--db-path', default="./out_SampComp.db", type=str,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
35 help='Path to the SampCompDB database path prefix.')
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
36 parser.add_argument('--annotation-bed', required=False, type=is_valid_file,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
37 help='BED file containing the annotation of the transcriptome used as reference when mapping')
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
38 parser.add_argument('--pvalue-types', type=str,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
39 default='GMM_logit_pvalue,KS_dwell_pvalue,KS_intensity_pvalue',
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
40 help='path to the annotations')
2
f657489960f7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit a57b5b935a80db12a0ed867b4219453a2854dd32"
iuc
parents: 1
diff changeset
41 parser.add_argument('--bedgraph', action="store_true",
1
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
42 help='write output in BEDGRAPH format instead of BED')
3
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
43 parser.add_argument('--pvalue-threshold', type=float, default=1.0,
1
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
44 help='Maximum reported p-value.')
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
45 parser.add_argument('--out-dir', default="./", type=is_valid_directory,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
46 help='path the plotting output directory.')
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
47
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
48 args = parser.parse_args()
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
49
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
50 db = SampCompDB(args.db_path, fasta_fn=args.ref_fasta,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
51 bed_fn=args.annotation_bed)
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
52 print(db)
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
53 print("DB read ids:", db.ref_id_list)
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
54
3
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
55 if args.bedgraph:
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
56 file_ext = 'bedgraph'
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
57 else:
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
58 file_ext = 'bed'
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
59
1
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
60 if args.annotation_bed:
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
61 for pt in args.pvalue_types.split(','):
3
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
62 out_path = '{}/{}.{}'.format(args.out_dir, pt, file_ext)
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
63 print("%s output for p-value type:" % out_path, pt)
e7d7d87dc1a3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit e0d4688a59e6eeba33adcfe803ac43d0bc2863e7"
iuc
parents: 2
diff changeset
64 db.save_to_bed(output_fn=out_path,
1
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
65 pvalue_field=pt, pvalue_thr=args.pvalue_threshold,
e5b88f0e7960 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/nanocompore commit 652a872061324ba1073bfa286777ffeefa352671"
iuc
parents:
diff changeset
66 bedgraph=args.bedgraph)