2010-03-27 02:55:12 +08:00
|
|
|
package org.broadinstitute.sting.gatk.walkers;
|
2009-09-24 23:16:11 +08:00
|
|
|
|
|
|
|
|
import org.broadinstitute.sting.WalkerTest;
|
|
|
|
|
import org.junit.Test;
|
|
|
|
|
|
|
|
|
|
import java.util.List;
|
|
|
|
|
import java.util.ArrayList;
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* @author aaron
|
|
|
|
|
* <p/>
|
|
|
|
|
* Class VariantsToVCFIntegrationTest
|
|
|
|
|
* <p/>
|
|
|
|
|
* test(s) for the VariantsToVCF walker.
|
|
|
|
|
*/
|
|
|
|
|
public class VariantsToVCFIntegrationTest extends WalkerTest {
|
|
|
|
|
|
|
|
|
|
|
2010-03-25 12:53:31 +08:00
|
|
|
@Test
|
2009-09-24 23:16:11 +08:00
|
|
|
public void testVariantsToVCFUsingGeliInput() {
|
|
|
|
|
List<String> md5 = new ArrayList<String>();
|
2010-03-27 02:34:59 +08:00
|
|
|
md5.add("211be63cf93cddf021e5b0eb9341b386");
|
2009-09-24 23:16:11 +08:00
|
|
|
|
|
|
|
|
WalkerTest.WalkerTestSpec spec = new WalkerTest.WalkerTestSpec(
|
2009-12-29 15:40:48 +08:00
|
|
|
"-R " + oneKGLocation + "reference/human_b36_both.fasta" +
|
2010-03-25 12:53:31 +08:00
|
|
|
" -B variant,Variants," + validationDataLocation + "NA12878.1kg.p2.chr1_10mb_11_mb.SLX.lod5.variants.geli.calls" +
|
2009-09-24 23:16:11 +08:00
|
|
|
" -T VariantsToVCF" +
|
|
|
|
|
" -L 1:10,000,000-11,000,000" +
|
2010-03-25 12:53:31 +08:00
|
|
|
" -sample NA123AB" +
|
|
|
|
|
" -o %s",
|
2009-09-24 23:16:11 +08:00
|
|
|
1, // just one output file
|
|
|
|
|
md5);
|
2010-03-25 12:53:31 +08:00
|
|
|
executeTest("testVariantsToVCFUsingGeliInput #1", spec).getFirst();
|
2009-09-24 23:16:11 +08:00
|
|
|
}
|
|
|
|
|
|
2010-03-25 12:53:31 +08:00
|
|
|
@Test
|
2009-09-24 23:16:11 +08:00
|
|
|
public void testGenotypesToVCFUsingGeliInput() {
|
|
|
|
|
List<String> md5 = new ArrayList<String>();
|
2010-03-27 02:34:59 +08:00
|
|
|
md5.add("b31446fa91b8ed82ad73a5a5a72700a7");
|
2009-09-24 23:16:11 +08:00
|
|
|
|
|
|
|
|
WalkerTest.WalkerTestSpec spec = new WalkerTest.WalkerTestSpec(
|
2009-12-29 15:40:48 +08:00
|
|
|
"-R " + oneKGLocation + "reference/human_b36_both.fasta" +
|
2010-03-25 12:53:31 +08:00
|
|
|
" -B variant,Variants," + validationDataLocation + "NA12878.1kg.p2.chr1_10mb_11_mb.SLX.lod5.genotypes.geli.calls" +
|
2009-09-24 23:16:11 +08:00
|
|
|
" -T VariantsToVCF" +
|
|
|
|
|
" -L 1:10,000,000-11,000,000" +
|
2010-03-25 12:53:31 +08:00
|
|
|
" -sample NA123AB" +
|
|
|
|
|
" -o %s",
|
2009-09-24 23:16:11 +08:00
|
|
|
1, // just one output file
|
|
|
|
|
md5);
|
2010-03-25 12:53:31 +08:00
|
|
|
executeTest("testVariantsToVCFUsingGeliInput #2", spec).getFirst();
|
2009-09-24 23:16:11 +08:00
|
|
|
}
|
|
|
|
|
|
2010-03-27 02:34:59 +08:00
|
|
|
@Test
|
|
|
|
|
public void testGenotypesToVCFUsingHapMapInput() {
|
|
|
|
|
List<String> md5 = new ArrayList<String>();
|
|
|
|
|
md5.add("03ff126faf5751a83bd7ab9e020bce7e");
|
|
|
|
|
|
|
|
|
|
WalkerTest.WalkerTestSpec spec = new WalkerTest.WalkerTestSpec(
|
|
|
|
|
"-R " + oneKGLocation + "reference/human_b36_both.fasta" +
|
2010-04-11 09:43:07 +08:00
|
|
|
" -B variant,HapMap," + validationDataLocation + "rawHapMap.yri.chr1.txt" +
|
2010-03-27 02:34:59 +08:00
|
|
|
" -T VariantsToVCF" +
|
|
|
|
|
" -L 1:1-1,000,000" +
|
|
|
|
|
" -o %s",
|
|
|
|
|
1, // just one output file
|
|
|
|
|
md5);
|
|
|
|
|
executeTest("testVariantsToVCFUsingHapMapInput", spec).getFirst();
|
|
|
|
|
}
|
2009-09-24 23:16:11 +08:00
|
|
|
}
|