Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Improper #225

Open
BirdmanRidesAgain opened this issue Dec 23, 2022 · 1 comment
Open

Improper #225

BirdmanRidesAgain opened this issue Dec 23, 2022 · 1 comment

Comments

@BirdmanRidesAgain
Copy link

Hi,
I've been attempting to run racon on a set of draft assemblies. Racon has worked on other assemblies in the past, but when I attempt to run the following line of code:

racon -t 60 ../../../ONT-data/ArdeotisONT/ardeotisONT_trimreads.fq.gz ardeotisONT_flye-assembly_remapped-SO.sam ../polished-flye/ardeotisONT_flye-assembly.fasta > ardeotisONT_assembly-racon.fasta

...this is the output, along with an empty file:

[racon::Polisher::initialize] loaded target sequences 7.326854 s


terminate called after throwing an instance of 'std::invalid_argument'
  what():  [bioparser::SamParser] error: invalid file format
Aborted (core dumped)

I'm not certain which file is the problematic one, but I used the following code to generate the overlaps file, as I think that's the most likely to be problematic.

minimap2 -t 60 ../polished-flye/ardeotisONT_flye-assembly.fasta -ax map-ont ../../../ONT-data/ArdeotisONT/ardeotisONT_trimreads.fq.gz | samtools sort > houbaraONT_flye-assembly_remapped-SO.sam

Any insight you have into what might be wrong with these commands would be greatly appreciated.

@rvaser
Copy link
Collaborator

rvaser commented Jan 17, 2023

Which racon/minimap2 version are you using? Did you by any chance forgot option -a while running minimap2? Invalid file format usually indicates mismatch between the file extension and the file format.

Best regards,
Robert

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants