Repository 'bcftools_consensus'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/bcftools_consensus

Changeset 17:95c6f3f4c77f (2022-09-20)
Previous changeset 16:d37e313a41be (2022-07-23) Next changeset 18:f38ad709f79f (2022-09-26)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bcftools commit db275932cbb485cb44ae91e0b421d6f57698db49
modified:
bcftools_consensus.xml
macros.xml
test-data/summary.pdf
added:
test-data/mpileup.AD.vcf
b
diff -r d37e313a41be -r 95c6f3f4c77f bcftools_consensus.xml
--- a/bcftools_consensus.xml Sat Jul 23 13:32:35 2022 +0000
+++ b/bcftools_consensus.xml Tue Sep 20 12:48:09 2022 +0000
[
@@ -1,5 +1,5 @@
 <?xml version='1.0' encoding='utf-8'?>
-<tool name="bcftools @EXECUTABLE@" id="bcftools_@EXECUTABLE@" version="@TOOL_VERSION@+galaxy1">
+<tool name="bcftools @EXECUTABLE@" id="bcftools_@EXECUTABLE@" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
     <description>Create consensus sequence by applying VCF variants to a reference fasta file</description>
     <macros>
         <token name="@EXECUTABLE@">consensus</token>
@@ -30,6 +30,18 @@
   --mask '${section.mask}'
 #end if
 
+#if $section.mark_del
+    --mark-del '$section.mark_del'
+#end if
+
+#if $section.mark_ins
+    --mark-ins $section.mark_ins
+#end if
+
+#if $section.mark_snv
+    --mark-snv $section.mark_snv
+#end if
+
 #if $section.select_haplotype:
   --haplotype '${section.select_haplotype}'
 #end if
@@ -40,7 +52,11 @@
 @EXCLUDE@
 
 #if $chain:
-  --chain '$chain_file'
+    --chain '$chain_file'
+#end if
+
+#if $absent
+    --absent '$absent'
 #end if
 
 ## Primary Input/Outputs
@@ -66,9 +82,51 @@
                 <option value="1">1</option>
                 <option value="2">2</option>
             </param>
+            <param argument="--mark-del" type="text" value="" optional="true" label="Mark deletions" help="Instead of removing sequence, insert CHAR for deletions">
+                <sanitizer invalid_char="">
+                    <valid initial="string.letters,string.digits">
+                        <add value="_" />
+                    </valid>
+                </sanitizer>
+                <validator type="regex">[0-9a-zA-Z_]+</validator>
+            </param>
+            <param argument="--mark-ins" type="select" optional="true" label="Mark insertions" help="Highlight insertions in uppercase (uc) or lowercase (lc), leaving the rest as is">
+                <option value="uc">Uppercase</option>
+                <option value="lc">Lowercase</option>
+            </param>
+            <param argument="--mark-snv" type="select" optional="true" label="Mark substitutions" help="Highlight substitutions in uppercase (uc) or lowercase (lc), leaving the rest as is">
+                <option value="uc">Uppercase</option>
+                <option value="lc">Lowercase</option>
+            </param>
+            <conditional name="conditional_mask">
+                <param name="selector" type="select" label="Mask file option">
+                    <option value="disabled">Disabled</option>
+                    <option value="enabled">Enabled</option>
+                </param>
+                <when value="disabled"/>
+                <when value="enabled">
+                    <param argument="--mask" type="data" format="tabular" label="Mask" help="Replace regions according to the next --mask-with option" />
+                    <param argument="--mask-with" type="text" value="N" optional="true" label="Mask with" help="Replace with CHAR (skips overlapping variants); change to uppercase (uc) or lowercase (lc)">
+                        <sanitizer invalid_char="">
+                            <valid initial="string.letters,string.digits">
+                                <add value="_" />
+                            </valid>
+                        </sanitizer>
+                        <validator type="regex">[0-9a-zA-Z_]+</validator>
+                    </param>
+                </when>
+            </conditional>
         </section>
         <param name="chain" type="boolean" truevalue="yes" falsevalue="no" label="Write a chain file for liftover" />
         <param name="rename" type="boolean" truevalue="yes" falsevalue="no" label="Set output FASTA ID from name of VCF" />
+        <param argument="--absent" type="text" value="" label="Absent" optional="true" help="It allows to set positions with no supporting evidence to N (or any other character)">
+            <sanitizer invalid_char="">
+                <valid initial="string.letters,string.digits,string.punctuation">
+                    <remove value="@" />
+                    <remove value="'" />
+                </valid>
+            </sanitizer>
+        </param>
         <section name="sec_restrict" expanded="false" title="Restrict to">    
             <expand macro="macro_include" />
             <expand macro="macro_exclude" />
@@ -140,6 +198,46 @@
                 </assert_contents>
             </output>
         </test>
+        <!--Test absent option-->
+        <test>
+            <expand macro="test_using_reference" ref="consensus.fa" />
+            <param name="input_file" ftype="vcf" value="consensus.vcf" />
+            <section name="sec_restrict">
+                <param name="include" value='TYPE="snp"' />
+            </section>
+            <param name="absent" value="W"/>
+            <output name="output_file">
+                <assert_contents>
+                    <has_text text="WWWAWAWWAWWWWWWWWCWWWWWWWW" />
+                </assert_contents>
+            </output>
+            <assert_command>
+                <has_text text="--absent" />
+            </assert_command>
+        </test>
+        <!--Test mask options -->
+        <test>
+            <expand macro="test_using_reference" ref="consensus.fa" />
+            <param name="input_file" ftype="vcf" value="consensus.vcf" />
+            <section name="sec_restrict">
+                <param name="include" value='TYPE="snp"' />
+            </section>
+            <section name="sec_default">
+                <param name="mark_del" value="DEL"/>
+                <param name="mark_ins" value="uc"/>
+                <param name="mark_snv" value="uc"/>
+            </section>
+            <output name="output_file">
+                <assert_contents>
+                    <has_text text="TACAAAATATGACATATCAAAAAGAACATAACCTACGTATCAACTAAAGTGGTTGTTTGA" />
+                </assert_contents>
+            </output>
+            <assert_command>
+                <has_text text="--mark-del" />
+                <has_text text="--mark-ins" />
+                <has_text text="--mark-snv" />
+            </assert_command>
+        </test>
     </tests>
     <help><![CDATA[
 =====================================
b
diff -r d37e313a41be -r 95c6f3f4c77f macros.xml
--- a/macros.xml Sat Jul 23 13:32:35 2022 +0000
+++ b/macros.xml Tue Sep 20 12:48:09 2022 +0000
[
@@ -1,5 +1,7 @@
 <macros>
-  <token name="@TOOL_VERSION@">1.10</token>
+  <token name="@TOOL_VERSION@">1.15.1</token>
+  <token name="@VERSION_SUFFIX@">0</token>
+  <token name="@PROFILE@">20.01</token>
   <xml name="bio_tools">
       <xrefs>
           <xref type="bio.tools">bcftools</xref>
@@ -8,15 +10,15 @@
   <xml name="requirements">
     <requirements>
       <requirement type="package" version="@TOOL_VERSION@">bcftools</requirement>
-      <requirement type="package" version="1.10">htslib</requirement>
+      <requirement type="package" version="1.15.1">htslib</requirement>
       <yield />
     </requirements>
   </xml>
   <xml name="samtools_requirement">
-      <requirement type="package" version="1.10">samtools</requirement>
+      <requirement type="package" version="1.15.1">samtools</requirement>
   </xml>
   <xml name="matplotlib_requirement">
-      <requirement type="package" version="3.4.3">matplotlib</requirement>
+      <requirement type="package" version="3.5.3">matplotlib</requirement>
   </xml>
   <xml name="version_command">
     <version_command>bcftools 2&gt;&amp;1 | grep 'Version:'</version_command>
@@ -39,6 +41,10 @@
 ]]>
   </token>
   <xml name="macro_input">
+    <!-- 
+    REQUIRES https://github.com/galaxyproject/galaxy/pull/14605/files
+    <param name="input_files" type="data" format="vcf,vcf.gz,vcf_bgzip,bcf" label="Other VCF/BCF Datasets" multiple="True" />
+    -->
     <param name="input_file" type="data" format="vcf,vcf_bgzip,bcf" label="VCF/BCF Data" />
   </xml>
   <token name="@PREPARE_INPUT_FILE@">
@@ -48,6 +54,7 @@
 #if $input_file.is_of_type('vcf')
   bgzip -c '$input_file' > $input_vcf &&
   bcftools index $input_vcf &&
+##elif $input_file.is_of_type('vcf_bgzip') or $input_file.is_of_type('vcf.gz')
 #elif $input_file.is_of_type('vcf_bgzip')
   ln -s '$input_file' $input_vcf &&
   #if $input_file.metadata.tabix_index:
@@ -71,6 +78,10 @@
   </token>
 
   <xml name="macro_inputs">
+    <!-- 
+    REQUIRES https://github.com/galaxyproject/galaxy/pull/14605/files
+    <param name="input_files" type="data" format="vcf,vcf.gz,vcf_bgzip,bcf" label="Other VCF/BCF Datasets" multiple="True" />
+    -->
     <param name="input_files" type="data" format="vcf,vcf_bgzip,bcf" label="Other VCF/BCF Datasets" multiple="True" />
   </xml>
   <token name="@PREPARE_INPUT_FILES@">
@@ -83,7 +94,8 @@
   #if $input_file.is_of_type('vcf')
     bgzip -c '$input_file' > $input_vcf &&
     bcftools index $input_vcf &&
-  #elif $input_file.is_of_type('vcf_bgzip')
+  ##elif $input_file.is_of_type('vcf_bgzip') or $input_file.is_of_type('vcf.gz') -> REQUIRES  https://github.com/galaxyproject/galaxy/pull/14605
+  #elif $input_file.is_of_type('vcf_bgzip') or $input_file.is_of_type('vcf.gz')
     ln -s '$input_file' $input_vcf &&
     #if $input_file.metadata.tabix_index:
       ln -s '${input_file.metadata.tabix_index}' ${input_vcf}.tbi &&
@@ -318,6 +330,11 @@
             <yield />
         </when>        
     </conditional>
+    <param argument="--@TYPE@s-overlap" type="select" optional="true" label="@LABEL_TYPE@ overlap" help="Include if POS in the region (0), record overlaps (1), variant overlaps (2)">
+        <option value="0">0: POS in the region</option>
+        <option value="1">1: Record overlaps</option>
+        <option value="2">2: Variant overlaps</option>
+    </param>
   </xml>
 
   <token name="@PARSE_INTERVALS@">
@@ -337,6 +354,28 @@
 ]]>
   </token>
   
