MetaQC is calculating uniquely mapping ratio, NRF, PBC...
Samtools quickcheck detects that /home/leelee/biodata/shTRA2A/markdup/shTRA2A-IP.markdup.bam is complete!
Samtools quickcheck detects that /home/leelee/biodata/shTRA2A/markdup/shTRA2A-input.markdup.bam is complete!
Time taken: 1.2238 min
MetaQC is calculating rRNA & tRNA contamination ratio...
Time taken: 1.0155 min
MetaQC is calculating saturation...
Time taken: 2.9656 min
Metaqc is calling peaks...
Peak calling step may take a few hours...
Process Process-2:
Traceback (most recent call last):
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/metaqc/modules/peakcalling/callpeakNew.py", line 564, in callPeak
readCounts(gtfile=gtffile, sample=item, path=path)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/metaqc/modules/peakcalling/callpeakNew.py", line 203, in readCounts
single_tx_new = single_tx_new.astype(data_type)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/generic.py", line 5790, in astype
col.astype(dtype=dtype[col_name], copy=copy, errors=errors)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/generic.py", line 5806, in astype
new_data = self._mgr.astype(dtype=dtype, copy=copy, errors=errors)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/internals/managers.py", line 414, in astype
return self.apply("astype", dtype=dtype, copy=copy, errors=errors)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/internals/managers.py", line 327, in apply
applied = getattr(b, f)(**kwargs)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/internals/blocks.py", line 592, in astype
new_values = astype_array_safe(values, dtype, copy=copy, errors=errors)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/dtypes/cast.py", line 1300, in astype_array_safe
new_values = astype_array(values, dtype, copy=copy)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/dtypes/cast.py", line 1248, in astype_array
values = astype_nansafe(values, dtype, copy=copy)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/pandas/core/dtypes/cast.py", line 1165, in astype_nansafe
return lib.astype_intsafe(arr, dtype)
File "pandas/_libs/lib.pyx", line 679, in pandas._libs.lib.astype_intsafe
ValueError: invalid literal for int() with base 10: '0.0'
shTRA2A-IP.markdup read length: 50
shTRA2A-input.markdup read length: 50
Time taken: 4.9873 min
MetaQC is calculating read distribution dispersion...
Process Process-3:
Traceback (most recent call last):
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/multiprocessing/process.py", line 315, in _bootstrap
self.run()
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/multiprocessing/process.py", line 108, in run
self._target(*self._args, **self._kwargs)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/metaqc/modules/stdev/sd.py", line 102, in sd
IP_readCounts_stdev, INPUT_readCounts_stdev = stdev_value(sample=item[0], path=path)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/metaqc/modules/stdev/sd.py", line 44, in stdev_value
readCounts_INPUT_dict[txID] = [int(readCountsINPUT)]
ValueError: invalid literal for int() with base 10: '0.0'
Time taken: 0.0002 min
MetaQC is calculating FRiP value...
sort: cannot read: metaqc_test/tmp/shTRA2A-IP.markdup.treated.peak.bed: No such file or directory
Time taken: 0.6047 min
MetaQC is plotting metagene...
<metagene> metaqc_test/peak/shTRA2A-IP.markdup_peak.xls No exists!
Time taken: 0.0002 min
MetaQC is calculating peak length distribution...
<peakLengthDis> metaqc_test/peak/shTRA2A-IP.markdup_peak.xlsNo exists!
Time taken: 0.0002 min
MetaQC is calculating the fraction of ovlerlapping-peak length between replicates...
No Replicate Samples!
Time taken: 0.0001 min
Traceback (most recent call last):
File "/home/leelee/miniconda3/envs/metaqc/bin/metaqc", line 476, in <module>
main()
File "/home/leelee/miniconda3/envs/metaqc/bin/metaqc", line 468, in main
makeHtml(path=dirname, htmlTemplate= html_File, replicateNumber=replicate_number)
File "/home/leelee/miniconda3/envs/metaqc/lib/python3.9/site-packages/metaqc/modules/html/displayResults.py", line 173, in makeHtml
fh_RDDinfo = open(path + '/file/stdevInfo.txt')
FileNotFoundError: [Errno 2] No such file or directory: 'metaqc_test/file/stdevInfo.txt'
I got the ideal result when running with the file in the test folder. I guess it might be a gtf file, but the URL provided is "ftp://igenome:G3nom3s4u@ussd-ftp.illumina.com/Homo_sapiens/UCSC/hg19/Homo_sapiens_UCSC_hg19.tar.gz" cannot be opened. This troubles me very much, I don’t know if you have any suggestions to solve it.
In addition, how should I cite the metaqc tool? I didn't find any documents in github that can be cited.
thanks,
LeeLee
Dear sir or madam,
Thanks for providing such an excellent tool, but I ran into some troubles during the operation, the following is my operation log:
I got the ideal result when running with the file in the test folder. I guess it might be a gtf file, but the URL provided is "ftp://igenome:G3nom3s4u@ussd-ftp.illumina.com/Homo_sapiens/UCSC/hg19/Homo_sapiens_UCSC_hg19.tar.gz" cannot be opened. This troubles me very much, I don’t know if you have any suggestions to solve it.
In addition, how should I cite the metaqc tool? I didn't find any documents in github that can be cited.
thanks,
LeeLee