|
|
@@ -415,6 +415,7 @@ impl Dorado {
|
|
|
if !bam_path.exists() {
|
|
|
info!("Creating new bam file");
|
|
|
self.basecall_align(&dorado_bin)?;
|
|
|
+ self.index()?;
|
|
|
} else {
|
|
|
// check if merge before call
|
|
|
let new_bam_path = bam_path
|
|
|
@@ -430,7 +431,6 @@ impl Dorado {
|
|
|
self.merge_bam(&new_bam_path)?;
|
|
|
}
|
|
|
|
|
|
- self.index()?;
|
|
|
self.run_cramino()?;
|
|
|
self.run_modkit()?;
|
|
|
// self.create_fastq()?;
|