This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SnpSift case control command for multiple samples

We are trying to write a script which can analyse multiple samples for SnpSift CaseControl. There is "+" symbol in the command to represent the no. of samples. How to define this option of the command in the script. Here is the script, we are trying to run, but it is not working.

!/bin/bash

for f1 in *vcf do java -jar SnpSift.jar caseControl "$f1+" $f1 > $f1_CC.vcf done

Please help with this.

snp next-gen rna-seq snpsift casecontrol
• 1,545 views • link
written by kumari.indu31 • kumari.indu31

Just a few remarks:

  • if you use for loop in one line, the syntax should be different and contain semicolons. Otherwise, you could use indentations.
  • if I am not mistaken, the usual input for caseControl is a vcf file with multiple samples in it (meaning multiple sample columns in the vcf file). In this case, the symbols +, -, 0 specify whether the case is cancer, control, or neutral.
• 0 views • link
written by e.rempel e.rempel

My question is regarding how can I define no. of samples in the command in script so that there is no need to add "+" symbol manually.

• 0 views • link
written by kumari.indu31 • kumari.indu31

You can't run caseControl without a multi sample vcf.

• 0 views • link
written by Arup Ghosh Arup Ghosh

0 answers

No answers yet.

Log in to answer this question.

More posts like this