Skip to main content
added 190 characters in body
Source Link
gringer
  • 15.1k
  • 5
  • 24
  • 83

I am trying to use dorado basecaller to process ONT data. I need to make the move table and I am doing it with --emit-move.

dorado basecaller [email protected]/ /pod5_files/ --recursive --emit-moves --modified-bases 5mCG_5hmCG -v > /bams/mod_stride1.bam

The main purpose of the move table for me is to calculate the number of ONT device measurements before basecaller emits one base (calculate translocation time per base). But the stride value in config.toml file is fixed to 5 (stride=5). I need to calculate the instantaneous translocation time per base (in seconds), because I am looking to find a measure similar to IPD values in PacBio, when we can get the time for sequencing each base.

I was wondering to change the stride value in config.toml file in the model folder, from stride=5 to stride=1, to directly infer the number of measurements without using Tombo for correcting and re-squiggling.

Unfortunately, even though this command works perfectly fine when the stride=5, when I set stride=1 it always gives me the following error: "CUDA out of memory", even when I triple the size of RAM for the job.

My questions are:

  1. Does setting stride=1 make sense for this purpose? Or do I have to do something else? (I know we have some K-mer models too)
  2. Is this even possible to change the stride value on config.toml file?
  3. How should I figure out how much RAM I need for the job?

I am trying to use dorado basecaller to process ONT data. I need to make the move table and I am doing it with --emit-move.

dorado basecaller [email protected]/ /pod5_files/ --recursive --emit-moves --modified-bases 5mCG_5hmCG -v > /bams/mod_stride1.bam

The main purpose of the move table for me is to calculate the number of ONT device measurements before basecaller emits one base (calculate translocation time per base). But the stride value in config.toml file is fixed to 5 (stride=5)

I was wondering to change the stride value in config.toml file in the model folder, from stride=5 to stride=1, to directly infer the number of measurements without using Tombo for correcting and re-squiggling.

Unfortunately, even though this command works perfectly fine when the stride=5, when I set stride=1 it always gives me the following error: "CUDA out of memory", even when I triple the size of RAM for the job.

My questions are:

  1. Does setting stride=1 make sense for this purpose? Or do I have to do something else? (I know we have some K-mer models too)
  2. Is this even possible to change the stride value on config.toml file?
  3. How should I figure out how much RAM I need for the job?

I am trying to use dorado basecaller to process ONT data. I need to make the move table and I am doing it with --emit-move.

dorado basecaller [email protected]/ /pod5_files/ --recursive --emit-moves --modified-bases 5mCG_5hmCG -v > /bams/mod_stride1.bam

The main purpose of the move table for me is to calculate the number of ONT device measurements before basecaller emits one base (calculate translocation time per base). But the stride value in config.toml file is fixed to 5 (stride=5). I need to calculate the instantaneous translocation time per base (in seconds), because I am looking to find a measure similar to IPD values in PacBio, when we can get the time for sequencing each base.

I was wondering to change the stride value in config.toml file in the model folder, from stride=5 to stride=1, to directly infer the number of measurements without using Tombo for correcting and re-squiggling.

Unfortunately, even though this command works perfectly fine when the stride=5, when I set stride=1 it always gives me the following error: "CUDA out of memory", even when I triple the size of RAM for the job.

My questions are:

  1. Does setting stride=1 make sense for this purpose? Or do I have to do something else? (I know we have some K-mer models too)
  2. Is this even possible to change the stride value on config.toml file?
  3. How should I figure out how much RAM I need for the job?
Source Link
Marjan
  • 309
  • 1
  • 6

How should I use dorado basecaller to calculate translocation time? Can I change models' config files?

I am trying to use dorado basecaller to process ONT data. I need to make the move table and I am doing it with --emit-move.

dorado basecaller [email protected]/ /pod5_files/ --recursive --emit-moves --modified-bases 5mCG_5hmCG -v > /bams/mod_stride1.bam

The main purpose of the move table for me is to calculate the number of ONT device measurements before basecaller emits one base (calculate translocation time per base). But the stride value in config.toml file is fixed to 5 (stride=5)

I was wondering to change the stride value in config.toml file in the model folder, from stride=5 to stride=1, to directly infer the number of measurements without using Tombo for correcting and re-squiggling.

Unfortunately, even though this command works perfectly fine when the stride=5, when I set stride=1 it always gives me the following error: "CUDA out of memory", even when I triple the size of RAM for the job.

My questions are:

  1. Does setting stride=1 make sense for this purpose? Or do I have to do something else? (I know we have some K-mer models too)
  2. Is this even possible to change the stride value on config.toml file?
  3. How should I figure out how much RAM I need for the job?