Mercurial > repos > yhoogstrate > crossmap
changeset 0:bfcdabc062ca draft
planemo upload for repository https://github.com/ErasmusMC-Bioinformatics/crossmap_galaxy_wrapper commit d4e9fe61901a612c78d9f26f172537b27fd2ddbb
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/README.rst Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,7 @@ +CrossMap wrapper for Galaxy +=========================== + +CrossMap is versatile tool to convert genome coordinates or annotation files between genome +assemblies. It supports mostly commonly used file types, including BAM, BED,BigWig, GFF, +GTF, SAM, Wiggle, and VCF formats. For large plain text file types, such as BED, GFF, GTF +and VCF, reading from remote servers and file compression are supported.
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/crossmap.xml Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,455 @@ +<tool id="crossmap" name="CrossMap" version="0.2.a"> + <description>Convert genome coordinates or annotation files between genome assemblies</description> + + <requirements> + <requirement type="package" version="312">ucsc_tools</requirement> + <requirement type="package" version="0.2">crossmap</requirement> + </requirements> + + <stdio> + <regex match="Usage: CrossMap.py" source="stdout" level="fatal"/> + <regex match=".*" source="both" level="log"/> + </stdio> + + <version_command>CrossMap.py 2>&1 | head -n 1 | grep -E --only-matching 'CrossMap.*'</version_command> + + <command> + CrossMap.py + ${multiple.input_format.replace("sam","bam")} + + #if $multiple.input_format == "vcf" and $multiple.seq_source.index_source == "cached" + <!-- This is the 2nd dbkey, and the corresponding value has to be looked up --> + "${filter(lambda x: str( x[1] ) == str($multiple.seq_source.input_chain ), $__app__.tool_data_tables['liftOver'].get_fields())[0][2] }" + #else + "$multiple.seq_source.input_chain" + #end if + + #if $multiple.input_format in ["bam", "sam"] + -m $multiple.insert_size + -s $multiple.insert_size_stdev + -t $multiple.insert_size_fold + #end if + + "$multiple.seq_source.input" + + #if $multiple.input_format == "vcf" + "$multiple.seq_source.input_fasta" + #end if + + #if str($multiple.include_fails) == "True" + > + #end if + + "$output" + + #if $multiple.input_format in ["bam", "sam"] + ; mv "${output}.${multiple.input_format}" "$output" + ; mv "${output}.unmap.${multiple.input_format}" "$output_unmapped" + #else if $multiple.input_format in ["vcf"] + ; mv "${output}" "$output" + ; mv "${output}.unmap" "$output_unmapped" + #else if $multiple.input_format in ["wig", "bigwig"] + ; mv "${output}.bw" "$output" + ; mv "${output}.sorted.bgr" "$output2" + #end if + </command> + + <inputs> + <conditional name="multiple"> + <param name="input_format" type="select" label="Convert a file of the following format"> + <option value="bam">BAM</option> + <option value="sam">SAM</option> + <option value="bed">BED or BED-like</option> + <option value="bigwig">BigWig</option> + <option value="gff">GFF or GTF</option> + <option value="vcf">VCF</option> + <option value="wig">Wiggle or bedGraph</option> + </param> + <when value="bam"> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data (chain file)"> + <option value="cached">Local data (in galaxy)</option> + <option value="history">From History</option> + </param> + <when value="cached"> + <param type="data" format="bam" name="input" label="BAM/SAM file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + <param name="input_chain" type="select" label="Lift Over To"> + <options from_file="liftOver.loc"> + <column name="name" index="1"/> + <column name="value" index="2"/> + <column name="dbkey" index="0"/> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + </when> + <when value="history"> + <param type="data" format="bam,sam" name="input" label="BAM/SAM file" /> + <param type="data" format="csv" name="input_chain" label="LiftOver chain file" /> + </when> + </conditional> + + <param name="insert_size" type="float" value="200.0" label="Insert size (-m)" help="Average insert size of pair-end sequencing (bp) [default=200.0]" /> + <param name="insert_size_stdev" type="float" value="30.0" label="Insert size std. dev (-s)" help="Stanadard deviation of insert size. [default=30.0]" /> + <param name="insert_size_fold" type="float" value="3.0" label="Insert size std. dev foldchange (-t)" help="A mapped pair is considered as 'proper pair' if both ends mapped to different strand and the distance between them is less then '-t' * stdev from the mean. [default=3.0]" /> + + <param name="include_fails" type="hidden" tvalue="False" /> + </when> + <when value="sam"><!-- BAM and SAM are exactly the same conditions, but they need to be separate to get the proper output format --> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data (chain file)"> + <option value="cached">Local data (in galaxy)</option> + <option value="history">From History</option> + </param> + <when value="cached"> + <param type="data" format="sam" name="input" label="BAM/SAM file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + <param name="input_chain" type="select" label="Lift Over To"> + <options from_file="liftOver.loc"> + <column name="name" index="1"/> + <column name="value" index="2"/> + <column name="dbkey" index="0"/> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + </when> + <when value="history"> + <param type="data" format="bam,sam" name="input" label="BAM/SAM file" /> + <param type="data" format="csv" name="input_chain" label="LiftOver chain file" /> + </when> + </conditional> + + <param name="insert_size" type="float" value="200.0" label="Insert size (-m)" help="Average insert size of pair-end sequencing (bp) [default=200.0]" /> + <param name="insert_size_stdev" type="float" value="30.0" label="Insert size std. dev (-s)" help="Stanadard deviation of insert size. [default=30.0]" /> + <param name="insert_size_fold" type="float" value="3.0" label="Insert size std. dev foldchange (-t)" help="A mapped pair is considered as 'proper pair' if both ends mapped to different strand and the distance between them is less then '-t' * stdev from the mean. [default=3.0]" /> + + <param name="include_fails" type="hidden" tvalue="False" /> + </when> + <when value="bed"> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data (chain file)"> + <option value="cached">Local data (in galaxy)</option> + <option value="history">From History</option> + </param> + <when value="cached"> + <param format="bed" name="input" type="data" label="BED file" help="BED format file must have at least 3 columns (chrom, start, end) and no more than 12 columns."> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + <param name="input_chain" type="select" label="Lift Over To"> + <options from_file="liftOver.loc"> + <column name="name" index="1"/> + <column name="value" index="2"/> + <column name="dbkey" index="0"/> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + </when> + <when value="history"> + <param type="data" format="bed" name="input" label="BED file" help="BED format file must have at least 3 columns (chrom, start, end) and no more than 12 columns." /> + <param type="data" format="csv" name="input_chain" label="LiftOver chain file" /> + </when> + </conditional> + + <param name="include_fails" type="boolean" truevalue="True" checked="false" falsevalue="False" label="Include failed liftovers" help="If a coordinate can not be lift over, do you want to include it in the output (it is still being marked 'fail')" /> + </when> + <when value="bigwig"> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data (chain file)"> + <option value="cached">Local data (in galaxy)</option> + <option value="history">From History</option> + </param> + <when value="cached"> + <param format="bigwig" name="input" type="data" label="BigWig file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + <param name="input_chain" type="select" label="Lift Over To"> + <options from_file="liftOver.loc"> + <column name="name" index="1"/> + <column name="value" index="2"/> + <column name="dbkey" index="0"/> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + </when> + <when value="history"> + <param type="data" format="bigwig" name="input" label="BigWig file" /> + <param type="data" format="csv" name="input_chain" label="LiftOver chain file" /> + </when> + </conditional> + + <param name="include_fails" type="hidden" tvalue="False" /> + </when> + <when value="gff"> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data (chain file)"> + <option value="cached">Local data (in galaxy)</option> + <option value="history">From History</option> + </param> + <when value="cached"> + <param format="gtf,gff,gff3" name="input" type="data" label="GTF/GFF file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + <param name="input_chain" type="select" label="Lift Over To"> + <options from_file="liftOver.loc"> + <column name="name" index="1"/> + <column name="value" index="2"/> + <column name="dbkey" index="0"/> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + </when> + <when value="history"> + <param type="data" format="gtf,gff,gff3" name="input" label="GTF/GFF file" /> + <param type="data" format="csv" name="input_chain" label="LiftOver chain file" /> + </when> + </conditional> + + <param name="include_fails" type="boolean" truevalue="True" checked="false" falsevalue="False" label="Include failed liftovers" help="If a coordinate can not be lift over, do you want to include it in the output (it is still being marked 'fail')" /> + </when> + <when value="vcf"> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data"> + <option value="cached">Local data (in galaxy)</option> + <option value="history_chain">Chain file from History</option> + <option value="history_all">Chain & FASTA files from History</option> + </param> + + <when value="cached"> + <param type="data" format="vcf" name="input" label="VCF file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + + <!-- automatically fetch a FASTA file from the same DBKEY as the chain file --> + <param name="input_chain" type="select" label="Lift Over To (Chain file)"> + <options from_file="liftOver.loc"> + <column name="name" index="1" /> + <column name="value" index="1" /><!-- It is not possible to send the *.chain file as value, and obtain the 2nd dbkey as parameter via a filter --> + <column name="dbkey" index="0" /> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + + <param name="input_fasta" type="select" label="Lift Over To (FASTA file)" help="The FASTA file must be on the same build (dbkey) as the LiftOver chain file"> + <options from_file="all_fasta.loc"> + <column name="name" index="2"/> + <column name="value" index="3"/> + <column name="dbkey" index="1"/> + <filter type="param_value" ref="input_chain" column="1" /> + </options> + </param> + </when> + + <when value="history_chain"> + <param type="data" format="vcf" name="input" label="VCF file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="all_fasta.loc" metadata_name="dbkey" metadata_column="1" message="LiftOver mapping (FASTA file) is not available for the specified build." /> + </param> + <param type="data" format="csv" name="input_chain" multiple="false" label="LiftOver chain file" /> + <param name="input_fasta" type="select" label="Lift Over To (FASTA file)" help="The FASTA file must be on the same build (dbkey) as the LiftOver chain file"> + <options from_file="all_fasta.loc"> + <column name="name" index="2"/> + <column name="value" index="3"/> + <column name="dbkey" index="1"/> + </options> + </param> + </when> + + <when value="history_all"> + <param type="data" format="vcf" name="input" label="VCF file" /> + <param type="data" format="csv" name="input_chain" multiple="false" label="LiftOver chain file" /> + <param type="data" format="fasta" name="input_fasta" multiple="false" label="Full genome FASTA file" /> + </when> + </conditional> + + <param name="include_fails" type="hidden" tvalue="False" /> + </when> + <when value="wig"> + <conditional name="seq_source"> + <param name="index_source" type="select" label="Source for LiftOver Data (chain file)"> + <option value="cached">Local data (in galaxy)</option> + <option value="history">From History</option> + </param> + <when value="cached"> + <param format="wig" name="input" type="data" label="Wiggle file"> + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_file" filename="liftOver.loc" metadata_name="dbkey" metadata_column="0" message="LiftOver mapping (chain file) is not available for the specified build." /> + </param> + <param name="input_chain" type="select" label="Lift Over To"> + <options from_file="liftOver.loc"> + <column name="name" index="1"/> + <column name="value" index="2"/> + <column name="dbkey" index="0"/> + <filter type="data_meta" ref="input" key="dbkey" column="0" /> + </options> + </param> + </when> + <when value="history"> + <param format="wig" name="input" type="data" label="Wiggle file" /> + <param type="data" format="csv" name="input_chain" label="LiftOver chain file" /> + </when> + </conditional> + + <param name="include_fails" type="hidden" tvalue="False" /> + </when> + </conditional> + </inputs> + + <outputs> + <data format="text" name="output" label="${tool.name} on "> + <change_format> + <when input="multiple.input_format" value="bam" format="bam" /> + <when input="multiple.input_format" value="sam" format="sam" /> + <when input="multiple.input_format" value="bed" format="bed" /> + <when input="multiple.input_format" value="bigwig" format="bigwig" /> + <when input="multiple.input_format" value="gff" format="gff" /> + <when input="multiple.input_format" value="vcf" format="vcf" /> + <when input="imultiple.nput_format" value="wig" format="bigwig" /> + </change_format> + </data> + + <data format="text" name="output_unmapped" label="${tool.name} unmapped on "> + <filter>input_format in ["bam" , "vcf"]</filter> + <change_format> + <when input="multiple.input_format" value="bam" format="bam" /> + <when input="multiple.input_format" value="vcf" format="vcf" /> + </change_format> + </data> + <data format="text" name="output2" label="${tool.name} on "> + <filter>input_format in ["wig"]</filter> + <change_format> + <when input="multiple.input_format" value="wig" format="bedgraph" /> + </change_format> + </data> + </outputs> + + <tests> + <!-- BAM/SAM --> + <test> + <param name="input_format" value="sam"/> + <param name="index_source" value="history"/> + <param name="input" value="test_bam_01_input_a.sam" ftype="sam"/> + <param name="input_chain" value="aToB.over.chain" ftype="csv"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_bam_01_output_a.sam" compare="diff" lines_diff="4"/> + <output name="output_unmapped" file="test_bam_01_output_a.unmap.sam"/> + </test> + + <!-- BED --> + <test> + <param name="input_format" value="bed"/> + <param name="index_source" value="history"/> + <param name="input" value="test_bed_01_input_a.bed" ftype="bed"/> + <param name="input_chain" value="aToB.over.chain" ftype="csv"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_bed_01_output_a__only-matches.bed"/> + </test> + <test> + <param name="input_format" value="bed"/> + <param name="index_source" value="history"/> + <param name="input" value="test_bed_01_input_a.bed" ftype="bed"/> + <param name="input_chain" value="aToB.over.chain" ftype="txt"/> + <param name="include_fails" value="True"/> + + <output name="output" file="test_bed_01_output_a__all.bed"/> + </test> + <test> + <param name="input_format" value="bed"/> + <param name="index_source" value="history"/> + <param name="input" value="test_bed_02_input_a.bed" ftype="bed"/> + <param name="input_chain" value="aToB.over.chain" ftype="txt"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_bed_02_output_a__only-matches.bed"/> + </test> + <test> + <param name="input_format" value="bed"/> + <param name="index_source" value="history"/> + <param name="input" value="test_bed_02_input_a.bed" ftype="bed"/> + <param name="input_chain" value="aToB.over.chain" ftype="txt"/> + <param name="include_fails" value="True"/> + + <output name="output" file="test_bed_02_output_a__all.bed"/> + </test> + + <!-- BigWig --><!-- Malfuncioning in CrossMap 0.2, but patched via galaxy toolshed installer --> + <test> + <param name="input_format" value="bigwig"/> + <param name="index_source" value="history"/> + <param name="input" value="test_bigwig_01_input_a.bw" ftype="bigwig"/> + <param name="input_chain" value="aToB.over.chain" ftype="csv"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_bigwig_01_output_a.bw"/> + <output name="output2" file="test_bigwig_01_output_a.sorted.bgr"/> + </test> + <!-- GFF --> + <test> + <param name="input_format" value="gff"/> + <param name="index_source" value="history"/> + <param name="input" value="test_gff_01_input_a.gtf" ftype="gtf"/> + <param name="input_chain" value="aToB.over.chain" ftype="csv"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_gff_01_output_a__only-matches.gtf"/> + </test> + <test> + <param name="input_format" value="gff"/> + <param name="index_source" value="history"/> + <param name="input" value="test_gff_01_input_a.gtf" ftype="gtf"/> + <param name="input_chain" value="aToB.over.chain" ftype="csv"/> + <param name="include_fails" value="True"/> + + <output name="output" file="test_gff_01_output_a__all.gtf"/> + </test> + + <!-- VCF --> + <test> + <param name="input_format" value="vcf"/> + <param name="index_source" value="history_all"/> + <param name="input" value="test_vcf_01_input.vcf" ftype="vcf"/> + <param name="input_chain" value="test_vcf_01.over.chain" ftype="csv"/> + <param name="input_fasta" value="test_vcf_01.fasta" ftype="fasta"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_vcf_01_output.vcf" compare="diff" lines_diff="4"/> + <output name="output_unmapped" file="test_vcf_01_output.vcf.unmap"/> + </test> + </tests> + + <!-- WIG - Doesn't understand fixedStep --> + <test> + <param name="input_format" value="wig"/> + <param name="index_source" value="history"/> + <param name="input" value="test_wig_01_input_a.wig" ftype="wig"/> + <param name="input_chain" value="aToB.over.chain" ftype="csv"/> + <param name="include_fails" value="False"/> + + <output name="output" file="test_wig_01_output_a.bw"/> + <output name="output2" file="test_wig_01_output_a.sorted.bgr"/> + </test> + <help> +CrossMap is versatile tool to convert genome coordinates or annotation files between genome +assemblies. It supports mostly commonly used file types, including BAM, BED,BigWig, GFF, +GTF, SAM, Wiggle, and VCF formats. For large plain text file types, such as BED, GFF, GTF +and VCF, reading from remote servers and file compression are supported. + +CrossMap bed +------------ +BED format file must have at least 3 columns (chrom, start, end) and no more than 12 columns. +BED format: http://genome.ucsc.edu/FAQ/FAQformat.html#format1 + </help> + + <citations> + <citation type="doi">10.1093/bioinformatics/btt730</citation> + </citations> +</tool>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/aToB.over.chain Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,498 @@ +chain 21270171362 chr1 249250621 + 10000 249233096 chr1 247249719 + 0 247199719 2 +619 137 0 +166661 50000 50000 +40302 50000 50000 +153649 50000 50000 +1098479 1 1 +47 1 1 +73 117 114 +773 1 1 +43 1 1 +864369 2 2 +51 3 0 +104 13694 13694 +104 0 3 +51 2 2 +134936 50000 50000 +1161048 150000 60000 +1440092 27273 50000 +7590365 50000 50000 +116914 100000 50000 +237250 50000 50000 +3518496 50000 50000 +12702424 150000 50000 +16145012 0 1 +7772 1 0 +4705841 0 1 +52977198 50000 50000 +157344 21065 50000 +16604841 50000 50000 +189539 150000 50000 +398739 21050000 20290000 +195588 50000 50000 +186739 150000 50000 +175055 50000 50000 +201709 100000 50000 +126477 130183 50000 +381 0 3 +315 0 2 +62 1 1 +45 1 0 +19 0 1 +8 1 1 +1158 1 0 +314 12 13 +2849 3 0 +5615 1 1 +37 1 1 +3172 6 4 +190 1 1 +34 1 1 +380 0 1 +2099 3 0 +765 2 2 +366 0 4 +1186 1 0 +460 0 1 +1242 28 28 +574 1 1 +21 1 1 +1460 1 1 +105 1 1 +701 3 0 +239 0 1 +970 11 11 +2365 1 1 +21 1 1 +384 8 8 +996 2 0 +10383 2 0 +713 0 1 +5188 1 1 +30 1 1 +1233 1 0 +132 1 1 +44 1 1 +1123 22 22 +1810 1 0 +512 1 1 +39 1 1 +1096 0 16 +743 1 0 +9028 0 1 +2506 0 6 +8446 2 0 +5505 0 3 +7541 1 0 +109864 50000 50000 +78698 50000 50000 +127263 50000 50000 +170669 50000 50000 +38311 100000 100000 +1022394 50000 50000 +281532 289973 50000 +1648 1539 0 +76 3225 3 +1018 34 34 +1716 2 0 +159 4 0 +1600 0 6 +1385 1 0 +2066 1 1 +32 1 3 +1556908 150000 50000 +185320 150000 50000 +172789 50000 50000 +220313 50000 50000 +455185 50000 50000 +22047237 1 0 +34365824 150000 50000 +259514 150000 50000 +17265625 50000 50000 +11394365 50000 50000 +13665999 150000 50000 +174886 + +chain 22466185064 chr2 243199373 + 10000 243102476 chr2 242951149 + 0 242751149 1 +1201236 0 1 +2057 1 1 +26 1 1 +4043 22 0 +82 7 101 +33 84 0 +26 1 23 +67 48 0 +176 318 0 +40 92 0 +33 0 74 +62 0 6 +22 0 118 +387 2446 0 +843 2450 0 +167 2304 0 +21342 1 1 +50 4 0 +2040 0 4 +9053 2 0 +88 3 53 +1078 0 1 +4122 1 0 +1608 28 1 +4029 0 9 +967 4547 1000 +420 27 0 +2237126 0 1 +2458 17 17 +4595 4 0 +3531 14 14 +3954 51129 50000 +1426129 160449 100025 +11087286 0 3 +773 103977 50000 +1576 0 1 +240 0 1 +318 2 0 +255 0 16 +2851 10 10 +1404 1 1 +27 1 1 +703 13 0 +375 0 1 +4783388 35000 25000 +848 0 2 +1225 0 2 +391 1 0 +3653 1 1 +28 1 1 +49 0 4 +4514 0 9 +9628 1 9 +1483 7 7 +1990 0 1 +1949 1 0 +2716230 1 0 +7799848 851 851 +38932068 5 9 +9477928 0 3 +7529698 72396 0 +1891551 294073 150000 +17672 1 0 +5792 1 0 +9017 0 266 +2944 0 1 +9388 2 0 +37828 4 0 +23436 0 4 +290931 1273578 1000000 +731068 3000000 3000000 +2558421 0 1 +1197 10 0 +685 42 46 +2309 0 1 +1332 22 23 +344 0 1 +1645 4 4 +4093 25 25 +2446 1 0 +6329 1 0 +769 0 1 +873 7 0 +1299 0 1 +8250 15 15 +4291 10626 12548 +1063 103 0 +75 0 218 +5897 1 1 +33 1 1 +9905 51 37 +1711 1 1 +45 1 1 +1597 4 4 +61 1 1 +869 13 0 +1633 15 15 +1179 3 0 +2819 0 2 +156 88 465 +100 0 25 +4753 0 8 +4903 6 6 +3847 0 2 +108126 0 145 +3542 1 0 +10276 0 1 +402 1 1 +47 1 1 +63 2 0 +583 1 1 +34 1 1 +8974 0 46 +379 1 1 +70 1 1 +2474 0 102 +1756 0 1 +8426 0 1 +11986570 151150 142000 +14207 1 1 +49 1 1 +1781 0 12 +409 1 0 +1874 1 1 +33 1 1 +4795 2 0 +3686 1 1 +15 1 1 +164 2 0 +7287 5 5 +702023 72796 150000 +69 1 0 +1598 0 1 +1340 0 1 +1030 0 4 +1166 10 10 +6406 1 1 +34 1 1 +3283 2 4 +496 1 23 +25 10 0 +40 1 1 +21282 2 1 +1406 0 44 +937 17 0 +18255 1 0 +172366 0 281526 +40517 0 406 +58864 1 2 +1843 0 2 +14994 0 1 +12992 2 0 +2571146 1 0 +35688061 108224 100000 +29671719 1 0 +39665349 1 0 +14876089 64197 20000 +401 0 3 +1831 1 1 +61 1 1 +654 11 1 +127 0 28 +116 25 29 +74 9 9 +3247 2 0 +1355 1 0 +816 1 0 +21 3 0 +526 1 1 +60 1 1 +773 1 0 +1494 1 0 +1683 1 0 +21 1 1 +596 1 0 +1582 1 0 +425 0 3 +1139 0 8 +138 144 813 +71 1 0 +385638 1 1 +219 1 1 +51 7 7 +54 15205 15205 +54 7 7 +51 1 1 +110 1 1 +5198249 0 1 +115 80 0 +60 40 0 +80 100 20 +442 1 1 +35 1 1 +119 9 129 +14 0 40 +61 1 1 +85 0 40 +1331 0 1 +1194 12 14 +1773 0 5 +1623 1 1 +25 1 1 +2422 21 22 +2315 0 1 +31 1 1 +1508 0 3 +50 0 4 +970 6 6 +7234 4 0 +1228 1 0 +1233 8 1 +4207 1 1 +41 1 1 +1715 50 50 +1532 0 2 +1714 2 0 +2299 2 0 +2143 1 1 +18 1 1 +917 1 0 +1082 1 1 +49 1 1 +9098 0 4 +5529 0 1 +3266 0 8 +894 6 0 +1440 1 1 +2095 2 0 +4179 18 13 +1785 0 1 +1674 0 6 +1217 32348 33000 +287 17 23 +40 173 0 +20 161 0 +13 66 0 +12 71 1 +12548 30000 30000 +952154 41011 25000 +1914 1 0 +225 1 1 +57 1 1 +2910 0 1 +471 1 0 +1956 1 1 +43 1 1 +991 0 92 +2863 3 0 +1838 306 0 +1006 30 0 +206 1 1 +27 0 3 +221 17 17 +1525 8 0 +413 1 1 +32 2 2 +1573 12 12 +2258676 + +chain 18404394637 chr3 198022430 + 60000 197962430 chr3 199501827 + 35000 199446827 3 +13552792 0 1 +585083 0 1 +2092 1 2 +71 0 1 +31067262 0 1 +90 1 0 +5628461 0 15 +788 1 0 +3676 0 2 +35 1 1 +4682 1 0 +1295 0 5 +437 53 42 +326 0 2 +3309 1 1 +40 1 1 +998 0 3 +3584 5 0 +930 5 7 +845 1 1 +28 2 1 +303 0 1 +507 0 5 +5088 1 0 +2185 10046 10048 +949 4 4 +35 1 1 +2874 0 1 +6717 0 1 +2823 1 0 +703 0 4 +5143 1 0 +2490 10 10 +162 7 7 +6308 0 15 +15177535 152352 260003 +2532 1 0 +1633 2 2 +24207544 3000000 4400000 +60453436 0 4 +40076813 20999 20000 +8195 1 1 +44 1 1 +1408 1 1 +33 1 1 +172 3 4 +25 1 1 +538 13 13 +173 1 1 +27 1 1 +636 1 1 +28 1 1 +470 3 3 +112 1 1 +375 1 1 +26 26 0 +115 78 0 +52 23 0 +104 26 0 +15 231 0 +56 1 1 +20 1 1 +815 16 16 +484 1 1 +49 18 0 +509 4 3 +23 1 1 +391 1 1 +44 1 1 +1278 0 4 +32 1 1 +263 4 0 +27 1 1 +2416 0 2 +3184 43 43 +1150 13 13 +736 1 1 +47 1 1 +1416 1 1 +20 1 1 +1020 1 1 +25 1 1 +482 5 0 +716 1 0 +1697 1 1 +29 1 1 +1270261 23108 27000 +40930 1 0 +17455 1 1 +28 4 4 +43 0 491 +1679 1 1 +33 7 5 +29 1 3 +54863 53 16 +7653 2 2 +29 2 2 +15 0 30 +180 0 75 +215 45 0 +360 915 0 +6188 96 0 +17 288 0 +1388 1 1 +37 1 1 +7956 0 1 +374 0 1 +2442058 + +chain 310844 chr4 191154276 + 9241835 9252105 chr4 191273063 + 8941714 8947231 850 +2726 14 14 +243 9 9 +78 4 4 +157 213 211 +265 85 85 +524 4751 0 +379 283 283 +539 + +chain 9105 chr4 191154276 + 8982741 8983189 chr10 135374737 - 120304868 120305301 22041837 +59 380 365 +9
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bam_01_input_a.sam Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,27 @@ +@HD VN:1.0 SO:coordinate +@SQ SN:chr1 LN:1000051 +@SQ SN:chr2 LN:1000051 +@SQ SN:chr3 LN:1000051 +@SQ SN:chr4 LN:9250051 +@PG ID:- VN:1.0.0 CL:cmatrix +@CO Test data for CrossMap bam +read_001 0 chr1 100000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_002 0 chr2 100000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_003 0 chr3 100000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_004 0 chr4 9200000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_005 0 chr4 8940000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_006 0 chr1 1000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_007 0 chr2 1000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_008 0 chr3 1000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_009 0 chr4 9250000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_010 0 chr4 9000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_011 16 chr1 100000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_012 16 chr2 100000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_013 16 chr3 100000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_014 16 chr4 9200000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_015 16 chr4 8940000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_016 16 chr1 1000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_017 16 chr2 1000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_018 16 chr3 1000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_019 16 chr4 9250000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_020 16 chr4 9000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN *
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bam_01_output_a.sam Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,23 @@ +@HD VN:1.0 SO:coordinate +@SQ SN:chr1 LN:247249719 +@SQ SN:chr10 LN:135374737 +@SQ SN:chr2 LN:242951149 +@SQ SN:chr3 LN:199501827 +@SQ SN:chr4 LN:191273063 +@PG ID:- VN:1.0.0 CL:cmatrix +@PG ID:CrossMap VN:0.2 +@CO Test data for CrossMap bam +@CO Liftover from original BAM/SAM file: +@CO Liftover is based on the chain file: +read_001 0 chr1 89863 60 50M = 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_002 0 chr2 90000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_003 0 chr3 75000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_006 0 chr1 989863 60 50M = 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_007 0 chr2 990000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_008 0 chr3 975000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_011 16 chr1 89863 60 50M = 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_012 16 chr2 90000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_013 16 chr3 75000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_016 16 chr1 989863 60 50M = 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_017 16 chr2 990000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_018 16 chr3 975000 60 50M chr1 1 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN *
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bam_01_output_a.unmap.sam Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,15 @@ +@HD VN:1.0 SO:coordinate +@SQ SN:chr1 LN:1000051 +@SQ SN:chr2 LN:1000051 +@SQ SN:chr3 LN:1000051 +@SQ SN:chr4 LN:9250051 +@PG ID:- VN:1.0.0 CL:cmatrix +@CO Test data for CrossMap bam +read_004 0 chr4 9200000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_005 0 chr4 8940000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_009 0 chr4 9250000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_010 0 chr4 9000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_014 16 chr4 9200000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_015 16 chr4 8940000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_019 16 chr4 9250000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN * +read_020 16 chr4 9000000 60 50M * 0 0 NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN *
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bed_01_input_a.bed Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,5 @@ +chr1 100000 1000000 +chr2 100000 1000000 +chr3 100000 1000000 +chr4 9200000 9250000 +chr4 8940000 9000000
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bed_01_output_a__all.bed Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,14 @@ +chr1 100000 1000000 (split.1:chr1:100000:177417:+) chr1 89863 167280 +chr1 100000 1000000 (split.2:chr1:227417:267719:+) chr1 217280 257582 +chr1 100000 1000000 (split.3:chr1:317719:471368:+) chr1 307582 461231 +chr1 100000 1000000 (split.4:chr1:521368:1000000:+) chr1 511231 989863 +chr2 100000 1000000 -> chr2 90000 990000 +chr3 100000 1000000 -> chr3 75000 975000 +chr4 9200000 9250000 (split.1:chr4:9241835:9244561:+) chr4 8941714 8944440 +chr4 9200000 9250000 (split.2:chr4:9244575:9244818:+) chr4 8944454 8944697 +chr4 9200000 9250000 (split.3:chr4:9244827:9244905:+) chr4 8944706 8944784 +chr4 9200000 9250000 (split.4:chr4:9244909:9245066:+) chr4 8944788 8944945 +chr4 9200000 9250000 (split.5:chr4:9245279:9245544:+) chr4 8945156 8945421 +chr4 9200000 9250000 (split.6:chr4:9245629:9246153:+) chr4 8945506 8946030 +chr4 8940000 9000000 (split.1:chr4:8982741:8982800:+) chr10 15069810 15069869 +chr4 8940000 9000000 (split.2:chr4:8983180:8983189:+) chr10 15069436 15069445
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bed_01_output_a__only-matches.bed Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,14 @@ +chr1 89863 167280 +chr1 217280 257582 +chr1 307582 461231 +chr1 511231 989863 +chr2 90000 990000 +chr3 75000 975000 +chr4 8941714 8944440 +chr4 8944454 8944697 +chr4 8944706 8944784 +chr4 8944788 8944945 +chr4 8945156 8945421 +chr4 8945506 8946030 +chr10 15069810 15069869 +chr10 15069436 15069445
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bed_02_input_a.bed Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,5 @@ +chr1 100 10000 +chr2 100 10000 +chr3 100 10000 +chr4 8941700 8947200 +chr5 1 100000000 \ No newline at end of file
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bed_02_output_a__all.bed Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,5 @@ +chr1 100 10000 Fail +chr2 100 10000 Fail +chr3 100 10000 Fail +chr4 8941700 8947200 Fail +chr5 1 100000000 Fail
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_bigwig_01_output_a.sorted.bgr Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,3 @@ +chr1 889863 989862 110.0 +chr2 890000 989999 220.0 +chr3 875000 974999 330.0
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_gff_01_input_a.gtf Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,5 @@ +chr1 CrossMap_test_data CDS 100000 1000000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr2 CrossMap_test_data CDS 100000 1000000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr3 CrossMap_test_data CDS 100000 1000000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr4 CrossMap_test_data CDS 9200000 9250000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr5 CrossMap_test_data CDS 8940000 9000000 0.000000 - 0 gene_id "TEST_GENE_1";
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_gff_01_output_a__all.gtf Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,5 @@ +chr1 CrossMap_test_data CDS 100000 1000000 0.000000 - 0 gene_id "TEST_GENE_1"; fail (multpile match to target assembly) +chr2 CrossMap_test_data CDS 100000 1000000 0.000000 - 0 gene_id "TEST_GENE_1"; -> chr2 CrossMap_test_data CDS 90000 990000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr3 CrossMap_test_data CDS 100000 1000000 0.000000 - 0 gene_id "TEST_GENE_1"; -> chr3 CrossMap_test_data CDS 75000 975000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr4 CrossMap_test_data CDS 9200000 9250000 0.000000 - 0 gene_id "TEST_GENE_1"; fail (multpile match to target assembly) +chr5 CrossMap_test_data CDS 8940000 9000000 0.000000 - 0 gene_id "TEST_GENE_1"; fail (no match to target assembly)
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_gff_01_output_a__only-matches.gtf Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,2 @@ +chr2 CrossMap_test_data CDS 90000 990000 0.000000 - 0 gene_id "TEST_GENE_1"; +chr3 CrossMap_test_data CDS 75000 975000 0.000000 - 0 gene_id "TEST_GENE_1";
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_vcf_01.fasta Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,10 @@ +>chr1 +CAAAAAAGCAGTTGACAGTTGTTCGGGCTTGTTACGCATCTTAAGATCGAATAAGGAGAGAGGGTCTAAACG +AGGCGCACCCCGCCTATGGGTGATCGAGGTACTAGGGGTTGGTCGCAGGTCTGTATTACCGTTAGCGGTGCA +AGGGGATCTGATCGAGTGATTCACCTACTCATGTGGCGAGCACGCCGACGAAATACTCCTGGTCGTGTTATA +AAGCCCTGGTTTTCCTTTCC +>chr2 +CACAAAATGCACGTGGATGCAGGCATTTATCCAACCCACACTATTACGTTCACCAAATGTGTGGACCAACTG +CGGGACTAGGTAAGCTTGTCCTCAATGAGCGAAATTGATATTTCTCTACCGACTTGGGGTCGACTGGACGAG +TCAGCTGTGCAACAGCTCAGCCGGTTTCGATAAACCGAAACCTTGAATGTTTGGACTTGCGTCATGGCGAAC +AAAGATCGTTCATGTCGCA \ No newline at end of file
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_vcf_01.over.chain Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,11 @@ +chain 4900 chr1 236 + 1 236 chr2 236 - 1 237 1 + 9 1 0 + 10 0 5 + 61 4 0 + 16 0 4 + 42 3 0 + 16 0 8 + 14 1 0 + 3 7 0 + 48 +
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_vcf_01_input.vcf Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,8 @@ +##fileformat=VCFv4.2 +##FORMAT=<ID=GT,Number=1,Type=Integer,Description="Genotype"> +##FORMAT=<ID=GP,Number=G,Type=Float,Description="Genotype Probabilities"> +##FORMAT=<ID=PL,Number=G,Type=Float,Description="Phred-scaled Genotype Likelihoods"> +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SAMP001 SAMP002 +1 10 rs11449 G A . PASS . GT 0/0 0/1 +1 100 rs84825 C C . PASS . GT:GP 0/1:. 0/1:0.03,0.97,0 +1 200 rs84823 T G . PASS . GT:PL ./.:. 1/1:10,5,0
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_vcf_01_output.vcf Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,12 @@ +##fileformat=VCFv4.2 +##FORMAT=<ID=GT,Number=1,Type=Integer,Description="Genotype"> +##FORMAT=<ID=GP,Number=G,Type=Float,Description="Genotype Probabilities"> +##FORMAT=<ID=PL,Number=G,Type=Float,Description="Phred-scaled Genotype Likelihoods"> +##liftOverProgram=CrossMap(https://sourceforge.net/projects/crossmap/) +##liftOverFile=test_vcf_01.over.chain +##new_reference_genome=test_vcf_01.fasta +##liftOverTime=August13,2015 +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SAMP001 SAMP002 +2 227 rs11449 C A . PASS . GT 0/0 0/1 +2 137 rs84825 T C . PASS . GT:GP 0/1:. 0/1:0.03,0.97,0 +2 36 rs84823 C G . PASS . GT:PL ./.:. 1/1:10,5,0
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_vcf_01_output.vcf.unmap Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,5 @@ +##fileformat=VCFv4.2 +##FORMAT=<ID=GT,Number=1,Type=Integer,Description="Genotype"> +##FORMAT=<ID=GP,Number=G,Type=Float,Description="Genotype Probabilities"> +##FORMAT=<ID=PL,Number=G,Type=Float,Description="Phred-scaled Genotype Likelihoods"> +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SAMP001 SAMP002
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_wig_01_input_a.wig Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,10 @@ +variableStep chrom=chr1 span=900000 +100000 110 +variableStep chrom=chr2 span=900000 +100000 220 +variableStep chrom=chr3 span=900000 +100000 330 +variableStep chrom=chr4 span=50000 +9200000 400 +variableStep chrom=chr4 span=60000 +8940000 450
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_wig_01_output_a.sorted.bgr Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,2 @@ +chr2 89999 989999 220.0 +chr3 74999 974999 330.0
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/all_fasta.loc.sample Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,18 @@ +#This file lists the locations and dbkeys of all the fasta files +#under the "genome" directory (a directory that contains a directory +#for each build). The script extract_fasta.py will generate the file +#all_fasta.loc. This file has the format (white space characters are +#TAB characters): +# +#<unique_build_id> <dbkey> <display_name> <file_path> +# +#So, all_fasta.loc could look something like this: +# +#apiMel3 apiMel3 Honeybee (Apis mellifera): apiMel3 /path/to/genome/apiMel3/apiMel3.fa +#hg19canon hg19 Human (Homo sapiens): hg19 Canonical /path/to/genome/hg19/hg19canon.fa +#hg19full hg19 Human (Homo sapiens): hg19 Full /path/to/genome/hg19/hg19full.fa +# +#Your all_fasta.loc file should contain an entry for each individual +#fasta file. So there will be multiple fasta files for each build, +#such as with hg19 above. +#
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/liftOver.loc.sample Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,27 @@ +#This is a sample file distributed with Galaxy that is used by the +#liftOver tools. The liftOver.loc file has this format (white space +#characters are TAB characters): +# +#<FromSpecies> <ToSpecies> <PathToChainFile> +# +#So, for example, if you had the chain file to convert from anoCar1 to galGal3 +#located at /depot/data2/galaxy/anoCar1/liftOver/anoCar1ToGalGal3.over.chain, +#then the liftOver.loc entry would look like this: +# +#anoCar1 galGal3 /depot/data2/galaxy/anoCar1/liftOver/anoCar1ToGalGal3.over.chain +# +#and your /depot/data2/galaxy/anoCar1/liftOver directory would +#contain all of your "chain" files (e.g.): +# +#-rw-rw-r-- 1 gua110 galaxy 24046079 2008-01-16 14:20 anoCar1ToGalGal3.over.chain +#-rw-rw-r-- 1 gua110 galaxy 13216668 2008-01-16 14:20 anoCar1ToGasAcu1.over.chain +#-rw-rw-r-- 1 gua110 galaxy 29597067 2008-01-16 14:20 anoCar1ToHg18.over.chain +#...etc... +# +#Your liftOver.loc file should include an entry per line for each build you can +#convert. For example: +# +#anoCar1 galGal3 /depot/data2/galaxy/anoCar1/liftOver/anoCar1ToGalGal3.over.chain +#anoCar1 gasAcu1 /depot/data2/galaxy/anoCar1/liftOver/anoCar1ToGasAcu1.over.chain +#anoCar1 hg18 /depot/data2/galaxy/anoCar1/liftOver/anoCar1ToHg18.over.chain +#...etc...
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_dependencies.xml Wed Aug 26 11:04:47 2015 -0400 @@ -0,0 +1,9 @@ +<?xml version="1.0"?> +<tool_dependency> + <package name="ucsc_tools" version="312"><!-- For conversion of Wig and BigWig --> + <repository changeset_revision="2d6bafd63401" name="package_ucsc_tools_312" owner="iuc" toolshed="https://toolshed.g2.bx.psu.edu" /> + </package> + <package name="crossmap" version="0.2"> + <repository changeset_revision="805c8fbd3b4e" name="package_crossmap_0_2" owner="yhoogstrate" prior_installation_required="True" toolshed="https://toolshed.g2.bx.psu.edu" /> + </package> +</tool_dependency>