PacificBiosciences / pbbioconda

PacBio Secondary Analysis Tools on Bioconda. Contains list of PacBio packages available via conda.
BSD 3-Clause Clear License
251 stars 45 forks source link

CCS --hifi-kinetics throws an error #726

Closed looxon93 closed 1 week ago

looxon93 commented 1 week ago

Hi I am running ccs tool with --hifi-kinetics to prepare my BAM file for PacBio Jasmine. I am using ccs 6.4.0 (commit v6.4.0) on Linux.

I tried debug mode but the error I am getting isn't informative at all.

Command line: ccs -j 1 --hifi-kinetics --log-level DEBUG 9358-1.m84165_240624_120322_s1.hifi_reads.bam ./9358-1.m84165_240624_120322_s1.hifi_reads.ccs.bam

Error is shown below:

| 20241105 13:51:26.954 | INFO | Effective commandline     : ccs -j 1 --hifi-kinetics --log-level DEBUG 9358-1.m84165_240624_120322_s1.hifi_reads.bam ./9358-1.m84165_240624_120322_s1.hifi_reads.ccs.bam
| 20241105 13:51:26.954 | DEBUG | Polish is ON
| 20241105 13:51:26.954 | INFO | CCS Version               : 6.4.0 (commit v6.4.0)
| 20241105 13:51:26.963 | FATAL | ccs ERROR: map::at

This is conda environment I am using (micromamba actually, but that shouldn't be the problem):

  Name           Version  Build        Channel
────────────────────────────────────────────────────
  _libgcc_mutex  0.1      conda_forge  conda-forge
  _openmp_mutex  4.5      2_gnu        conda-forge
  libgcc         14.2.0   h77fa898_1   conda-forge
  libgomp        14.2.0   h77fa898_1   conda-forge
  libstdcxx      14.2.0   hc0a3c3a_1   conda-forge
  pbccs          6.4.0    h9ee0642_0   bioconda
  pbjasmine      2.4.0    h4ac6f70_0   bioconda

These are read tags available for one example read:

RG:Z:8585369d ac:B:i,6,0,8,0 ec:f:7.33613 ma:i:0 np:i:6 rq:f:0.999417 sn:B:f,12.1924,17.9366,4.5802,7.98824 we:i:5717341 ws:i:36860 zm:i:158667473 qs:i:7 qe:i:19244 MM:Z:C+m?,50,27,52,23,14,3,17,0,6,10,11,134,12,31,35,13,39,15,55,12,62,94,36,22,64,9,15,12,30,21,19,21,35,40,22,0,19,4,8,0,12,19,6,45,8,2,15,88,39,71,116,0,6,4,8,23,7,12,5,41,26,50,34,1,36,38,29,14,42,6,35,13,60,214,7,4,17,19,28,3,28,77,2,7,12,27,49,4,26,4,76,42,46,5,144,21,29,84,30,9,3,14,17,0,51,58,6,17,112,55,3,14,2,13; ML:B:C,215,231,195,135,170,255,80,155,249,118,246,136,194,4,58,214,152,95,221,61,173,11,156,131,255,104,254,248,240,240,206,139,200,190,187,251,251,210,255,180,163,185,252,118,196,239,32,115,222,149,214,127,185,247,116,255,206,109,23,120,253,193,185,224,104,249,169,248,138,87,164,179,251,204,149,94,219,235,111,112,192,170,36,160,222,120,235,73,250,191,160,103,226,188,104,224,213,161,51,36,192,232,227,203,96,169,205,98,197,162,207,196,162,197

Can you help, how to solve this?

Thanks, Luka.

armintoepfer commented 1 week ago

Your input is already HiFi.