Mercurial > repos > jdv > nanopolish

changeset 9:f1141f6a2d65 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/nanopolish commit 0ace87b59137f1ed770db97fbb33036e16205edf
author jdv
date Mon, 12 Feb 2018 23:15:26 -0500
parents b437c0a7ca04
children c00a942cfc0b
files nanopolish_variants.pl nanopolish_variants.xml test-data/consensus.fa test-data/consensus.hp.fa test-data/double.bam
diffstat 5 files changed, 56 insertions(+), 10 deletions(-) [+]
line wrap: on
line diff
--- a/nanopolish_variants.pl	Mon Feb 12 03:06:42 2018 -0500
+++ b/nanopolish_variants.pl	Mon Feb 12 23:15:26 2018 -0500
@@ -118,11 +118,11 @@
 open my $out_cons, '>', $fn_consensus
     or die "Failed to open output consensus: $!";
 for (@fa_files) {
-    open my $in, '<', $_;
-    while (my $line = <$in>) {
-        print {$out_cons} $line;
+    my $parser =  BioX::Seq::Stream->new($_);
+    while (my $seq = $parser->next_seq) {
+        $seq->id =~ s/^.+\K:\d+-\d+$//; # strip coordinates from ID
+        print {$out_cons} $seq->as_fasta;
     }
-    close $in;
 }
 close $out_cons;
 
--- a/nanopolish_variants.xml	Mon Feb 12 03:06:42 2018 -0500
+++ b/nanopolish_variants.xml	Mon Feb 12 23:15:26 2018 -0500
@@ -104,8 +104,12 @@
             <param name="fix_homopolymers" value="True"/>
             <output name="out_consensus" file="consensus.hp.fa" compare="diff" />
             <output name="out_variants" file="consensus.hp.vcf" compare="diff" />
+            <assert_command>
+                <has_text text="--fix-homopolymers" />
+                <not_has_text text="--calculate-all-support" />
+            </assert_command>
         </test>
-        <!-- test pre-calculated index input -->
+        <!-- test pre-calculated index input, fix-homopolymers command -->
         <test>
             <param name="input_reads" value="called.fa" ftype="fasta" />
             <param name="input_fast5" value="test.fast5.tar.gz" ftype="fast5.tar.gz" />
@@ -117,7 +121,19 @@
             <output name="out_variants" file="consensus.hp.vcf" compare="diff" />
             <assert_command>
                 <has_text text="--index" />
-                <has_text text="--fix-homopolymers" />
+            </assert_command>
+        </test>
+        <!-- test calculate-all-support command -->
+        <test>
+            <param name="input_reads" value="called.fa" ftype="fasta" />
+            <param name="input_fast5" value="test.fast5.tar.gz" ftype="fast5.tar.gz" />
+            <param name="input_bam"   value="called.bam" ftype="bam" />
+            <param name="input_ref"   value="ref.fa" ftype="fasta" />
+            <param name="input_index" value="index.tar" ftype="tar" />
+            <param name="calculate_all_support" value="True"/>
+            <assert_command>
+                <not_has_text text="--fix-homopolymers" />
+                <has_text text="--calculate-all-support" />
             </assert_command>
         </test>
     </tests>
--- a/test-data/consensus.fa	Mon Feb 12 03:06:42 2018 -0500
+++ b/test-data/consensus.fa	Mon Feb 12 23:15:26 2018 -0500
@@ -1,2 +1,17 @@
->ref:0-921
-ACGGTACTTCGTTCGAGGACACGTATTGCTGGTGCTGCATGCCAGCGACTACGCATTCTCATTTAACCTTTCTGTTGGTGCTGATATTGCGAGGTTACCTCTACCAAGCTGGAGAAGCCACACCATCCGCCAAATACAATCCTTTCGAAATGGGTTTATCGAAGACTGTGATCACCGTTTTCCTGAATCGTACTGACAACAGACAATCTATCCTGATTGCCCCAGTTGATTCCGCTGTCTTTAAAATTAGAAAAACACGGGTAGAAGAGTCTGGATCCCGGCTAACACATTTCTGGAAGGCGCAGTATGGAAAACCAAACGCTCACCGGGATCTTAGTACCTTGACACTGATCACTCGGATTGCACTTGATACTAAAGTTGTATCTGCCTGACAAGCTCTCTGACGGCAGGTCACTTGCACACCACCAGGGATTATAAGTAACAGGATAAAGGCAATCAATATACCTCATCTCGAACAGGGTCCAATCATGAACCAAGTTGCTCCCGGAGCCATGCCCAAGGACCAAAGAGGGCATGTGCAAGAGCCCCATTAGAGGCATAACAGAACACTGACCTCACTCCAAGCCCCTTGGCGATTCCATCCCACCAGATACAAAGTCTGTGCCCTTCTGGAGGAAAGGAGGCAAACATAGGTGCTGTTATCGGTAGTGTAGCTCTTGGGGTTTGCAACAGCGGCGCAGATAACAGCAGCTGCAGCCCTAATACAAGCCGCCGAGAATGCCGCCAGCATCCTCCGGCTTAAGGAGAGCATTGCTGCAACCAATGAAGCTGTGCATGAGGTCGCAACGGATTGTCACAATAGCGGTAGGGAGAAGATCTACAGCAATTTGTTAATGAGAAGATAGAGCGACAGGCAAGTGGGTTAAATGAATGCGTAGTCGCTGTACCTTTTCAACAATACGTAACTGA
+>ref
+ACGGTACTTCGTTCGAGGACACGTATTGCTGGTGCTGCATGCCAGCGACTACGCATTCTC
+ATTTAACCTTTCTGTTGGTGCTGATATTGCGAGGTTACCTCTACCAAGCTGGAGAAGCCA
+CACCATCCGCCAAATACAATCCTTTCGAAATGGGTTTATCGAAGACTGTGATCACCGTTT
+TCCTGAATCGTACTGACAACAGACAATCTATCCTGATTGCCCCAGTTGATTCCGCTGTCT
+TTAAAATTAGAAAAACACGGGTAGAAGAGTCTGGATCCCGGCTAACACATTTCTGGAAGG
+CGCAGTATGGAAAACCAAACGCTCACCGGGATCTTAGTACCTTGACACTGATCACTCGGA
+TTGCACTTGATACTAAAGTTGTATCTGCCTGACAAGCTCTCTGACGGCAGGTCACTTGCA
+CACCACCAGGGATTATAAGTAACAGGATAAAGGCAATCAATATACCTCATCTCGAACAGG
+GTCCAATCATGAACCAAGTTGCTCCCGGAGCCATGCCCAAGGACCAAAGAGGGCATGTGC
+AAGAGCCCCATTAGAGGCATAACAGAACACTGACCTCACTCCAAGCCCCTTGGCGATTCC
+ATCCCACCAGATACAAAGTCTGTGCCCTTCTGGAGGAAAGGAGGCAAACATAGGTGCTGT
+TATCGGTAGTGTAGCTCTTGGGGTTTGCAACAGCGGCGCAGATAACAGCAGCTGCAGCCC
+TAATACAAGCCGCCGAGAATGCCGCCAGCATCCTCCGGCTTAAGGAGAGCATTGCTGCAA
+CCAATGAAGCTGTGCATGAGGTCGCAACGGATTGTCACAATAGCGGTAGGGAGAAGATCT
+ACAGCAATTTGTTAATGAGAAGATAGAGCGACAGGCAAGTGGGTTAAATGAATGCGTAGT
+CGCTGTACCTTTTCAACAATACGTAACTGA
--- a/test-data/consensus.hp.fa	Mon Feb 12 03:06:42 2018 -0500
+++ b/test-data/consensus.hp.fa	Mon Feb 12 23:15:26 2018 -0500
@@ -1,2 +1,17 @@
->ref:0-921
-ACGGTACTTCGTTCGAGGACACGTATTGCTGGTGCTGCATGCCAGCGACTACGCATTCTCATTTAACCTTTCTGTTGGTGCTGATATTGCGAGGTTACCTCTACCAAGCTGGAGAAGCCACACCATCCGCCAAATACAATCCTTTCGAAATGGGGTTTATCGAAGACTGTGATCACCGTTTTTCCTGAATCGTACTGACAACAGACAATCTATCCTGATTGCCCCAGTTGATTCCGCTGTCTTTAAAATTAGAAAAACACGGGTAGAAGAGTCTGGATCCCGGCTAACACATTTCTGGAAGGCGCAGTATGGAAACCAAACGCTCACCGGGATCTTAGTACCTTGACACTGATCACTCGGATTGCACTTGATACTAAAGTTGTATCTGCCTGACAAGCTCTCTGACGGCAGGTCACTTGCACACCACCAGGGATTATAAGTAACAGGATAAAGGCAATCAATATACCTCATCTCGAACAGGGTCCAATCATGAACCAAGTTGCTCCCGGAGCCATGCCCAAGGACCAAAGAGGGCATGTGCAAGAGCCCCATTAGAGGCATAACAGAACACTGACCTCACTCCAAGCCCCTTGGCGATTCCATCCCACCAGATACAAAGTCTGTGCCCTTCTGGAGGAAAGGAGGCAAAACATAGGTGCTGTTATCGGTAGTGTAGCTCTTGGGGTTTGCAACAGCGGCGCAGATAACAGCAGCTGCAGCCCTAATACAAGCCGCCGAGAATGCCGCCAGCATCCTCCGGCTTAAGGAGAGCATTGCTGCAACCAATGAAGCTGTGCATGAGGTCGCAACGGATTGTCACAATAGCGGTAGGGAGAAGATCTACAGCAATTTGTTAATGAGAAGATAGAGCGACAGGCAAGTGGGTTAAATGAATGCGTAGTCGCTGTACCTTTTCAACAATACGTAACTGA
+>ref
+ACGGTACTTCGTTCGAGGACACGTATTGCTGGTGCTGCATGCCAGCGACTACGCATTCTC
+ATTTAACCTTTCTGTTGGTGCTGATATTGCGAGGTTACCTCTACCAAGCTGGAGAAGCCA
+CACCATCCGCCAAATACAATCCTTTCGAAATGGGGTTTATCGAAGACTGTGATCACCGTT
+TTTCCTGAATCGTACTGACAACAGACAATCTATCCTGATTGCCCCAGTTGATTCCGCTGT
+CTTTAAAATTAGAAAAACACGGGTAGAAGAGTCTGGATCCCGGCTAACACATTTCTGGAA
+GGCGCAGTATGGAAACCAAACGCTCACCGGGATCTTAGTACCTTGACACTGATCACTCGG
+ATTGCACTTGATACTAAAGTTGTATCTGCCTGACAAGCTCTCTGACGGCAGGTCACTTGC
+ACACCACCAGGGATTATAAGTAACAGGATAAAGGCAATCAATATACCTCATCTCGAACAG
+GGTCCAATCATGAACCAAGTTGCTCCCGGAGCCATGCCCAAGGACCAAAGAGGGCATGTG
+CAAGAGCCCCATTAGAGGCATAACAGAACACTGACCTCACTCCAAGCCCCTTGGCGATTC
+CATCCCACCAGATACAAAGTCTGTGCCCTTCTGGAGGAAAGGAGGCAAAACATAGGTGCT
+GTTATCGGTAGTGTAGCTCTTGGGGTTTGCAACAGCGGCGCAGATAACAGCAGCTGCAGC
+CCTAATACAAGCCGCCGAGAATGCCGCCAGCATCCTCCGGCTTAAGGAGAGCATTGCTGC
+AACCAATGAAGCTGTGCATGAGGTCGCAACGGATTGTCACAATAGCGGTAGGGAGAAGAT
+CTACAGCAATTTGTTAATGAGAAGATAGAGCGACAGGCAAGTGGGTTAAATGAATGCGTA
+GTCGCTGTACCTTTTCAACAATACGTAACTGA
Binary file test-data/double.bam has changed