Mercurial > repos > mheinzl > fsd_regions
diff fsd_regions.py @ 8:6c2608e8d094 draft
planemo upload for repository https://github.com/monikaheinzl/duplexanalysis_galaxy/tree/master/tools/fsd_regions commit 8833d1a8a49d7b6d4a9c849b0335d3260564b351-dirty
| author | mheinzl |
|---|---|
| date | Tue, 20 Nov 2018 09:51:47 -0500 |
| parents | 3b8a0e462021 |
| children | eabfdc012d7b |
line wrap: on
line diff
--- a/fsd_regions.py Tue Nov 20 09:45:06 2018 -0500 +++ b/fsd_regions.py Tue Nov 20 09:51:47 2018 -0500 @@ -71,7 +71,7 @@ seqDic_ab = dict(zip(all_ab, quant_ab)) seqDic_ba = dict(zip(all_ba, quant_ba)) - if re.search(r'_(\d)+_(\d)+$', str(mut_array[0,0])) is None: + if re.search(r'(\d)+_(\d)+$', str(mut_array[0,0])) is None: seq_mut, seqMut_index = numpy.unique(numpy.array(mut_array[:, 1]), return_index=True) group = mut_array[seqMut_index,0] mut_array = mut_array[seqMut_index,:] @@ -156,7 +156,7 @@ for i, count in zip(groupUnique, quantAfterRegion): index_of_current_region = numpy.where(group == i)[0] plt.text(0.55, 0.14 - s, "{}=\n".format(i), size=11, transform=plt.gcf().transFigure) - if re.search(r'_(\d)+_(\d)+$', str(mut_array[0, 0])) is None: + if re.search(r'(\d)+_(\d)+$', str(mut_array[0, 0])) is None: nr_tags_ab = len(numpy.unique(mut_array[index_of_current_region, 1])) else: nr_tags_ab = len(mut_array[index_of_current_region, 1])