PacificBiosciences / kineticsTools

Tools for detecting DNA modifications from single molecule, real-time sequencing data
19 stars 21 forks source link

help with reprocessMotifSites.py? #4

Closed llwestphal closed 10 years ago

llwestphal commented 10 years ago

Hi, I'm trying to use reprocessingMotifSites.py and I've run into this error message. Any help is appreciated! Thanks.

usc-secure-wireless-034-135:Scripts FinkelLabAir$ python reprocessMotifSites.py --reference path/to/rev_compPACBIO.fasta --motifs /path/to/-home-smrtanalysis-smrtanalysis-current-common-jobs-016-016796-data-motifs.gff.gz --motif_summary /path/to/-home-smrtanalysis-smrtanalysis-current-common-jobs-016-016796-data-motif_summary.csv --undetected --gff ~/Desktop/new.gff ~/Downloads/-home-smrtanalysis-smrtanalysis-current-common-jobs-016-016796-data-aligned_reads.cmp.h5 Traceback (most recent call last): File "reprocessMotifSites.py", line 552, in kt.start() File "/Library/Python/2.7/site-packages/pbcore-0.8.5-py2.7.egg/pbcore/util/ToolRunner.py", line 105, in start return self.run() File "reprocessMotifSites.py", line 280, in run ret = self._mainLoop() File "reprocessMotifSites.py", line 490, in _mainLoop self.loadReferenceAndModel(self.args.reference, self.args.infile) File "reprocessMotifSites.py", line 431, in loadReferenceAndModel contigs = ReferenceUtils.loadReferenceContigs(referencePath, cmpH5Path) File "/Library/Python/2.7/site-packages/kineticsTools-0.4.0-py2.7-macosx-10.9-intel.egg/kineticsTools/ReferenceUtils.py", line 60, in loadReferenceContigs contigDict[x.MD5].id = x.ID KeyError: 'ab1a06dc3709a6f554551cc5df814102'

dalexander commented 10 years ago

@obanerjee, please help here

obanerjee commented 10 years ago

Please write back to me at my email address if the R scripts I sent do not help.