Mercurial > repos > cpt > cpt_intersect_adj
comparison intersect_and_adjacent.py @ 3:10000414e916 draft
planemo upload commit 28919dd79f328dd09dae92661dc19ffa5757a765
author | cpt |
---|---|
date | Mon, 12 Aug 2024 04:23:34 +0000 |
parents | 923d2528480d |
children | 3e3b5ba626b9 |
comparison
equal
deleted
inserted
replaced
2:0dbc7fbad0ab | 3:10000414e916 |
---|---|
266 type=int, | 266 type=int, |
267 default=50, | 267 default=50, |
268 help="Allows features this far away to still be considered 'adjacent'", | 268 help="Allows features this far away to still be considered 'adjacent'", |
269 ) | 269 ) |
270 parser.add_argument( | 270 parser.add_argument( |
271 "-stranding", | 271 "stranding", |
272 action="store_true", | 272 nargs="?", |
273 help="Only allow adjacency for same-strand features", | 273 default="", |
274 help="Only allow adjacency for same-strand features if set to '-stranding'", | |
274 ) | 275 ) |
275 parser.add_argument("--oa", type=str, default="a_hits_near_b.gff") | 276 parser.add_argument("--oa", type=str, default="a_hits_near_b.gff") |
276 parser.add_argument("--ob", type=str, default="b_hits_near_a.gff") | 277 parser.add_argument("--ob", type=str, default="b_hits_near_a.gff") |
277 args = parser.parse_args() | 278 args = parser.parse_args() |
278 | 279 |
280 stranding = args.stranding == "-stranding" | |
279 b, a = intersect(args.a, args.b, args.window, args.stranding) | 281 b, a = intersect(args.a, args.b, args.window, args.stranding) |
280 | 282 |
281 with open(args.oa, "w") as handle: | 283 with open(args.oa, "w") as handle: |
282 for rec in a: | 284 for rec in a: |
283 gffWrite([rec], handle) | 285 gffWrite([rec], handle) |