Mercurial > repos > iuc > bcftools_plugin_color_chrs
changeset 8:47ee1925fe9e draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bcftools commit e648d86d550ddf2eb67237752320c390b3a780e5
author | iuc |
---|---|
date | Wed, 05 Jun 2019 12:59:21 -0400 |
parents | 78a435cb9900 |
children | 05928e15da21 |
files | bcftools_plugin_color_chrs.xml macros.xml test-data/cnv.vcf test-data/cnv_baf_only.vcf test-data/cnv_pairwise_summary.tab test-data/cnv_summary.tab test-data/color_chrs.dat test-data/color_chrs_unrelated.dat test-data/fasta_indexes.loc test-data/mpileup.ref.fa.fai test-data/plugin_frameshifts.bed test-data/stats.a.vcf test-data/test-cache/23andme.fa test-data/test-cache/23andme.fa.fai test-data/test-cache/consensus.fa test-data/test-cache/consensus.fa.fai test-data/test-cache/csq.fa test-data/test-cache/csq.fa.fai test-data/test-cache/gvcf.fa test-data/test-cache/gvcf.fa.fai test-data/test-cache/mpileup.ref.fa test-data/test-cache/mpileup.ref.fa.fai test-data/test-cache/norm.fa test-data/test-cache/norm.fa.fai tool_data_table_conf.xml.test |
diffstat | 25 files changed, 840 insertions(+), 103 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bcftools_plugin_color_chrs.xml Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,105 @@ +<?xml version='1.0' encoding='utf-8'?> +<tool name="bcftools @EXECUTABLE@" id="bcftools_plugin_@PLUGIN_ID@" version="@TOOL_VERSION@+galaxy1"> + <description>plugin Color shared chromosomal segments, requires phased GTs</description> + <macros> + <token name="@EXECUTABLE@">color-chrs</token> + <token name="@PLUGIN_ID@">color_chrs</token> + <import>macros.xml</import> + </macros> + <expand macro="requirements" /> + <expand macro="version_command" /> + <command detect_errors="aggressive"><![CDATA[ +@PREPARE_ENV@ +@PREPARE_INPUT_FILE@ +#set $section = $sec_restrict +@PREPARE_TARGETS_FILE@ +@PREPARE_REGIONS_FILE@ + +bcftools plugin @EXECUTABLE@ + +## VCF input section +#set $section = $sec_restrict +@REGIONS@ +@TARGETS@ +@INCLUDE@ +@EXCLUDE@ + +@THREADS@ + +## Primary Input/Outputs +@INPUT_FILE@ + +## Plugin options +#if $sample_relation.sample_rel_sel == 'trio': +-- --trio '${sample_relation.mother},${sample_relation.father},${sample_relation.child}' +#else: +-- --unrelated '${sample_relation.sample_a},${sample_relation.sample_b}' +#end if +-p output_tmp +&& color-chrs.pl output_tmp.dat -p output_tmp +&& mv output_tmp.dat '$output_data' +&& mv output_tmp.svg '$output_svg' +]]> + </command> + <inputs> + <expand macro="macro_input" /> + <conditional name="sample_relation"> + <param name="sample_rel_sel" type="select" label="Sample relationship"> + <option value="trio">trio - mother,father,child</option> + <option value="unrelated">unrelated sample pair</option> + </param> + <when value="trio"> + <param name="mother" type="text" value="" label="name of mother"/> + <param name="father" type="text" value="" label="name of father"/> + <param name="child" type="text" value="" label="name of child"/> + </when> + <when value="unrelated"> + <param name="sample_a" type="text" value="" label="name of sample A"/> + <param name="sample_b" type="text" value="" label="name of sample B"/> + </when> + </conditional> + <section name="sec_restrict" expanded="false" title="Restrict analysis to"> + <expand macro="macro_include" /> + <expand macro="macro_exclude" /> + <expand macro="macro_restrict" /> + <expand macro="macro_restrict" type="target" label_type="Target" /> + </section> + </inputs> + <outputs> + <data name="output_data" format="tabular" /> + <data name="output_svg" format="svg" /> + </outputs> + <tests> + <test expect_num_outputs="2"> + <param name="input_file" ftype="vcf" value="stats.a.vcf" /> + <conditional name="sample_relation"> + <param name="sample_rel_sel" value="trio" /> + <param name="mother" value="A" /> + <param name="father" value="C" /> + <param name="child" value="B" /> + </conditional> + <output name="output_data" ftype="tabular" file="color_chrs.dat" /> + </test> + <test expect_num_outputs="2"> + <param name="input_file" ftype="vcf" value="stats.a.vcf" /> + <conditional name="sample_relation"> + <param name="sample_rel_sel" value="unrelated" /> + <param name="sample_a" value="A" /> + <param name="sample_b" value="B" /> + </conditional> + <output name="output_data" ftype="tabular" file="color_chrs_unrelated.dat" /> + </test> + </tests> + <help><![CDATA[ +===================================== + bcftools @EXECUTABLE@ plugin +===================================== + +@REGIONS_HELP@ +@TARGETS_HELP@ +@EXPRESSIONS_HELP@ + +]]> +</help> + <expand macro="citations" /> +</tool>
--- a/macros.xml Thu Feb 21 15:55:03 2019 -0500 +++ b/macros.xml Wed Jun 05 12:59:21 2019 -0400 @@ -1,13 +1,5 @@ <macros> <token name="@TOOL_VERSION@">1.9</token> - <xml name="stdio"> - <stdio> - <exit_code range="1:" /> - <exit_code range=":-1" /> - <regex match="Error:" /> - <regex match="Exception:" /> - </stdio> - </xml> <xml name="requirements"> <requirements> <requirement type="package" version="@TOOL_VERSION@">bcftools</requirement> @@ -18,6 +10,9 @@ <xml name="samtools_requirement"> <requirement type="package" version="1.9">samtools</requirement> </xml> + <xml name="matplotlib_requirement"> + <requirement type="package" version="3.1.0">matplotlib</requirement> + </xml> <xml name="version_command"> <version_command>bcftools 2>&1 | grep 'Version:'</version_command> </xml> @@ -111,16 +106,45 @@ $vcfs_list_file </token> + <xml name="test_using_reference" token_select_from="history" token_ref=""> + <conditional name="reference_source"> + <param name="reference_source_selector" value="@SELECT_FROM@" /> + <param name="fasta_ref" ftype="fasta" value="@REF@" /> + </conditional> + </xml> + <xml name="macro_fasta_ref"> - <param name="fasta_ref" argument="--fasta-ref" type="data" format="data" optional="true" label="Reference sequence in FASTA format" /> + <conditional name="reference_source"> + <param name="reference_source_selector" type="select" label="Choose the source for the reference genome"> + <option value="cached">Use a built-in genome</option> + <option value="history">Use a genome from the history</option> + </param> + <when value="cached"> + <param name="fasta_ref" type="select" label="Reference genome"> + <options from_data_table="fasta_indexes"> + <filter type="data_meta" column="1" key="dbkey" ref="input_file" /> + <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file" /> + </options> + </param> + </when> + <when value="history"> + <param name="fasta_ref" type="data" format="fasta" label="Reference genome" /> + </when> + </conditional> </xml> <token name="@PREPARE_FASTA_REF@"> <![CDATA[ #set $input_fa_ref = None #if 'fasta_ref' in $section and $section.fasta_ref: - #set $input_fa_ref = 'ref.fa' - ln -s '$section.fasta_ref' $input_fa_ref && - samtools faidx $input_fa_ref && + #if 'reference_source_selector' in $section: + #if str($section.reference_source_selector) == "history": + #set $input_fa_ref = 'ref.fa' + ln -s '$section.fasta_ref' $input_fa_ref && + samtools faidx $input_fa_ref && + #else: + #set $input_fa_ref = str($section.fasta_ref.fields.path) + #end if + #end if #end if ]]> </token> @@ -132,27 +156,6 @@ #end if </token> - <xml name="macro_ref_fasta"> - <conditional name="reference_source"> - <param name="reference_source_selector" type="select" label="Choose the source for the reference genome"> - <option value="cached">Locally cached</option> - <option value="history">History</option> - </param> - <when value="cached"> - <param name="ref_file" type="select" label="Select reference genome"> - <options from_data_table="fasta_indexes"> - <!--<filter type="data_meta" key="dbkey" ref="input_bam" column="value"/>--> - </options> - <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/> - </param> - </when> - <when value="history"> <!-- FIX ME!!!! --> - <param name="ref_file" type="data" format="fasta" label="Using reference file" /> - </when> - </conditional> - </xml> - - <xml name="macro_AF_file"> <param name="AF_file" argument="--AF-file" type="data" format="tabular" optional="true" label="Allele frequencies file" help="Tab-delimited file containing the columns CHR,POS,REF,ALT,AF" /> </xml> @@ -268,25 +271,97 @@ </data> </xml> - <xml name="macro_regions"> - <conditional name="regions"> - <param name="regions_src" type="select" label="Regions"> - <option value="__none__">None</option> - <option value="regions">regions</option> - <option value="regions_file">regions-file</option> + <xml name="macro_invert_targets"> + <param name="invert_targets_file" type="boolean" truevalue="^" falsevalue="" + label="Invert Targets" + help="inverts the query/filtering applied by the targets" /> + </xml> + + <xml name="macro_restriction_spec" token_type="region" token_label_type="Region"> + <repeat name="@TYPE@s" title="@LABEL_TYPE@ Filter" default="1" min="1"> + <param name="chrom" type="text" label="@LABEL_TYPE@ chromosome"> + <validator type="expression" message="A chromosome identifier is required when specifying a @LABEL_TYPE@ filter">value.strip()</validator> + </param> + <param name="start" type="text" label="@LABEL_TYPE@ start position"> + <validator type="expression" message="an integer number is required">not value or value.isdigit()</validator> + </param> + <param name="stop" type="text" label="@LABEL_TYPE@ end position"> + <validator type="expression" message="an integer number is required">not value or value.isdigit()</validator> + </param> + <yield /> + </repeat> + </xml> + + <xml name="macro_restrictions_file" token_type="region" token_label_type="Region"> + <param name="@TYPE@s_file" type="data" format="tabular" label="@LABEL_TYPE@s File" help="restrict to @LABEL_TYPE@s listed in a file" /> + </xml> + + <xml name="macro_restrict" token_type="region" token_label_type="Region" > + <conditional name="@TYPE@s"> + <param name="@TYPE@s_src" type="select" label="@LABEL_TYPE@s"> + <option value="__none__">Do not restrict to @LABEL_TYPE@s</option> + <option value="@TYPE@s">Specify one or more @LABEL_TYPE@(s) directly</option> + <option value="@TYPE@s_file">Operate on @LABEL_TYPE@s specified in a history dataset</option> </param> <when value="__none__"/> - <when value="regions"> - <param name="regions" type="text" value="" optional="true" label="Restrict to comma-separated list of regions" - help="Each region is specifed as: chr or chr:pos or chr:from-to"> - <validator type="regex" message="">^(\w+(:\d+(-\d+)?)?(,\w+(:\d+(-\d+)?)?)*)?$</validator> - </param> + <when value="@TYPE@s"> + <expand macro="macro_restriction_spec" type="@TYPE@" label_type="@LABEL_TYPE@" /> + <yield /> </when> - <when value="regions_file"> - <param name="regions_file" type="data" format="vcf,bed,tabular" optional="true" label="Regions file" help="Restrict to regions listed in a file" /> - </when> + <when value="@TYPE@s_file"> + <expand macro="macro_restrictions_file" type="@TYPE@" label_type="@LABEL_TYPE@" /> + <yield /> + </when> </conditional> </xml> + + <token name="@PARSE_INTERVALS@"> +<![CDATA[ +#set $components = [] +#for $i in $intervals: + #set $chrom = str($i.chrom).strip() + #set $start = str($i.start).strip() + #set $stop = str($i.stop).strip() + #if $start or $stop: + $components.append($chrom + ':' + ($start or '0') + '-' + $stop) + #else: + $components.append($chrom) + #end if +#end for +#set $intervals_spec = ','.join($components) +]]> + </token> + + <token name="@REGIONS@"> +<![CDATA[ +#if $section.regions.regions_src == 'regions': + #set $intervals = $section.regions.regions + @PARSE_INTERVALS@ + --regions '$intervals_spec' +#elif $section.regions.regions_src == 'regions_file' and $section.regions.regions_file: + #if $regions_path is not None: + --regions-file '$regions_path' + #else: + --regions-file '$section.regions.regions_file' + #end if +#end if +]]> + </token> + + <token name="@TARGETS@"> +<![CDATA[ +#if $targets_path: + --targets-file "${section.targets.invert_targets_file}${targets_path}" +#elif $section.targets.targets_src == 'targets': + #set $intervals = $section.targets.targets + @PARSE_INTERVALS@ + --targets '${section.targets.invert_targets_file}$intervals_spec' +#elif $section.targets.targets_src == 'targets_file' and $section.targets.targets_file: + --targets-file "${section.targets.invert_targets_file}${section.targets.targets_file}" +#end if +]]> + </token> + <token name="@PREPARE_REGIONS_FILE@"> <![CDATA[ #set $regions_path = None @@ -300,23 +375,7 @@ #end if ]]> </token> - <token name="@REGIONS@"> -#if $section.regions.regions_src == 'regions' and $section.regions.regions != '': - --regions '$section.regions.regions' -#elif $section.regions.regions_src == 'regions_file' and $section.regions.regions_file: - #if $regions_path is not None: - --regions-file '$regions_path' - #else: - --regions-file '$section.regions.regions_file' - #end if -#end if - </token> - <xml name="macro_targets_file"> - <param name="targets_file" type="data" format="tabular" label="Targets File" help="restrict to targets listed in a file" > - <yield/> - </param> - <param name="invert_targets_file" type="boolean" truevalue="^" falsevalue="" label="Invert Targets" help="inverts the query/filtering applied by the target file selection" /> - </xml> + <token name="@PREPARE_TARGETS_FILE@"> <![CDATA[ #set $targets_path = None @@ -333,6 +392,7 @@ #end if ]]> </token> + <token name="@TARGETS_FILE@"> <![CDATA[ #if $targets_path is not None: @@ -343,41 +403,6 @@ ]]> </token> - <xml name="macro_targets"> - <conditional name="targets"> - <param name="targets_src" type="select" label="Targets"> - <option value="__none__">None</option> - <option value="targets">targets</option> - <option value="targets_file">targets-file</option> - </param> - <when value="__none__"/> - <when value="targets"> - <param name="targets" type="text" value="" optional="true" label="Restrict to comma-separated list of targets" - help="Each target is specifed as: chr or chr:pos or chr:from-to"> - <validator type="regex" message="">^(\w+(:\d+(-\d+)?)?(,\w+(:\d+(-\d+)?)?)*)?$</validator> - </param> - <param name="invert_targets_file" type="boolean" truevalue="^" falsevalue="" label="Invert Targets" help="inverts the query/filtering applied by the targets" /> - </when> - <when value="targets_file"> - <expand macro="macro_targets_file"> - </expand> - </when> - </conditional> - </xml> - <token name="@TARGETS@"> -<![CDATA[ -#if $targets_path: - --targets-file "${section.targets.invert_targets_file}${targets_path}" -#else: - #if $section.targets.targets_src == 'targets' and $section.targets.targets != '': - --targets '${section.targets.invert_targets_file}${section.targets.targets}' - #elif $section.targets.targets_src == 'targets_file' and $section.targets.targets_file: - --targets-file "${section.targets.invert_targets_file}${section.targets.targets_file}" - #end if -#end if -]]> - </token> - <xml name="macro_samples"> <param argument="--samples" type="text" value="" optional="true" label="Samples" help="Comma separated list of samples to annotate (or exclude)">
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/cnv.vcf Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,192 @@ +##fileformat=VCFv4.2 +##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> +##FORMAT=<ID=GC,Number=1,Type=Float,Description="GenCall score"> +##FORMAT=<ID=IA,Number=1,Type=Float,Description="Intensity of the A Allele"> +##FORMAT=<ID=IB,Number=1,Type=Float,Description="Intensity of the B Allele"> +##FORMAT=<ID=BAF,Number=1,Type=Float,Description="B Allele Frequency"> +##FORMAT=<ID=LRR,Number=1,Type=Float,Description="Log R Ratio"> +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT test +10 135656 rs10904561 T G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.0415:0.5176:0.543:0.661:0.9477 +10 135708 rs7917054 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.2441:0.501:0.742:0.659:0.9604 +10 178434 rs7089889 T G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1440:0.4833:0.836:0.827:0.8884 +10 188805 rs12146291 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.2144:0.4715:0.776:0.743:0.9610 +10 203471 rs10903451 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:0.4201:0.0000:0.619:0.022:0.9337 +10 252693 rs2379078 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.2673:0.0068:0.806:0.081:0.9018 +10 273344 rs2448378 A C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.7407:0.0000:0.497:0.012:0.9411 +10 279248 rs2496278 C A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.4228:0.6268:0.285:0.958:0.6772 +10 291134 rs10508201 T C . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1085:0.4723:0.317:0.444:0.9106 +10 293358 exm2249181 A C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0097:0.0044:1.507:0.061:0.3992 +10 313504 rs4880568 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.5836:0.9803:0.017:0.321:0.9547 +10 317501 rs2018975 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.3392:0:1.069:0.039:0.7852 +10 323283 rs3740304 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1836:0.5588:1.053:1.110:0.8640 +10 326894 rs3125027 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0517:0.0005:1.300:0.026:0.8916 +10 327162 exm805022 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1197:0:0.995:0.015:0.5328 +10 329493 rs3815985 C A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0480:0.5187:0.824:0.794:0.9109 +10 338553 rs3125031 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.3387:0.4696:0.777:1.004:0.7901 +10 354301 rs35198327 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.1338:0.4929:0.903:0.996:0.8546 +10 375475 variant.11341 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0747:0.000199999999999978:1.132:0.012:0.4773 +10 377161 rs7070654 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1388:0:0.582:0.020:0.9458 +10 387060 rs7904155 G A . . . GT:LRR:BAF:IA:IB:GC 1/1:0.0615:0.9945:0.080:1.883:0.8310 +10 390962 exm805090 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0087:9.9999999999989e-05:1.226:0.063:0.6950 +10 390983 exm805091 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.3606:0:0.398:0.035:0.4585 +10 394615 exm805105 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0170:0:1.345:0.020:0.4108 +10 402390 exm805131 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0341:0.0000:1.127:0.023:0.8735 +10 403792 exm805136 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1428:0.0000:0.908:0.019:0.8483 +10 410501 exm805166 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0189:0:1.162:0.050:0.4628 +10 413010 rs9787422 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0871:0.5334:0.908:0.961:0.8338 +10 423555 rs12245012 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1264:0.00149999999999995:1.658:0.040:0.8015 +10 423671 rs10904067 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.3622:0.5566:0.326:0.370:0.9306 +10 431161 exm2271240 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.0175:0.4838:0.935:1.004:0.8308 +10 435971 exm2249233 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.2716:0.0000:1.217:0.042:0.3992 +10 435990 exm805214 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1427:0:1.107:0.031:0.4295 +10 453567 rs4881254 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.1831:0.5079:1.032:1.425:0.8149 +10 462885 rs10904173 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.0930:0.4819:0.568:0.489:0.9332 +10 465046 exm805254 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1108:0:1.041:0.038:0.8704 +10 468599 rs12415961 T C . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1647:0.5021:1.192:0.825:0.8051 +10 486935 exm805280 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1369:0.0186:0.643:0.060:0.4804 +10 487973 rs4881313 T C . . . GT:LRR:BAF:IA:IB:GC 0/1:0.0644:0.5628:0.747:0.992:0.8463 +10 521431 exm2271370 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:0.0436:0.9952:0.028:1.491:0.8545 +10 521723 rs4881336 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.2056:1:0.055:1.228:0.8344 +10 547567 rs816599 A G . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.0079:1.0000:0.023:1.101:0.9490 +10 554186 rs12251997 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1537:0.00309999999999999:1.565:0.032:0.8228 +10 556129 rs10466270 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.0014:0.5142:0.709:0.732:0.9257 +10 566379 rs11252926 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0211:0.5024:0.559:0.493:0.9438 +10 588406 rs10904450 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0834:0.0006:1.837:0.066:0.8107 +10 592283 rs816646 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.3907:0.5006:0.268:0.158:0.8540 +10 601089 rs816650 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0897:0:1.053:0.009:0.9557 +10 615782 rs17221309 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.2800:0.00319999999999998:1.764:0.035:0.7820 +10 635857 rs7069611 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0450:0.0013:0.593:0.017:0.9595 +10 669358 rs7898821 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1257:0:1.694:0.067:0.7270 +10 673037 rs816563 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:0.1368:0.9065:0.167:0.495:0.4895 +10 680428 rs2124585 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.0590:0.5212:0.660:1.123:0.7133 +10 712617 rs17136372 A C . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.0174:0.9992:0.028:0.975:0.9556 +10 738630 rs1750792 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0269:0:1.575:0.025:0.8126 +10 749238 rs553326 G T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1758:0.00229999999999997:0.955:0.019:0.9300 +10 754033 rs2265090 C A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1347:0:1.181:0.025:0.9176 +10 757222 rs10904546 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:0.3035:0.0000:1.975:0.056:0.8047 +10 766105 rs7906313 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1672:0.0116:0.960:0.062:0.9381 +10 776027 rs11253377 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.1474:0.9964:0.046:1.017:0.9087 +10 787896 rs1769242 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.2049:0.0000:0.965:0.029:0.9107 +10 800213 rs2254501 G A . . . GT:LRR:BAF:IA:IB:GC 1/1:0.2422:0.9918:0.027:0.612:0.9374 +10 802445 rs1769215 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0188:0.00370000000000004:0.859:0.024:0.9303 +10 803721 rs2790381 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.0421:0.9798:0.029:0.485:0.9700 +10 811876 rs12414585 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0717:0.516:0.374:0.276:0.9351 +10 813426 rs11253424 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.5256:0:0.333:0.000:0.8993 +10 815149 rs7084027 A G . . . GT:LRR:BAF:IA:IB:GC 1/1:0.3183:0.9944:0.037:1.547:0.8419 +10 820868 rs10752019 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0749:0.0118:1.349:0.070:0.8882 +10 833528 exm2271241 T G . . . GT:LRR:BAF:IA:IB:GC 0/1:-0.0329:0.4896:0.569:0.681:0.9348 +10 834736 rs4881530 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0453:0.00249999999999995:1.578:0.040:0.8076 +10 836115 rs11253444 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:-0.1380:0.5199:0.596:0.670:0.9248 +10 838179 rs1536337 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0534:0.0000:1.318:0.095:0.7967 +10 858022 rs9124 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.5565:0:0.715:0.012:0.9397 +10 858924 exm805329 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0866:0:0.873:0.021:0.9056 +10 859047 exm805338 C A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1868:0:0.917:0.007:0.5296 +10 860687 exm805350 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1509:0:1.172:0.059:0.3777 +10 860726 exm805354 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0224:0.0000:2.255:0.138:0.3747 +10 860970 exm805360 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1069:0.0000:1.859:0.070:0.7352 +10 871110 exm805378 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0701:0.0000:1.373:0.002:0.4607 +10 871746 exm805386 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1743:0.00239999999999996:0.597:0.017:0.5312 +10 875350 exm805392 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1814:0:0.924:0.018:0.5191 +10 888899 exm805400 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0585:0.00470000000000004:1.182:0.028:0.5116 +10 888916 exm805401 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1408:0.0000:1.480:0.041:0.4128 +10 894838 rs12249828 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.2119:0.4839:1.145:1.256:0.8432 +10 903950 rs11253471 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:-0.0047:0.4747:0.562:0.549:0.9470 +10 909757 exm805413 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0614:0:1.538:0.027:0.4097 +10 909766 exm805415 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0578:0:1.502:0.023:0.4128 +10 910074 exm805420 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0503:0.0071:1.001:0.057:0.5002 +10 910081 exm805421 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1352:0.002:1.125:0.033:0.4254 +10 927331 rs11253489 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0370:0.5325:0.976:1.435:0.7729 +10 931618 exm805430 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0474:0:1.222:0.009:0.4871 +10 931631 exm805431 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0819:0.00390000000000001:0.827:0.023:0.5296 +11 135656 rs10904561 T G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.0415:0.5176:0.543:0.661:0.9477 +11 135708 rs7917054 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.2441:0.501:0.742:0.659:0.9604 +11 178434 rs7089889 T G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1440:0.4833:0.836:0.827:0.8884 +11 188805 rs12146291 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.2144:0.4715:0.776:0.743:0.9610 +11 203471 rs10903451 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:0.4201:0.0000:0.619:0.022:0.9337 +11 252693 rs2379078 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.2673:0.0068:0.806:0.081:0.9018 +11 273344 rs2448378 A C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.7407:0.0000:0.497:0.012:0.9411 +11 279248 rs2496278 C A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.4228:0.6268:0.285:0.958:0.6772 +11 291134 rs10508201 T C . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1085:0.4723:0.317:0.444:0.9106 +11 293358 exm2249181 A C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0097:0.0044:1.507:0.061:0.3992 +11 313504 rs4880568 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.5836:0.9803:0.017:0.321:0.9547 +11 317501 rs2018975 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.3392:0:1.069:0.039:0.7852 +11 323283 rs3740304 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1836:0.5588:1.053:1.110:0.8640 +11 326894 rs3125027 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0517:0.0005:1.300:0.026:0.8916 +11 327162 exm805022 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1197:0:0.995:0.015:0.5328 +11 329493 rs3815985 C A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0480:0.5187:0.824:0.794:0.9109 +11 338553 rs3125031 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.3387:0.4696:0.777:1.004:0.7901 +11 354301 rs35198327 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.1338:0.4929:0.903:0.996:0.8546 +11 375475 variant.11341 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0747:0.000199999999999978:1.132:0.012:0.4773 +11 377161 rs7070654 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1388:0:0.582:0.020:0.9458 +11 387060 rs7904155 G A . . . GT:LRR:BAF:IA:IB:GC 1/1:0.0615:0.9945:0.080:1.883:0.8310 +11 390962 exm805090 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0087:9.9999999999989e-05:1.226:0.063:0.6950 +11 390983 exm805091 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.3606:0:0.398:0.035:0.4585 +11 394615 exm805105 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0170:0:1.345:0.020:0.4108 +11 402390 exm805131 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0341:0.0000:1.127:0.023:0.8735 +11 403792 exm805136 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1428:0.0000:0.908:0.019:0.8483 +11 410501 exm805166 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0189:0:1.162:0.050:0.4628 +11 413010 rs9787422 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0871:0.5334:0.908:0.961:0.8338 +11 423555 rs12245012 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1264:0.00149999999999995:1.658:0.040:0.8015 +11 423671 rs10904067 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.3622:0.5566:0.326:0.370:0.9306 +11 431161 exm2271240 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.0175:0.4838:0.935:1.004:0.8308 +11 435971 exm2249233 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.2716:0.0000:1.217:0.042:0.3992 +11 435990 exm805214 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1427:0:1.107:0.031:0.4295 +11 453567 rs4881254 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.1831:0.5079:1.032:1.425:0.8149 +11 462885 rs10904173 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.0930:0.4819:0.568:0.489:0.9332 +11 465046 exm805254 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1108:0:1.041:0.038:0.8704 +11 468599 rs12415961 T C . . . GT:LRR:BAF:IA:IB:GC 0/1:0.1647:0.5021:1.192:0.825:0.8051 +11 486935 exm805280 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1369:0.0186:0.643:0.060:0.4804 +11 487973 rs4881313 T C . . . GT:LRR:BAF:IA:IB:GC 0/1:0.0644:0.5628:0.747:0.992:0.8463 +11 521431 exm2271370 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:0.0436:0.9952:0.028:1.491:0.8545 +11 521723 rs4881336 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.2056:1:0.055:1.228:0.8344 +11 547567 rs816599 A G . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.0079:1.0000:0.023:1.101:0.9490 +11 554186 rs12251997 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1537:0.00309999999999999:1.565:0.032:0.8228 +11 556129 rs10466270 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:0.0014:0.5142:0.709:0.732:0.9257 +11 566379 rs11252926 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0211:0.5024:0.559:0.493:0.9438 +11 588406 rs10904450 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0834:0.0006:1.837:0.066:0.8107 +11 592283 rs816646 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.3907:0.5006:0.268:0.158:0.8540 +11 601089 rs816650 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0897:0:1.053:0.009:0.9557 +11 615782 rs17221309 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.2800:0.00319999999999998:1.764:0.035:0.7820 +11 635857 rs7069611 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0450:0.0013:0.593:0.017:0.9595 +11 669358 rs7898821 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1257:0:1.694:0.067:0.7270 +11 673037 rs816563 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:0.1368:0.9065:0.167:0.495:0.4895 +11 680428 rs2124585 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:-0.0590:0.5212:0.660:1.123:0.7133 +11 712617 rs17136372 A C . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.0174:0.9992:0.028:0.975:0.9556 +11 738630 rs1750792 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0269:0:1.575:0.025:0.8126 +11 749238 rs553326 G T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1758:0.00229999999999997:0.955:0.019:0.9300 +11 754033 rs2265090 C A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1347:0:1.181:0.025:0.9176 +11 757222 rs10904546 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:0.3035:0.0000:1.975:0.056:0.8047 +11 766105 rs7906313 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1672:0.0116:0.960:0.062:0.9381 +11 776027 rs11253377 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.1474:0.9964:0.046:1.017:0.9087 +11 787896 rs1769242 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.2049:0.0000:0.965:0.029:0.9107 +11 800213 rs2254501 G A . . . GT:LRR:BAF:IA:IB:GC 1/1:0.2422:0.9918:0.027:0.612:0.9374 +11 802445 rs1769215 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0188:0.00370000000000004:0.859:0.024:0.9303 +11 803721 rs2790381 C T . . . GT:LRR:BAF:IA:IB:GC 1/1:-0.0421:0.9798:0.029:0.485:0.9700 +11 811876 rs12414585 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0717:0.516:0.374:0.276:0.9351 +11 813426 rs11253424 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.5256:0:0.333:0.000:0.8993 +11 815149 rs7084027 A G . . . GT:LRR:BAF:IA:IB:GC 1/1:0.3183:0.9944:0.037:1.547:0.8419 +11 820868 rs10752019 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0749:0.0118:1.349:0.070:0.8882 +11 833528 exm2271241 T G . . . GT:LRR:BAF:IA:IB:GC 0/1:-0.0329:0.4896:0.569:0.681:0.9348 +11 834736 rs4881530 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0453:0.00249999999999995:1.578:0.040:0.8076 +11 836115 rs11253444 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:-0.1380:0.5199:0.596:0.670:0.9248 +11 838179 rs1536337 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0534:0.0000:1.318:0.095:0.7967 +11 858022 rs9124 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.5565:0:0.715:0.012:0.9397 +11 858924 exm805329 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0866:0:0.873:0.021:0.9056 +11 859047 exm805338 C A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1868:0:0.917:0.007:0.5296 +11 860687 exm805350 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1509:0:1.172:0.059:0.3777 +11 860726 exm805354 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0224:0.0000:2.255:0.138:0.3747 +11 860970 exm805360 A G . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1069:0.0000:1.859:0.070:0.7352 +11 871110 exm805378 G C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0701:0.0000:1.373:0.002:0.4607 +11 871746 exm805386 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1743:0.00239999999999996:0.597:0.017:0.5312 +11 875350 exm805392 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1814:0:0.924:0.018:0.5191 +11 888899 exm805400 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0585:0.00470000000000004:1.182:0.028:0.5116 +11 888916 exm805401 T C . . . GT:LRR:BAF:IA:IB:GC 0/0:0.1408:0.0000:1.480:0.041:0.4128 +11 894838 rs12249828 C T . . . GT:LRR:BAF:IA:IB:GC 1/0:0.2119:0.4839:1.145:1.256:0.8432 +11 903950 rs11253471 A G . . . GT:LRR:BAF:IA:IB:GC 0/1:-0.0047:0.4747:0.562:0.549:0.9470 +11 909757 exm805413 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0614:0:1.538:0.027:0.4097 +11 909766 exm805415 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0578:0:1.502:0.023:0.4128 +11 910074 exm805420 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0503:0.0071:1.001:0.057:0.5002 +11 910081 exm805421 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.1352:0.002:1.125:0.033:0.4254 +11 927331 rs11253489 G A . . . GT:LRR:BAF:IA:IB:GC 1/0:0.0370:0.5325:0.976:1.435:0.7729 +11 931618 exm805430 G A . . . GT:LRR:BAF:IA:IB:GC 0/0:0.0474:0:1.222:0.009:0.4871 +11 931631 exm805431 C T . . . GT:LRR:BAF:IA:IB:GC 0/0:-0.0819:0.00390000000000001:0.827:0.023:0.5296
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/cnv_baf_only.vcf Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,188 @@ +##fileformat=VCFv4.2 +##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> +##FORMAT=<ID=BAF,Number=1,Type=Float,Description="B Allele Frequency"> +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT test +10 135656 rs10904561 T G . . . GT:BAF 0/1:0.5176 +10 135708 rs7917054 G A . . . GT:BAF 1/0:0.501 +10 178434 rs7089889 T G . . . GT:BAF 0/1:0.4833 +10 188805 rs12146291 A G . . . GT:BAF 0/1:0.4715 +10 203471 rs10903451 A G . . . GT:BAF 0/0:0.0000 +10 252693 rs2379078 T C . . . GT:BAF 0/0:0.0068 +10 273344 rs2448378 A C . . . GT:BAF 0/0:0.0000 +10 279248 rs2496278 C A . . . GT:BAF 1/0:0.6268 +10 291134 rs10508201 T C . . . GT:BAF 0/1:0.4723 +10 293358 exm2249181 A C . . . GT:BAF 0/0:0.0044 +10 313504 rs4880568 C T . . . GT:BAF 1/1:0.9803 +10 317501 rs2018975 C T . . . GT:BAF 0/0:0 +10 323283 rs3740304 A G . . . GT:BAF 0/1:0.5588 +10 326894 rs3125027 T C . . . GT:BAF 0/0:0.0005 +10 327162 exm805022 C T . . . GT:BAF 0/0:0 +10 329493 rs3815985 C A . . . GT:BAF 1/0:0.5187 +10 338553 rs3125031 C T . . . GT:BAF 1/0:0.4696 +10 354301 rs35198327 C T . . . GT:BAF 1/0:0.4929 +10 375475 variant.11341 C T . . . GT:BAF 0/0:0.000199999999999978 +10 377161 rs7070654 C T . . . GT:BAF 0/0:0 +10 387060 rs7904155 G A . . . GT:BAF 1/1:0.9945 +10 390962 exm805090 G A . . . GT:BAF 0/0:9.9999999999989e-05 +10 390983 exm805091 G A . . . GT:BAF 0/0:0 +10 394615 exm805105 G A . . . GT:BAF 0/0:0 +10 402390 exm805131 G C . . . GT:BAF 0/0:0.0000 +10 403792 exm805136 G C . . . GT:BAF 0/0:0.0000 +10 410501 exm805166 C T . . . GT:BAF 0/0:0 +10 413010 rs9787422 C T . . . GT:BAF 1/0:0.5334 +10 423555 rs12245012 G A . . . GT:BAF 0/0:0.00149999999999995 +10 423671 rs10904067 C T . . . GT:BAF 1/0:0.5566 +10 431161 exm2271240 C T . . . GT:BAF 1/0:0.4838 +10 435971 exm2249233 T C . . . GT:BAF 0/0:0.0000 +10 435990 exm805214 G A . . . GT:BAF 0/0:0 +10 453567 rs4881254 C T . . . GT:BAF 1/0:0.5079 +10 462885 rs10904173 G A . . . GT:BAF 1/0:0.4819 +10 465046 exm805254 G A . . . GT:BAF 0/0:0 +10 468599 rs12415961 T C . . . GT:BAF 0/1:0.5021 +10 486935 exm805280 G C . . . GT:BAF 0/0:0.0186 +10 487973 rs4881313 T C . . . GT:BAF 0/1:0.5628 +10 521431 exm2271370 C T . . . GT:BAF 1/1:0.9952 +10 521723 rs4881336 C T . . . GT:BAF 1/1:1 +10 547567 rs816599 A G . . . GT:BAF 1/1:1.0000 +10 554186 rs12251997 C T . . . GT:BAF 0/0:0.00309999999999999 +10 556129 rs10466270 A G . . . GT:BAF 0/1:0.5142 +10 566379 rs11252926 C T . . . GT:BAF 1/0:0.5024 +10 588406 rs10904450 A G . . . GT:BAF 0/0:0.0006 +10 592283 rs816646 G A . . . GT:BAF 1/0:0.5006 +10 601089 rs816650 C T . . . GT:BAF 0/0:0 +10 615782 rs17221309 G A . . . GT:BAF 0/0:0.00319999999999998 +10 635857 rs7069611 T C . . . GT:BAF 0/0:0.0013 +10 669358 rs7898821 G A . . . GT:BAF 0/0:0 +10 673037 rs816563 C T . . . GT:BAF 1/1:0.9065 +10 680428 rs2124585 G A . . . GT:BAF 1/0:0.5212 +10 712617 rs17136372 A C . . . GT:BAF 1/1:0.9992 +10 738630 rs1750792 C T . . . GT:BAF 0/0:0 +10 749238 rs553326 G T . . . GT:BAF 0/0:0.00229999999999997 +10 754033 rs2265090 C A . . . GT:BAF 0/0:0 +10 757222 rs10904546 A G . . . GT:BAF 0/0:0.0000 +10 766105 rs7906313 T C . . . GT:BAF 0/0:0.0116 +10 776027 rs11253377 C T . . . GT:BAF 1/1:0.9964 +10 787896 rs1769242 T C . . . GT:BAF 0/0:0.0000 +10 800213 rs2254501 G A . . . GT:BAF 1/1:0.9918 +10 802445 rs1769215 C T . . . GT:BAF 0/0:0.00370000000000004 +10 803721 rs2790381 C T . . . GT:BAF 1/1:0.9798 +10 811876 rs12414585 G A . . . GT:BAF 1/0:0.516 +10 813426 rs11253424 C T . . . GT:BAF 0/0:0 +10 815149 rs7084027 A G . . . GT:BAF 1/1:0.9944 +10 820868 rs10752019 C T . . . GT:BAF 0/0:0.0118 +10 833528 exm2271241 T G . . . GT:BAF 0/1:0.4896 +10 834736 rs4881530 G A . . . GT:BAF 0/0:0.00249999999999995 +10 836115 rs11253444 A G . . . GT:BAF 0/1:0.5199 +10 838179 rs1536337 T C . . . GT:BAF 0/0:0.0000 +10 858022 rs9124 C T . . . GT:BAF 0/0:0 +10 858924 exm805329 G A . . . GT:BAF 0/0:0 +10 859047 exm805338 C A . . . GT:BAF 0/0:0 +10 860687 exm805350 C T . . . GT:BAF 0/0:0 +10 860726 exm805354 T C . . . GT:BAF 0/0:0.0000 +10 860970 exm805360 A G . . . GT:BAF 0/0:0.0000 +10 871110 exm805378 G C . . . GT:BAF 0/0:0.0000 +10 871746 exm805386 G A . . . GT:BAF 0/0:0.00239999999999996 +10 875350 exm805392 G A . . . GT:BAF 0/0:0 +10 888899 exm805400 C T . . . GT:BAF 0/0:0.00470000000000004 +10 888916 exm805401 T C . . . GT:BAF 0/0:0.0000 +10 894838 rs12249828 C T . . . GT:BAF 1/0:0.4839 +10 903950 rs11253471 A G . . . GT:BAF 0/1:0.4747 +10 909757 exm805413 G A . . . GT:BAF 0/0:0 +10 909766 exm805415 G A . . . GT:BAF 0/0:0 +10 910074 exm805420 C T . . . GT:BAF 0/0:0.0071 +10 910081 exm805421 C T . . . GT:BAF 0/0:0.002 +10 927331 rs11253489 G A . . . GT:BAF 1/0:0.5325 +10 931618 exm805430 G A . . . GT:BAF 0/0:0 +10 931631 exm805431 C T . . . GT:BAF 0/0:0.00390000000000001 +11 135656 rs10904561 T G . . . GT:BAF 0/1:0.5176 +11 135708 rs7917054 G A . . . GT:BAF 1/0:0.501 +11 178434 rs7089889 T G . . . GT:BAF 0/1:0.4833 +11 188805 rs12146291 A G . . . GT:BAF 0/1:0.4715 +11 203471 rs10903451 A G . . . GT:BAF 0/0:0.0000 +11 252693 rs2379078 T C . . . GT:BAF 0/0:0.0068 +11 273344 rs2448378 A C . . . GT:BAF 0/0:0.0000 +11 279248 rs2496278 C A . . . GT:BAF 1/0:0.6268 +11 291134 rs10508201 T C . . . GT:BAF 0/1:0.4723 +11 293358 exm2249181 A C . . . GT:BAF 0/0:0.0044 +11 313504 rs4880568 C T . . . GT:BAF 1/1:0.9803 +11 317501 rs2018975 C T . . . GT:BAF 0/0:0 +11 323283 rs3740304 A G . . . GT:BAF 0/1:0.5588 +11 326894 rs3125027 T C . . . GT:BAF 0/0:0.0005 +11 327162 exm805022 C T . . . GT:BAF 0/0:0 +11 329493 rs3815985 C A . . . GT:BAF 1/0:0.5187 +11 338553 rs3125031 C T . . . GT:BAF 1/0:0.4696 +11 354301 rs35198327 C T . . . GT:BAF 1/0:0.4929 +11 375475 variant.11341 C T . . . GT:BAF 0/0:0.000199999999999978 +11 377161 rs7070654 C T . . . GT:BAF 0/0:0 +11 387060 rs7904155 G A . . . GT:BAF 1/1:0.9945 +11 390962 exm805090 G A . . . GT:BAF 0/0:9.9999999999989e-05 +11 390983 exm805091 G A . . . GT:BAF 0/0:0 +11 394615 exm805105 G A . . . GT:BAF 0/0:0 +11 402390 exm805131 G C . . . GT:BAF 0/0:0.0000 +11 403792 exm805136 G C . . . GT:BAF 0/0:0.0000 +11 410501 exm805166 C T . . . GT:BAF 0/0:0 +11 413010 rs9787422 C T . . . GT:BAF 1/0:0.5334 +11 423555 rs12245012 G A . . . GT:BAF 0/0:0.00149999999999995 +11 423671 rs10904067 C T . . . GT:BAF 1/0:0.5566 +11 431161 exm2271240 C T . . . GT:BAF 1/0:0.4838 +11 435971 exm2249233 T C . . . GT:BAF 0/0:0.0000 +11 435990 exm805214 G A . . . GT:BAF 0/0:0 +11 453567 rs4881254 C T . . . GT:BAF 1/0:0.5079 +11 462885 rs10904173 G A . . . GT:BAF 1/0:0.4819 +11 465046 exm805254 G A . . . GT:BAF 0/0:0 +11 468599 rs12415961 T C . . . GT:BAF 0/1:0.5021 +11 486935 exm805280 G C . . . GT:BAF 0/0:0.0186 +11 487973 rs4881313 T C . . . GT:BAF 0/1:0.5628 +11 521431 exm2271370 C T . . . GT:BAF 1/1:0.9952 +11 521723 rs4881336 C T . . . GT:BAF 1/1:1 +11 547567 rs816599 A G . . . GT:BAF 1/1:1.0000 +11 554186 rs12251997 C T . . . GT:BAF 0/0:0.00309999999999999 +11 556129 rs10466270 A G . . . GT:BAF 0/1:0.5142 +11 566379 rs11252926 C T . . . GT:BAF 1/0:0.5024 +11 588406 rs10904450 A G . . . GT:BAF 0/0:0.0006 +11 592283 rs816646 G A . . . GT:BAF 1/0:0.5006 +11 601089 rs816650 C T . . . GT:BAF 0/0:0 +11 615782 rs17221309 G A . . . GT:BAF 0/0:0.00319999999999998 +11 635857 rs7069611 T C . . . GT:BAF 0/0:0.0013 +11 669358 rs7898821 G A . . . GT:BAF 0/0:0 +11 673037 rs816563 C T . . . GT:BAF 1/1:0.9065 +11 680428 rs2124585 G A . . . GT:BAF 1/0:0.5212 +11 712617 rs17136372 A C . . . GT:BAF 1/1:0.9992 +11 738630 rs1750792 C T . . . GT:BAF 0/0:0 +11 749238 rs553326 G T . . . GT:BAF 0/0:0.00229999999999997 +11 754033 rs2265090 C A . . . GT:BAF 0/0:0 +11 757222 rs10904546 A G . . . GT:BAF 0/0:0.0000 +11 766105 rs7906313 T C . . . GT:BAF 0/0:0.0116 +11 776027 rs11253377 C T . . . GT:BAF 1/1:0.9964 +11 787896 rs1769242 T C . . . GT:BAF 0/0:0.0000 +11 800213 rs2254501 G A . . . GT:BAF 1/1:0.9918 +11 802445 rs1769215 C T . . . GT:BAF 0/0:0.00370000000000004 +11 803721 rs2790381 C T . . . GT:BAF 1/1:0.9798 +11 811876 rs12414585 G A . . . GT:BAF 1/0:0.516 +11 813426 rs11253424 C T . . . GT:BAF 0/0:0 +11 815149 rs7084027 A G . . . GT:BAF 1/1:0.9944 +11 820868 rs10752019 C T . . . GT:BAF 0/0:0.0118 +11 833528 exm2271241 T G . . . GT:BAF 0/1:0.4896 +11 834736 rs4881530 G A . . . GT:BAF 0/0:0.00249999999999995 +11 836115 rs11253444 A G . . . GT:BAF 0/1:0.5199 +11 838179 rs1536337 T C . . . GT:BAF 0/0:0.0000 +11 858022 rs9124 C T . . . GT:BAF 0/0:0 +11 858924 exm805329 G A . . . GT:BAF 0/0:0 +11 859047 exm805338 C A . . . GT:BAF 0/0:0 +11 860687 exm805350 C T . . . GT:BAF 0/0:0 +11 860726 exm805354 T C . . . GT:BAF 0/0:0.0000 +11 860970 exm805360 A G . . . GT:BAF 0/0:0.0000 +11 871110 exm805378 G C . . . GT:BAF 0/0:0.0000 +11 871746 exm805386 G A . . . GT:BAF 0/0:0.00239999999999996 +11 875350 exm805392 G A . . . GT:BAF 0/0:0 +11 888899 exm805400 C T . . . GT:BAF 0/0:0.00470000000000004 +11 888916 exm805401 T C . . . GT:BAF 0/0:0.0000 +11 894838 rs12249828 C T . . . GT:BAF 1/0:0.4839 +11 903950 rs11253471 A G . . . GT:BAF 0/1:0.4747 +11 909757 exm805413 G A . . . GT:BAF 0/0:0 +11 909766 exm805415 G A . . . GT:BAF 0/0:0 +11 910074 exm805420 C T . . . GT:BAF 0/0:0.0071 +11 910081 exm805421 C T . . . GT:BAF 0/0:0.002 +11 927331 rs11253489 G A . . . GT:BAF 1/0:0.5325 +11 931618 exm805430 G A . . . GT:BAF 0/0:0 +11 931631 exm805431 C T . . . GT:BAF 0/0:0.00390000000000001
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/cnv_pairwise_summary.tab Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,6 @@ +# This file was produced by: bcftools cnv* +# The command line was: bcftools cnv .+ test .+ test * +# +# RG, Regions * +RG 10 135656 931631 2 2 21.7 92 27 92 27 +RG 11 135656 931631 2 2 21.7 92 27 92 27
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/cnv_summary.tab Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,7 @@ +# RG, Regions * +# This file was produced by: bcftools cnv* +# The command line was: bcftools cnv * +# +# RG, Regions * +RG 10 135656 931631 2 22.6 92 27 +RG 11 135656 931631 2 22.6 92 27
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/color_chrs.dat Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,5 @@ +# SG, shared segment [2]Chromosome [3]Start [4]End [5]B:1 [6]B:2 +# SW, number of switches [3]Sample [4]Chromosome [5]nHets [5]nSwitches [6]switch rate +SG 1 1 999 A:1 C:1 +SW A 1 0 0 0.000000 +SW C 1 0 0 0.000000
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/color_chrs_unrelated.dat Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,5 @@ +# SG, shared segment [2]Chromosome [3]Start [4]End [5]B:1 [6]B:2 +# SW, number of switches [3]Sample [4]Chromosome [5]nHets [5]nSwitches [6]switch rate +SG 1 1 999 A:1 - +SW A 1 0 0 0.000000 +SW (null) 1 0 0 0.000000
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/fasta_indexes.loc Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,6 @@ +23andme ? test genome ${__HERE__}/test-cache/23andme.fa +consensus ? test genome ${__HERE__}/test-cache/consensus.fa +csq ? test genome ${__HERE__}/test-cache/csq.fa +gvcf ? test genome ${__HERE__}/test-cache/gvcf.fa +mpileup ? test genome ${__HERE__}/test-cache/mpileup.ref.fa +norm ? test genome ${__HERE__}/test-cache/norm.fa
--- a/test-data/mpileup.ref.fa.fai Thu Feb 21 15:55:03 2019 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,1 +0,0 @@ -17 4200 14 60 61
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/plugin_frameshifts.bed Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,1 @@ +1 3106150 3106160
--- a/test-data/stats.a.vcf Thu Feb 21 15:55:03 2019 -0500 +++ b/test-data/stats.a.vcf Wed Jun 05 12:59:21 2019 -0400 @@ -3,6 +3,6 @@ ##contig=<ID=1,assembly=b37,length=249250621> ##reference=file:///lustre/scratch105/projects/g1k/ref/main_project/human_g1k_v37.fasta #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A B C -1 1000 . G A . PASS . GT 0/0 0/1 1/1 -1 1001 . G A . PASS . GT 0/0 0/1 1/1 -1 1002 . G A . PASS . GT 0/0 0/1 1/1 +1 1000 . G A . PASS . GT 0|0 0|1 1|1 +1 1001 . G A . PASS . GT 0/0 0|1 1|1 +1 1002 . G A . PASS . GT 0|0 0|1 1|1
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/23andme.fa Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,14 @@ +>1 +CACGTNACGGCTGAAGTCCAAGGTAC +CGTATCGAGTTCACAGTCGATAGCTC +GATCGATAGCATCGCTAGCNNNACTA +CGATCGATCGCTCTCCGTAACACTCA +AAAACGATCGATCGACTGCTCTTTAG +CGATGACTTTAGGGGAAAAA +>2 +CGCTCAGCCGTACAGCCGAGCAGGAC +ACGCTATTTTAGATCGACTGGCTNNG +CGCTAGCTACGCTTTAGCACGAGAA +>Y +NNNGCATACGTGTCCATCACGATGAT +AGCGATGATCGATC
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/23andme.fa.fai Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,3 @@ +1 150 3 26 27 +2 77 162 26 27 +Y 40 245 26 27
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/consensus.fa Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,20 @@ +>1:2-501 +TACCATATGTGACATATAAAAAAGAACATAACCTACGTATCAACTAAAGTGGTTGTTTGC +AGAAAAGGAAGACTTAAAAAGAGTCAGTACTAACCTACATAATATATACAATGTTCATTA +AATAATAAAATGAGCTCATCATACTTAGGTCATCATAAATATATCTGAAATTCACAAATA +TTGATCAAATGGTAAAATAGACAAGTAGATTTTAATAGGTTAAACAATTACTGATTCTCT +TGAAAGAATAAATTTAATATGAGACCTATTTCATTATAATGAACTCACAAATTAGAAACT +TCACACTGGGGGCTGGAGAGATGGCTCAGTAGTTAAGAACACTGACTGCTCTTCTGAAGG +TCCTGAGTTCAAATCCCAGCAACCACATGGTGACTTACAACCATCTGTAATGACATCTGA +TGCCCTCTGGTGTGTCTGAAGACAGCTACAGTGTACTTACATAAAATAATAAATAAATCT +TTAAAAACAAAAAAAAAGAA +>2 +GAAGATCTTTTCCTTATTAAGGATCTGAAGCTCTGTAGATTTGTATTCTATTAAACATGG +AGAGATTAGTGATTTTCCATATTCTTTAAGTCATTTTAGAGTAATGTGTTCTTAAGATAA +ATCAGAAAAACAAAAACTTGTGCTTTCCTGTTTGAAAAACAAACAGCTGTGGGGAATGGT +GTCGGGACAGCCTTTTTATAAAATTTTTCTAAATAATGTTGAGGCTTTGATACGTCAAAG +TTATATTTCAAATGGAATCACTTAGACCTCGTTTCTGAGTGTCAATGGCCATATTGGGGA +TTTGCTGCTGCCAATGACAGCACACCCTGGGAATGCCCCAACTACTTACTACAAAGCAGT +GTTACATGGAGAAGATCTTCAAGAGTCTTTTTGCTAGATCTTTCCTTGGCTTTTGATGTG +ACTCCTCTCAATAAAATCCACAGTAATATAGTGAGTGGTCTCCTGCTCCAAACCAGTATT +TCAGACACAGTTAATCCAGAC
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/consensus.fa.fai Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,2 @@ +1:2-501 500 9 60 61 +2 501 521 60 61
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/csq.fa Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,34 @@ +>1 +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +>2 +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +>3 +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +ACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGTACGT +
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/csq.fa.fai Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,3 @@ +1 600 3 60 61 +2 600 616 60 61 +3 600 1229 60 61
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/gvcf.fa Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,9 @@ +>22 +NNNNNNNNCCTTGGCCAAGTCACTTCCTCCTTCAGGAACATTGCAGTGGGCCTAAGTGCC +TCCTCTCGGGACTGGTATGGGGACGGTCATGCAATCTGGACAACATTCACCTTTAAAAGT +TTATTGATCTTTTGTGACATGCACGTGGGTTCCCAGTAGCAAGAAACTAAAGGGTCGCAG +GCCGGTTTCTGCTAATTTCTTTAATTCCAAGACAGTCTCAAATATTTTCTTATTAACTTC +CTGGAGGGAGGCTTATCATTCTCTCTTTTGGATGATTCTAAGTACCAGCTAAAATACAGC +TATCATTCATTTTCCTTGATTTGGGAGCCTAATTTCTTTAATTTAGTATGCAAGAAAACC +AATTTGGAAATATCAACTGTTTTGGAAACCTTAGACCTAGGTCATCCTTAGTAAGATCTT +CCCATTTATATAAATACTTGCAAGTAGTAGTGCCATAATT
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/gvcf.fa.fai Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,1 @@ +22 460 4 60 61
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/mpileup.ref.fa Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,71 @@ +>17 17:1-4200 +AAGCTTCTCACCCTGTTCCTGCATAGATAATTGCATGACAATTGCCTTGTCCCTGCTGAA +TGTGCTCTGGGGTCTCTGGGGTCTCACCCACGACCAACTCCCTGGGCCTGGCACCAGGGA +GCTTAACAAACATCTGTCCAGCGAATACCTGCATCCCTAGAAGTGAAGCCACCGCCCAAA +GACACGCCCATGTCCAGCTTAACCTGCATCCCTAGAAGTGAAGGCACCGCCCAAAGACAC +GCCCATGTCCAGCTTATTCTGCCCAGTTCCTCTCCAGAAAGGCTGCATGGTTGACACACA +GTGCCTGCGACAAAGCTGAATGCTATCATTTAAAAACTCCTTGCTGGTTTGAGAGGCAGA +AAATGATATCTCATAGTTGCTTTACTTTGCATATTTTAAAATTGTGACTTTCATGGCATA +AATAATACTGGTTTATTACAGAAGCACTAGAAAATGCATGTGGACAAAAGTTGGGATTAG +GAGAGAGAAATGAAGACATATGTCCACACAAAAACCTGTTCATTGCAGCTTTCTACCATC +ACCAAAAATTGCAAACAACCACACGCCCTTCAACTGGGGAACTCATCAACAACAAACTTG +TGGTTTACCCACACAATGGAAGACCACTTAGCAACAAAAAGGACCAAACTCCTGGTACAT +GCAACTGACAGATGAATCTCAAACGCATTCCTCCGTGTGAAAGAAGCCGGACTCACAGGG +CAACACACTATCTGACTGTTTCATGGGAAAGTCTGGAAACGGCAACACCATTGAGACAGA +AAACAGGTGAGTGGTTGCCTGGGGCCAGGGAACTTTCTGGGGTCATATTCTCTGTGTTGA +TTCTGGTGGTGGAAACAAGACTGTCCCAGCCTGGGTGATACAGCGAGACCCCATCTCTAC +CAAAAAATTAAAAATTAGCTGGGCATGGTGGTGCATGCCTGTAGTCCCAGCTATTCACAG +TGCTGAGGTGGGAAGATGCTTGAGCCCAGGAGTTCAAGGCTGCAATGAGCTATGATTGCG +CCACTGCACTTTGGCCTGGACAACAGAGCAAAACCCTGTCTCTAAAAAAAGAAAAGAAAA +GAAAAACTCACTGGATATGAATGATACAGGTTGAGGATCCATTATCTGAAATGCTTGGAC +CAGATGTTTTGAATTTTGGATTTTTTCATATTTTGTAATCTTTGCAGTATATTTACCAGT +TCAGCATCCCTAACTCAAAAATTCAAAAATCTGAAATCCCAAACGCGCCAATAAGCATTC +CCTTTGAGCGTCATGTCGGTGCTTGGAATGTTTGGGGTTTTGGATTTACAGCTTTGGGAC +GCTCAACCTGTACCTCAATAAACCTGATTTTAAAAAAGTTTGGGGGGATTCCCCTAAGCC +CGCCACCCGGAGACAGCGGATTTCCTTAGTTACTTACTATGCTCCTTGGCCATTTCTCTA +GGTATTGGTATATTGTGTCTGCTGTGAACTGTCCTTGGCCTGTTTGGTGACGGGTGAGGA +GCAGGGACAGAAGGGTCCTGCGTGCCCTGCCTTCACAAGCCCCTGGAAGGAAAGTTGTTT +TGGGATCTCTGCACCCTCAGCCTGGACAACTTGTGCCCATCTGGTGACCCCTCACTCAGC +CACCAGACTTCCACGACAGGCTCCAGCCTCGGCACCTTCAGCCATGGACAGTTCCGCCAG +CGTTGCCCTCTGTTCTGCTGTTTTCTCTACCAGAAGTGCCCTTCCCTCCTCACCTGACCA +CTCTGGGGAAATCCCTCAGCACCCTCCCTGAGCATACCCTACTCTGGCACAAGCCCACCC +TGCAAAGCCCCTGAGGCCCGCCCTGTGGCGTCTCTCCCTCCCTTGCTGTCAGGACAGTGG +TCCTGGCCACCGGGGCTCACGGAGCCGCCCTGTGCCGTGTACCTCTGAGCCCTCTGCACA +GTGCCTTCTGCTTGCCTGTGGCTTTGAGAAGAAACCCCTTCTGGTTATACATAAGACAGC +CAGAGAAGGGAGTTGCCCAGGGTGGCACAGCACGTTGCTGCCAGTTACTGCCATTTTCAC +GGGCATGAAATGGAGATAACAACAGGAGCGACCGCACAGGCTGCTGAGCGCGTCACACGC +AGCCATCGCGCAGCTCAGGGATATTACGTGTAACTCGACATGTCAGCGATTGTCACAGGC +ACTGCTACTCCTGGGGTTTTCCATCAAACCCTCAAGAGCTGGGCCTGGGGTCAACTTCCG +GCCTGGGGAAACTGGGGCAAGTATCACCAGAGATGAGCTTTATAAAAATAATGGTGCTAG +CTGGGCATGGTGGCTTGCACCTGTAATCCCAGCACTTTGGGAGGCCGAGCTAGGAGGATC +GTTTGAGTCCAGCAGTTTGAGACCAGCCTGGCCAATACGGCAAAACCCAGTCTCTACAAA +AAATACAAAAAACAACTAGCCAGGCGTGGTGGTGCACACCTGTAGTCCCAGCTACTCAGG +AGGCTGAGGGGGAAGGACTGCTTGAGCCCAGGAGTTTGAGGCTGCTGTGAGCTGTGATCG +CATCACTGCATTCCAGCCCGGTGACAGAGTGAGTCACTGTCTCAAAAAAGAAAGGAAGAA +ATAAAGAAAACAAATAAAAATAATAGTGCAGACAAAAGGCCTTGACCCATCTAGCTTTGG +CCCTCAGCATCAACCGCTAGATACGTCCCTCCCTTTCTTCTGGGGCACAGGTCACACTCT +CTTCCAGGTCTAGGATGCAGCTGAGGGGTGCCCCTCTTACCATCTAATCTGTGCCCTTAT +TTCCTCTGCTTTAGTGAGGAAGAGGCCCCTGGTCCATGAAGGGGCCTTTCAGAGACGGGG +ACCCCTGAGGAGCCCCGAGCAGCAGCCGTCGTGTCTCACCCAGGGTGTCTGAAACAGATG +TGGAGGTCTCGGGTGAGGCGTGGCTCAGATACAGGGAGTGGCCCACAGCTCGGCCTGTCT +TTGAAAGGCCACGTGACCTGGCCCACGGCTGGCAGGTGGGACCCAGCTGCAGGGGTCCAG +CAGCACCCACAGCAGCCACCTGTGGCAGGGAGGAGCTTGTGGTACAGTGGACAGGCCCTG +CCCAGATGGCCCCCCGCCTGCCTGTGGAAGTTGACCAGACCATCTGTCACAGCAGGTAAG +ACTCTGCTTTCTGGGCAACCCAGCAGGTGACCCTGGAATTCCTGTCCATCTGGCAGGTGG +GCATTGAAACTGGTTTAAAAATGTCACACCATAGGCCGGGCACAGTGGCTCACGCCTGTA +ATCCCAGCCCTTTGGGAGGCCAGGGTGGGTGGATCACTTGAGGTCAGGAGTTCAAGACCA +GCCTGGCCAACATGGTGAAACCCCGTCTACTAAAAATACAAAAATTAGCCTGGCGTGGTG +GCGCATGCCTGTAATCCCAGCTACTTGGGAAGCTGAGGGATGAGAACTGCTTGAACCTGG +GAGGCAGACGTTGCAGTGAGCTGAGATCACGCCACTGCACTCCAGCCTGGGCAACAGAGT +AAGACTCTGTCTCAAAAAAAAAAAAATCACACCATTTTGGCTTCAGATTGCATATCCTCC +TGCAAGGATATATACGCGTGAAATTCAAGTCAATGACAAATCAGAAGAAAAAACATATAT +ATACGCAAACCAGTATCCTACTGTGTGTGTCGTTTGTTGTGTTTTCGACAGCTGTCCGTG +TTATAATAATTCCTCTAGTTCAAATTTATTCATTTTTAACTTCATAGTACCACATTCTAC +ACACTGCCCATGTCCCCTCAAGCTTCCCCTGGCTCCTGCAACCACAAATCTACTCTCTGC +CTCTGTGGGTTGACCTATTCTGGACACGTCATAGAAATAGAGTCCTGCAACACGTGGCCG +TCTGTGTCTGGCTTCTCTCGCTTAGCATCTTGTTTCCAAGGTCCTCCCACAGTGTAGCAT +GCACCTGCTACACTCCTTCTTAGGGCTGATATTCCACGCACCTGCTACACTCCTTCTTAT +GGCTGATATTCCACGCACCTGCTACACTCCTTCTTAGGGCTGATATTCCACACACCCGCT +ACACTCCTTCTTAGGGCTGATATTCCACGCACCCGCTACACTCCTTCTTAGGGCTGATAT +TCCACGCACCTGCTACACTCCTTCTTAGGGCTGATATTCCACGCACCTGCTACACTCCTT +CTTAGGGCTGATATTCCACGCACCTGCTACACTCCTTCTTAGGGCTGATATTCCACGCAC
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/mpileup.ref.fa.fai Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,1 @@ +17 4200 14 60 61
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/norm.fa Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,26 @@ +>20 20:1339000-1339300 +AGGATGGGGCTCATTAATAGAGCTCCACTTGTCTCCAGAATCACTGGTGAGGAAGGGGAG +TGTTGCCCCCACATTCGTGCACAGCAGGGATGGTTCACCGAACTCCACACCAGTCTCTGC +AGAGCCTGTTGGGGAGAGGAGGGCTGTGGTTTCTTTGATGGTGTTCACCTGGAGTAGAGC +AAGTATTGTCAAAAGGGTCATCCTCGGAGGTTGCAGTGAGCCGAGATCGCACCATTGCAC +TGCAGCCTGGGAGACAGAGCAAGACTCCATCTCAAAAAAAAAAAAAAAAAAAAAGGCCAT +C +>1 1:10143-10443 +CTAACCCCTAACCCTAACCCTAACCCTAACCCTAACCTAACCCTAACCCTAACCCTAACC +CTAACCCTAACCCTAACCCTAACCCTAACCCCTAACCCTAACCCTAAACCCTAAACCCTA +ACCCTAACCCTAACCCTAACCCTAACCCCAACCCCAACCCCAACCCCAACCCCAACCCCA +ACCCTAACCCCTAACCCTAACCCTAACCCTACCCTAACCCTAACCCTAACCCTAACCCTA +ACCCTAACCCCTAACCCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCC +>2 1:1382388-1382602 +GGGCGTCTCATAGCTGGAGCAATGGCGAGCGCCTGGACAAGGGAGGGGAAGGGGTTCTTA +TTACTGACGCGGGTAGCCCCTACTGCTGTGTGGTTCCCCTATTTTTTTTTTTTTCTTTTT +GAGACGGAGTCTCGCTCTGTCACCCAGGCTGGAGTGCAGTGGCACAATCTCGGCTCACTG +CAAGCTCCACCTCCTGGGTTCACGCCATTCTCCTG +>3 madeup +ACTGGACACGTGGACACACACACACACACACACACACACACAGTCAAACCACCTACCAGA +>4 20:8917026-8917085 +TCCCCTCTTGACCTCTCTCTATTTTTTTTTTTTTTTCTGAGATGGATTTTTGCTCTTGTT +>5 20:18724313-18724343 +GTCTCAAAAAAAAAAAAAAAAAAAAGAAAAG +>21 +TTTATTATTATTATTATTAAATTGAATTTATTTAGTGTACATACATTCATGTGTATTGTG
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test-cache/norm.fa.fai Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,7 @@ +20 301 23 60 61 +1 300 347 60 61 +2 215 673 60 61 +3 60 902 60 61 +4 60 985 60 61 +5 31 1070 31 32 +21 60 1106 60 61
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_data_table_conf.xml.test Wed Jun 05 12:59:21 2019 -0400 @@ -0,0 +1,7 @@ +<tables> + <!-- Locations of FASTA index ffiles for testing --> + <table name="fasta_indexes" comment_char="#"> + <columns>value, dbkey, name, path</columns> + <file path="${__HERE__}/test-data/fasta_indexes.loc" /> + </table> +</tables>