gravatar for wes

3 hours ago by

I had downloaded RSII data from NCBI which consists of three *.bax.h5.1 files and one bas.h5.1 file.

Next, I would like to obtain subreads file using pbh5tools.

By referring to guideline below, the input is bas.h5 file but there is error. May I know if the cmd below is correct? or should I use baxh5.1 as input file which I tried as well (shown below) but error too.

lira.no-ip.org:8080/doc/python-pbh5tools/html/#installation

bash5tools.py --outFilePrefix m150803_212153_42216_c100858342550000001823192601241650_s1_p0.1.subreads --readType subreads --outType fastq m150803_212153_42216_c100858342550000001823192601241650_s1_p0.bas.h5.1

Traceback (most recent call last):
File "/home/cbr01/anaconda3/envs/py2/bin/bash5tools.py", line 166, in
sys.exit(BasH5ToolsRunner().start())
File "/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/util/ToolRunner.py", line 85, in start
return self.run()
File "/home/cbr01/anaconda3/envs/py2/bin/bash5tools.py", line 121, in run
inBasH5 = BasH5Reader(self.args.inFile)
File "/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py", line 660, in init
for fn in self.file["/MultiPart/Parts"] ]
File "/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py", line 304, in init
raise IOError("Invalid or nonexistent bax/bas file %s" % filename)
IOError: Invalid or nonexistent bax/bas file /media/cbr01/analysis/WEE/DNASeq/SRX2718652/m150803_212153_42216_c100858342550000001823192601241650_s1_p0.1.bax.h5

bash5tools.py --outFilePrefix m150803_212153_42216_c100858342550000001823192601241650_s1_p0.1.subreads --readType subreads --outType fastq m150803_212153_42216_c100858342550000001823192601241650_s1_p0.1.bax.h5.1
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:273: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
numEvent = h5Group["ZMW/NumEvent"].value
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:274: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
holeNumber = h5Group["ZMW/HoleNumber"].value
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:349: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
holeNumbers = self._mainBasecallsGroup["ZMW/HoleNumber"].value
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:356: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
fh["/PulseData/Regions"].value)
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:379: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
holeStatus = self._mainBasecallsGroup["ZMW/HoleStatus"].value
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:388: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
(self._mainBasecallsGroup["ZMW/NumEvent"].value > 0) &
/home/cbr01/anaconda3/envs/py2/lib/python2.7/site-packages/pbcore/io/BasH5IO.py:565: H5pyDeprecationWarning: dataset.value has been deprecated. Use dataset[()] instead.
self.__metricCache[name] = self._mainBasecallsGroup[k].value

link

written
3 hours ago
by

wes20



Source link