+  <token name="@MASK@">
+<![CDATA[
+#if $section.conditional_soft_filter.selector == 'enabled' and $section.conditional_soft_filter.soft_filter
+  #if $section.conditional_soft_filter.masks.masks_src == 'regions':
+    #set $intervals = $section.conditional_soft_filter.masks.masks
+    @PARSE_INTERVALS@
+    --mask '$intervals_spec'
+  #elif $section.conditional_soft_filter.masks.masks_src == 'masks_file' and $section.conditional_soft_filter.masks.masks_file:
+    #if $masks_path is not None:
+      --mask-file '$masks_path'
+    #else:
+      --mask-file '$section.conditional_soft_filter.masks.masks_file'
+    #end if
+  #end if
+  #if $section.conditional_soft_filter.masks_overlap
+    --mask-overlap $section.conditional_soft_filter.masks_overlap
+  #end if
+#end if
+
+]]>
+  </token>
+
   <token name="@REGIONS@">
 <![CDATA[
 #if $section.regions.regions_src == 'regions':
@@ -350,6 +389,10 @@
     --regions-file '$section.regions.regions_file'
   #end if
 #end if
+#if $section.regions_overlap
+  --regions-overlap $section.regions_overlap
+#end if
+
 ]]>
   </token>
   
@@ -364,6 +407,9 @@
 #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
