# HG changeset patch # User iuc # Date 1700518663 0 # Node ID 6be888be75f9aa75aa02ea984d12c9892c7d5deb # Parent 5826298f6a73c6b22be40e44cbdd77ab8a0fe284 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tool_collections/samtools/samtools_view commit cd62639660bef041ba14ecff337fb98e84e75d8a diff -r 5826298f6a73 -r 6be888be75f9 macros.xml --- a/macros.xml Mon Aug 15 09:19:43 2022 +0000 +++ b/macros.xml Mon Nov 20 22:17:43 2023 +0000 @@ -5,7 +5,14 @@ + 1.15.1 + 2 20.05 - - @misc{SAM_def, - title={Definition of SAM/BAM format}, - url = {https://samtools.github.io/hts-specs/},} - - 10.1093/bioinformatics/btp352 - 10.1093/bioinformatics/btr076 - 10.1093/bioinformatics/btr509 - - @misc{Danecek_et_al, - Author={Danecek, P., Schiffels, S., Durbin, R.}, - title={Multiallelic calling model in bcftools (-m)}, - url = {http://samtools.github.io/bcftools/call-m.pdf},} - - - @misc{Durbin_VCQC, - Author={Durbin, R.}, - title={Segregation based metric for variant call QC}, - url = {http://samtools.github.io/bcftools/rd-SegBias.pdf},} - - - @misc{Li_SamMath, - Author={Li, H.}, - title={Mathematical Notes on SAMtools Algorithms}, - url = {http://www.broadinstitute.org/gatk/media/docs/Samtools.pdf},} - - - @misc{SamTools_github, - title={SAMTools GitHub page}, - url = {https://github.com/samtools/samtools},} - + 10.1093/gigascience/giab008 diff -r 5826298f6a73 -r 6be888be75f9 samtools_view.xml --- a/samtools_view.xml Mon Aug 15 09:19:43 2022 +0000 +++ b/samtools_view.xml Mon Nov 20 22:17:43 2023 +0000 @@ -1,4 +1,4 @@ - + - reformat, filter, or subsample SAM, BAM or CRAM macros.xml @@ -154,11 +154,11 @@ ## not dealing with all of the reads in the indexed ## file. We have to do an extra pass over the input to ## count the reads to subsample. - sample_fragment=`samtools view -c $std_filters infile $reg_filters | awk '{s=\$1} END {frac=s/${mode.subsample_config.subsampling_mode.target}; printf("%.8f\n", frac > 1 ? $seed+1/frac : ".0")}'` && + sample_fragment=`samtools view -c $std_filters infile $reg_filters | awk '{s=\$1} END {fac=s/${mode.subsample_config.subsampling_mode.target}; printf("%.8f\n", fac > 1 ? 1/fac : 1)}'` && #else: ## We can get the count of reads to subsample using ## an inexpensive call to idxstats. - sample_fragment=`samtools idxstats infile | awk '{s+=\$4+\$3} END {frac=s/${mode.subsample_config.subsampling_mode.target}; printf("%.8f\n", frac > 1 ? $seed+1/frac : ".0")}'` && + sample_fragment=`samtools idxstats infile | awk '{s+=\$4+\$3} END {fac=s/${mode.subsample_config.subsampling_mode.target}; printf("%.8f\n", fac > 1 ? 1/fac : 1)}'` && #end if #end if #end if @@ -172,12 +172,13 @@ $std_filters #if $with_subsampling: + --subsample-seed $seed #if str($mode.subsample_config.subsampling_mode.select_subsample) == "target": ##this is calculated at execution time before the main samtools command - -s \${sample_fragment} + --subsample \${sample_fragment} #else: - #set $fraction = $seed + 1 / float($mode.subsample_config.subsampling_mode.factor) - -s $fraction + #set $fraction = 1 / float($mode.subsample_config.subsampling_mode.factor) + --subsample $fraction #end if #end if @@ -398,13 +399,13 @@ - - + + - + @@ -413,7 +414,7 @@ - + @@ -422,7 +423,7 @@ - + @@ -431,7 +432,7 @@ - + @@ -440,7 +441,7 @@ - + @@ -456,7 +457,7 @@ - + @@ -465,7 +466,7 @@ - + @@ -481,7 +482,7 @@ - + @@ -490,7 +491,7 @@ - + @@ -506,7 +507,7 @@ - + @@ -519,7 +520,7 @@ - + @@ -532,7 +533,7 @@ - + @@ -545,7 +546,7 @@ - + @@ -556,7 +557,7 @@ - + @@ -578,7 +579,7 @@ - + @@ -600,7 +601,7 @@ - + @@ -622,7 +623,7 @@ - + @@ -645,7 +646,7 @@ - + @@ -668,7 +669,7 @@ - + @@ -691,7 +692,7 @@ - + @@ -710,7 +711,7 @@ - + @@ -729,7 +730,7 @@ - + @@ -749,7 +750,7 @@ - + @@ -769,7 +770,7 @@ - + @@ -789,7 +790,7 @@ - + @@ -809,7 +810,7 @@ - + @@ -829,7 +830,7 @@ - + @@ -849,7 +850,7 @@ - + @@ -870,7 +871,7 @@ - + @@ -889,7 +890,7 @@ - +