Skip to main content

How to use a for loop to access inputs in snakemake?

I want to use the gatk MakeVcfSampleNameMap tool to generate a tab-delimited file that maps sample names to the corresponding VCF files in a Snakemake workflow. These vcfs are in a the directory input/

output - > (sample_name_map.txt)

input/vcf1.vcf.gz   control5-a
input/vcf2.vcf.gz   control5-c

MakeVcfSampleNameMap accepts a single --INPUT argument at a time, so I have to submit multiple VCF files to it in the form:

  • --INPUT {input.vcfs[0]}
  • --INPUT {input.vcfs[1]}, etc.

snakemake rule:

rule MakeVcfSampleNameMap:
    input:
        vcfs=["input/vcf1.vcf.gz", "input/vcf2.vcf.gz"],
    output:
        "sample_name_map.txt",
    shell:
        "gatk MakeVcfSampleNameMap \
        --INPUT {input.vcfs[0]} \
        --INPUT {input.vcfs[1]} \
        --OUTPUT {output}"

This becomes cumbersome when there are many input files. How can use a loop in the Snakemake rule to access all the VCF files?



source https://stackoverflow.com/questions/75700654/how-to-use-a-for-loop-to-access-inputs-in-snakemake

Comments

Popular posts from this blog

Why is my reports service not connecting?

I am trying to pull some data from a Postgres database using Node.js and node-postures but I can't figure out why my service isn't connecting. my routes/index.js file: const express = require('express'); const router = express.Router(); const ordersCountController = require('../controllers/ordersCountController'); const ordersController = require('../controllers/ordersController'); const weeklyReportsController = require('../controllers/weeklyReportsController'); router.get('/orders_count', ordersCountController); router.get('/orders', ordersController); router.get('/weekly_reports', weeklyReportsController); module.exports = router; My controllers/weeklyReportsController.js file: const weeklyReportsService = require('../services/weeklyReportsService'); const weeklyReportsController = async (req, res) => { try { const data = await weeklyReportsService; res.json({data}) console...

How to show number of registered users in Laravel based on usertype?

i'm trying to display data from the database in the admin dashboard i used this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count(); echo $users; ?> and i have successfully get the correct data from the database but what if i want to display a specific data for example in this user table there is "usertype" that specify if the user is normal user or admin i want to user the same code above but to display a specific usertype i tried this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count()->WHERE usertype =admin; echo $users; ?> but it didn't work, what am i doing wrong? source https://stackoverflow.com/questions/68199726/how-to-show-number-of-registered-users-in-laravel-based-on-usertype

The number of nodes and edges I calculated manually and by a method provided by networkx library dont match

I am trying to analyse a graph using networkx library in Python. I have implemented a function that returns the nodes and edges of that graph. After that I counted the number of nodes and edges by implementing the len() function for nodes and edges, which gave me 17 for both nodes and edges. However, when I used the method graph.number_of_nodes() of networkx library it gave me 53. My code: import pandas as pd import networkx as nx import matplotlib.pyplot as plt header_list = ["a","b","w"] E =pd.read_csv("aves-barn-swallow-contact-network.edges",sep=" ",header= None,names=header_list) # print(E.head()) # prints first few rows of the csv file G = nx.from_pandas_edgelist(E,"a","b",["w"]) nx.draw(G) plt.show() class GraphAnalyse: #finding nodes and edges def nodeEdge(self,Graph): nodes = Graph.nodes() edges = Graph.edges() ...