+#if $section.targets_overlap
+  --targets-overlap $section.targets_overlap
+#end if
 ]]>
   </token>
 
@@ -410,8 +456,8 @@
 
   <xml name="macro_samples">
       <param argument="--samples" type="text" value="" optional="true" label="Samples"
-             help="Comma separated list of samples to annotate (or exclude)">
-          <validator type="regex" message="">^(\w+(,\w+)*)?$</validator>
+             help="Comma-separated list of samples to annotate (or exclude) or - to include all samples">
+          <validator type="regex" message="Comma-separated list of samples or - to include all samples">^(-|\w+(,\w+)*)?$</validator>
       </param>
       <param name="invert_samples" type="boolean" truevalue="^" falsevalue="" checked="false" label="Invert Samples"
              help="Inverts the query/filtering applied by Samples (adds &quot;^&quot; prefix to exclude)" />
@@ -476,7 +522,7 @@
   <xml name="macro_columns">
     <param name="columns" type="text" value="" optional="true" label="Columns"
             help="List of columns in the annotation file, e.g. CHROM,POS,REF,ALT,-,INFO/TAG. See man page for details">
-        <validator type="regex" message="COLUMN names  separated by commas">^([^,]+(,[^,]+)*)?$</validator>
+        <validator type="regex" message="COLUMN names separated by commas">^([^,]+(,[^,]+)*)?$</validator>
     </param>
   </xml>
   <token name="@COLUMNS@">
