# HG changeset patch # User mheinzl # Date 1612428935 0 # Node ID ced1a529e7cd251270ed4379c829df03e2ed39da # Parent ded0dc6a20d34203aab89e459049918e7f1965c7 planemo upload for repository https://github.com/Single-Molecule-Genetics/VariantAnalyzerGalaxy/tree/master/tools/variant_analyzer commit ee4a8e6cf290e6c8a4d55f9cd2839d60ab3b11c8 diff -r ded0dc6a20d3 -r ced1a529e7cd mut2read.py --- a/mut2read.py Mon Jan 25 13:21:55 2021 +0000 +++ b/mut2read.py Thu Feb 04 08:55:35 2021 +0000 @@ -74,14 +74,16 @@ stop_pos = variant.start #chrom_stop_pos = str(chrom) + "#" + str(stop_pos) ref = variant.REF - alt = variant.ALT[0] + if len(variant.ALT) > 0: + alt = variant.ALT[0] + else: + alt = variant.ALT chrom_stop_pos = str(chrom) + "#" + str(stop_pos) + "#" + ref + "#" + alt dcs_len = [] if len(ref) == len(alt): for pileupcolumn in bam.pileup(chrom, stop_pos - 1, stop_pos + 1, max_depth=100000000): - if pileupcolumn.reference_pos == stop_pos: count_alt = 0 count_ref = 0 diff -r ded0dc6a20d3 -r ced1a529e7cd mut2sscs.py --- a/mut2sscs.py Mon Jan 25 13:21:55 2021 +0000 +++ b/mut2sscs.py Thu Feb 04 08:55:35 2021 +0000 @@ -68,7 +68,10 @@ stop_pos = variant.start #chrom_stop_pos = str(chrom) + "#" + str(stop_pos) ref = variant.REF - alt = variant.ALT[0] + if len(variant.ALT) > 0: + alt = variant.ALT[0] + else: + alt = variant.ALT chrom_stop_pos = str(chrom) + "#" + str(stop_pos) + "#" + ref + "#" + alt if len(ref) == len(alt): diff -r ded0dc6a20d3 -r ced1a529e7cd read2mut.py --- a/read2mut.py Mon Jan 25 13:21:55 2021 +0000 +++ b/read2mut.py Thu Feb 04 08:55:35 2021 +0000 @@ -132,7 +132,10 @@ stop_pos = variant.start #chrom_stop_pos = str(chrom) + "#" + str(stop_pos) ref = variant.REF - alt = variant.ALT[0] + if len(variant.ALT) > 0: + alt = variant.ALT[0] + else: + alt = variant.ALT chrom_stop_pos = str(chrom) + "#" + str(stop_pos) + "#" + ref + "#" + alt if len(ref) == len(alt):