<tool id="t-vs-vnormal_junctions" name="Virtual Normal Correction SVs" version="1.8">
	<description> Filter SVs based on presence in VN set </description>
	
	<requirements>		
		<requirement type="package" version="1.7">cgatools</requirement>
	</requirements>
	
	<command interpreter="bash">
	JunctionDiff-vs-background.sh 
		--variants $variants 
		--reference ${reference.fields.reference_crr_cgatools} 
		--VN_junctions ${reference.fields.VN_genomes_junctionfile_list}${VNset}
		--cgatools_binary cgatools
		--outputfile_filtered $output_filtered 
		--scoreThresholdA $scoreThresholdA 
		--scoreThresholdB $scoreThresholdB 
		--distance $distance 
		--minlength $minlength			
	</command>
	
	<inputs>
		<param name="variants" type="data" format="tabular" label="CG Junctions file"/>	
		
		<!--select build-->
		<param name="reference" type="select" label="Select Build">
			<options from_data_table="virtual_normal_correction" />
		</param>

		<param name="VNset" type="select" label="Select Virtual Normal set to use" help="1000Genomes set can only be used for hg19 samples, for hg18 54 genomes will be used.">
			<option value="46_diversity.txt" > CG Diversity Panel (46 Genomes) </option>
			<option value="433_1000g.txt" > CG 1000G project genomes (433 Genomes) (hg19 only) </option>
			<option value="479_diversity_1000g.txt" > Diversity and 1000G (479 genomes) (hg19 only) </option>
                        <option value="10_tutorial.txt" > Small VN for tutorial (10 Genomes) </option>
		</param>
		
			
		<param name="scoreThresholdA" type="text" value="10" label="scoreThreshold" help="The minimum number of discordant mate pair alignments supporting the junction from input genome"/>
		<param name="scoreThresholdB" type="text" value="10" label="scoreThreshold" help="The minimum number of discordant mate pair alignments supporting the junction from background genomes"/>
		<param name="distance" type="text" value="200" label="Maximum distance between coordinates of potentially compatible junctions."/>		
		<param name="minlength" type="text" value="500" label="Minimum deletion junctions length to be included into the difference file."/>
		<param name="report" type="select"  label="Generate report file?">
			<option value="N" selected="true"> No </option>
			<option value="Y"> Yes </option>
		</param>
		<param name="fname" type="text" value="" label="Prefix for your output file" help="Optional"/>		
	</inputs>

  <outputs>
       <data format="tabular" name="output_filtered" label="${fname} Filtered junctions for ${tool.name} on  ${on_string}"/>		
	   <data format="tabular" name="output_report" from_work_dir= "output_reports.tsv" label="${fname} report for ${tool.name} on  ${on_string}">
		   <filter> report == "Y" </filter>
	   </data>	
  </outputs>

	<help> 
**What it does**



**Input Files**
Complete Genomics Junctions file

**Output Files**
Junctions remaining after filtering


	</help>

</tool>