@@ -496,7 +542,18 @@
 ${section.vcf_ids}
   </token>
 
+<xml name="macro_output_tags">
+  <param name="output_tags" argument="--annotate" type="select" optional="true" multiple="True" display="checkboxes" label="Optional tags to output" help="--annotate">
+      <yield />
+  </param>
+</xml>
+
+<xml name="macro_overlap" token_argument="" token_label="">
+</xml>
+
+
   <token name="@OUTPUT_HELP@">
+
       <![CDATA[
 Output Type
 -----------
b
diff -r d37e313a41be -r 95c6f3f4c77f test-data/mpileup.AD.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mpileup.AD.vcf Tue Sep 20 12:48:09 2022 +0000
b
b'@@ -0,0 +1,4131 @@\n+##fileformat=VCFv4.2\n+##FILTER=<ID=PASS,Description="All filters passed">\n+##bcftoolsVersion=1.11-21-g0987715+htslib-1.11-3-g450c912\n+##bcftoolsCommand=mpileup -O v -a AD,DP,INFO/ADR,INFO/ADF,DP4,INFO/DPR -f xxx//mpileup.ref.fa -b xxx//bamlist.txt\n+##reference=file://mpileup.ref.fa\n+##contig=<ID=17,length=4200>\n+##ALT=<ID=*,Description="Represents allele(s) other than observed.">\n+##INFO=<ID=INDEL,Number=0,Type=Flag,Description="Indicates that the variant is an INDEL.">\n+##INFO=<ID=IDV,Number=1,Type=Integer,Description="Maximum number of raw reads supporting an indel">\n+##INFO=<ID=IMF,Number=1,Type=Float,Description="Maximum fraction of raw reads supporting an indel">\n+##INFO=<ID=DP,Number=1,Type=Integer,Description="Raw read depth">\n+##INFO=<ID=VDB,Number=1,Type=Float,Description="Variant Distance Bias for filtering splice-site artefacts in RNA-seq data (bigger is better)",Version="3">\n+##INFO=<ID=RPB,Number=1,Type=Float,Description="Mann-Whitney U test of Read Position Bias (bigger is better)">\n+##INFO=<ID=MQB,Number=1,Type=Float,Description="Mann-Whitney U test of Mapping Quality Bias (bigger is better)">\n+##INFO=<ID=BQB,Number=1,Type=Float,Description="Mann-Whitney U test of Base Quality Bias (bigger is better)">\n+##INFO=<ID=MQSB,Number=1,Type=Float,Description="Mann-Whitney U test of Mapping Quality vs Strand Bias (bigger is better)">\n+##INFO=<ID=SGB,Number=1,Type=Float,Description="Segregation based metric.">\n+##INFO=<ID=MQ0F,Number=1,Type=Float,Description="Fraction of MQ0 reads (smaller is better)">\n+##INFO=<ID=I16,Number=16,Type=Float,Description="Auxiliary tag used for calling, see description of bcf_callret1_t in bam2bcf.h">\n+##INFO=<ID=QS,Number=R,Type=Float,Description="Auxiliary tag used for calling">\n+##FORMAT=<ID=PL,Number=G,Type=Integer,Description="List of Phred-scaled genotype likelihoods">\n+##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Number of high-quality bases">\n+##INFO=<ID=DPR,Number=R,Type=Integer,Description="Number of high-quality bases observed for each allele">\n+##FORMAT=<ID=DP4,Number=4,Type=Integer,Description="Number of high-quality ref-fwd, ref-reverse, alt-fwd and alt-reverse bases">\n+##FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic depths (high-quality bases)">\n+##INFO=<ID=ADF,Number=R,Type=Integer,Description="Total allelic depths on the forward strand (high-quality bases)">\n+##INFO=<ID=ADR,Number=R,Type=Integer,Description="Total allelic depths on the reverse strand (high-quality bases)">\n+#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\tHG00100\tHG00101\tHG00102\tHG99999\n+17\t1\t.\tA\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t2\t.\tA\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t3\t.\tG\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4\t.\tC\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t5\t.\tT\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t6\t.\tT\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t7\t.\tC\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0\tPL:DP:DP4:AD\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t8\t.\tT\t<*>\t0\t.\tDP=14;ADF=0,0;ADR=0,0;DPR=0,0;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,'..b'L:DP:DP4:AD\t0,9,103:3:2,1,0,0:3,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4083\t.\tC\t<*>\t0\t.\tDP=3;ADF=2,0;ADR=1,0;DPR=3,0;I16=2,1,0,0,104,3648,0,0,134,6338,0,0,22,340,0,0;QS=1,0;MQSB=1;MQ0F=0\tPL:DP:DP4:AD\t0,9,98:3:2,1,0,0:3,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4084\t.\tA\t<*>\t0\t.\tDP=2;ADF=1,0;ADR=1,0;DPR=2,0;I16=1,1,0,0,78,3050,0,0,97,4969,0,0,20,298,0,0;QS=1,0;MQSB=1;MQ0F=0\tPL:DP:DP4:AD\t0,6,74:2:1,1,0,0:2,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4085\t.\tC\t<*>\t0\t.\tDP=2;ADF=1,0;ADR=1,0;DPR=2,0;I16=1,1,0,0,62,1940,0,0,97,4969,0,0,18,260,0,0;QS=1,0;MQSB=1;MQ0F=0\tPL:DP:DP4:AD\t0,6,62:2:1,1,0,0:2,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4086\t.\tG\t<*>\t0\t.\tDP=2;ADF=1,0;ADR=1,0;DPR=2,0;I16=1,1,0,0,56,1640,0,0,97,4969,0,0,16,226,0,0;QS=1,0;MQSB=1;MQ0F=0\tPL:DP:DP4:AD\t0,6,56:2:1,1,0,0:2,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4087\t.\tC\t<*>\t0\t.\tDP=2;ADF=1,0;ADR=1,0;DPR=2,0;I16=1,1,0,0,69,2405,0,0,97,4969,0,0,14,196,0,0;QS=1,0;MQSB=1;MQ0F=0\tPL:DP:DP4:AD\t0,6,68:2:1,1,0,0:2,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4088\t.\tA\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,39,1521,0,0,37,1369,0,0,13,169,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,37:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4089\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,36,1296,0,0,37,1369,0,0,12,144,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,36:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4090\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,33,1089,0,0,37,1369,0,0,11,121,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,33:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4091\t.\tT\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,36,1296,0,0,37,1369,0,0,10,100,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,36:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4092\t.\tG\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,37,1369,0,0,37,1369,0,0,9,81,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,37:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4093\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,35,1225,0,0,37,1369,0,0,8,64,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,35:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4094\t.\tT\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,40,1600,0,0,37,1369,0,0,7,49,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,37:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4095\t.\tA\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,35,1225,0,0,37,1369,0,0,6,36,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,35:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4096\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,32,1024,0,0,37,1369,0,0,5,25,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,32:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4097\t.\tA\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,35,1225,0,0,37,1369,0,0,4,16,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,35:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4098\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,31,961,0,0,37,1369,0,0,3,9,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,31:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4099\t.\tT\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,32,1024,0,0,37,1369,0,0,2,4,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,32:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4100\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,27,729,0,0,37,1369,0,0,1,1,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,27:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n+17\t4101\t.\tC\t<*>\t0\t.\tDP=1;ADF=1,0;ADR=0,0;DPR=1,0;I16=1,0,0,0,26,676,0,0,37,1369,0,0,0,0,0,0;QS=1,0;MQ0F=0\tPL:DP:DP4:AD\t0,3,26:1:1,0,0,0:1,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\t0,0,0:0:0,0,0,0:0,0\n'
b
diff -r d37e313a41be -r 95c6f3f4c77f test-data/summary.pdf
b
Binary file test-data/summary.pdf has changed