umi_tools whitelist error

[mohonardost]

Hello everyone

I used the following command to make a whitelist using umi-tools (v=0.5.0) and python (v=2.7.13):

umi_tools whitelist --stdin hgmm_100_R1.fastq.gz \
--bc-pattern=CCCCCCCCCCCCCCCCNNNNNNNNNN \
--set-cell-number=100 \
--plot-prefix=100_cells_whitelist \
--log2stderr > whitelist.txt

It starts to work as follow, however it bump into an error:

# output generated by whitelist --stdin hgmm_100_R1.fastq.gz --bc-pattern=CCCCCCCCCCCCCCCCNNNNNNNNNN --set-cell-number=100 --plot-prefix=100_cells_whitelist --log2stderr
# job started at Wed Sep 6 22:28:00 2017 on n126 -- a096c517-be3e-4d56-9f0e-c4753dbae35e
# pid: 80367, system: Linux 2.6.32-431.el6.x86_64 #1 SMP Fri Nov 22 03:15:09 UTC 2013 x86_64
# blacklist_tsv : None
# cell_number : 100
# compresslevel : 6
# error_correct_threshold : 1
# expect_cells : False
# extract_method : string
# filter_cell_barcodes : False
# log2stderr : True
# loglevel : 1
# method : reads
# pattern : CCCCCCCCCCCCCCCCNNNNNNNNNN
# pattern2 : None
# plot_prefix : 100_cells_whitelist
# prime3 : None
# random_seed : None
# read2_in : None
# short_help : None
# stderr : <open file '<stderr>', mode 'w' at 0x2aaaaaf0d1e0>
# stdin : <gzip open file 'hgmm_100_R1.fastq.gz', mode 'rb' at 0x2aab4f798db0 0x2aaab49f3b90>
# stdlog : <open file '<stderr>', mode 'w' at 0x2aaaaaf0d1e0>
# stdout : <open file '<stdout>', mode 'w' at 0x2aaaaaf0d150>
# subset_reads : 100000000
# timeit_file : None
# timeit_header : None
# timeit_name : all
# whitelist_tsv : None
2017-09-06 22:28:00,838 INFO Starting barcode extraction
2017-09-06 22:28:00,882 INFO Parsed 0 reads
2017-09-06 22:28:04,640 INFO Parsed 100000 reads
2017-09-06 22:28:07,552 INFO Parsed 200000 reads
2017-09-06 22:28:10,455 INFO Parsed 300000 reads
2017-09-06 22:28:13,344 INFO Parsed 400000 reads
2017-09-06 22:28:16,238 INFO Parsed 500000 reads
2017-09-06 22:28:19,126 INFO Parsed 600000 reads
2017-09-06 22:28:22,027 INFO Parsed 700000 reads
2017-09-06 22:28:24,921 INFO Parsed 800000 reads
2017-09-06 22:28:27,818 INFO Parsed 900000 reads
2017-09-06 22:28:30,720 INFO Parsed 1000000 reads
2017-09-06 22:28:33,617 INFO Parsed 1100000 reads
2017-09-06 22:28:36,510 INFO Parsed 1200000 reads
2017-09-06 22:28:39,408 INFO Parsed 1300000 reads
2017-09-06 22:28:42,302 INFO Parsed 1400000 reads
2017-09-06 22:28:45,195 INFO Parsed 1500000 reads
2017-09-06 22:28:48,087 INFO Parsed 1600000 reads
2017-09-06 22:28:50,980 INFO Parsed 1700000 reads
2017-09-06 22:28:53,870 INFO Parsed 1800000 reads
2017-09-06 22:28:56,771 INFO Parsed 1900000 reads
2017-09-06 22:28:59,670 INFO Parsed 2000000 reads
2017-09-06 22:29:02,566 INFO Parsed 2100000 reads
2017-09-06 22:29:05,458 INFO Parsed 2200000 reads
2017-09-06 22:29:08,354 INFO Parsed 2300000 reads
2017-09-06 22:29:11,248 INFO Parsed 2400000 reads
2017-09-06 22:29:14,140 INFO Parsed 2500000 reads
2017-09-06 22:29:17,031 INFO Parsed 2600000 reads
2017-09-06 22:29:19,919 INFO Parsed 2700000 reads
2017-09-06 22:29:22,819 INFO Parsed 2800000 reads
2017-09-06 22:29:25,720 INFO Parsed 2900000 reads
2017-09-06 22:29:28,614 INFO Parsed 3000000 reads
2017-09-06 22:29:31,503 INFO Parsed 3100000 reads
2017-09-06 22:29:34,395 INFO Parsed 3200000 reads
2017-09-06 22:29:37,287 INFO Parsed 3300000 reads
2017-09-06 22:29:40,185 INFO Parsed 3400000 reads
2017-09-06 22:29:43,078 INFO Parsed 3500000 reads
2017-09-06 22:29:45,970 INFO Parsed 3600000 reads
2017-09-06 22:29:48,876 INFO Parsed 3700000 reads
2017-09-06 22:29:51,772 INFO Parsed 3800000 reads
2017-09-06 22:29:54,667 INFO Parsed 3900000 reads
2017-09-06 22:29:57,561 INFO Parsed 4000000 reads
2017-09-06 22:30:00,458 INFO Parsed 4100000 reads
2017-09-06 22:30:03,348 INFO Parsed 4200000 reads
2017-09-06 22:30:06,241 INFO Parsed 4300000 reads
2017-09-06 22:30:09,135 INFO Parsed 4400000 reads
2017-09-06 22:30:12,028 INFO Parsed 4500000 reads
2017-09-06 22:30:14,929 INFO Parsed 4600000 reads
2017-09-06 22:30:17,826 INFO Parsed 4700000 reads
2017-09-06 22:30:20,728 INFO Parsed 4800000 reads
2017-09-06 22:30:23,625 INFO Parsed 4900000 reads
2017-09-06 22:30:26,524 INFO Parsed 5000000 reads
2017-09-06 22:30:29,420 INFO Parsed 5100000 reads
2017-09-06 22:30:32,320 INFO Parsed 5200000 reads
2017-09-06 22:30:35,221 INFO Parsed 5300000 reads
2017-09-06 22:30:38,132 INFO Parsed 5400000 reads
2017-09-06 22:30:41,055 INFO Parsed 5500000 reads
2017-09-06 22:30:43,952 INFO Parsed 5600000 reads
2017-09-06 22:30:46,861 INFO Parsed 5700000 reads
2017-09-06 22:30:49,757 INFO Parsed 5800000 reads
2017-09-06 22:30:52,654 INFO Parsed 5900000 reads
2017-09-06 22:30:55,553 INFO Parsed 6000000 reads
2017-09-06 22:30:58,451 INFO Parsed 6100000 reads
2017-09-06 22:31:01,347 INFO Parsed 6200000 reads
2017-09-06 22:31:04,246 INFO Parsed 6300000 reads
2017-09-06 22:31:07,151 INFO Parsed 6400000 reads
2017-09-06 22:31:10,047 INFO Parsed 6500000 reads
2017-09-06 22:31:12,946 INFO Parsed 6600000 reads
2017-09-06 22:31:15,837 INFO Parsed 6700000 reads
2017-09-06 22:31:18,732 INFO Parsed 6800000 reads
2017-09-06 22:31:21,622 INFO Parsed 6900000 reads
2017-09-06 22:31:24,517 INFO Parsed 7000000 reads
2017-09-06 22:31:27,415 INFO Parsed 7100000 reads
2017-09-06 22:31:30,242 INFO Starting - whitelist determination
Traceback (most recent call last):
File "/home/honardoostma/.conda/envs/pysam-bzlinkerr/bin/umi_tools", line 11, in <module>
sys.exit(main())
File "/home/honardoostma/.conda/envs/pysam-bzlinkerr/lib/python2.7/site-packages/umi_tools/umi_tools.py", line 50, in main
module.main(sys.argv)
File "/home/honardoostma/.conda/envs/pysam-bzlinkerr/lib/python2.7/site-packages/umi_tools/whitelist.py", line 371, in main
options.plot_prefix)
File "/home/honardoostma/.conda/envs/pysam-bzlinkerr/lib/python2.7/site-packages/umi_tools/umi_methods.py", line 399, in getCellWhitelist
cell_barcode_counts, expect_cells, cell_number, plotfile_prefix)
File "/home/honardoostma/.conda/envs/pysam-bzlinkerr/lib/python2.7/site-packages/umi_tools/umi_methods.py", line 264, in getKneeEstimate
if local_min is not None:
UnboundLocalError: local variable 'local_min' referenced before assignment

does anyone know what is the problem and how I can solve it?


Thanks