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

Unequal lengths in sequence #66

Open
bioramg opened this issue Mar 11, 2022 · 3 comments
Open

Unequal lengths in sequence #66

bioramg opened this issue Mar 11, 2022 · 3 comments

Comments

@bioramg
Copy link

bioramg commented Mar 11, 2022

I am trying to polish ONT contigs with Illumina reads. But getting an error while using the Illumina pair-end reads.

bwa mem -t 40 -x ont2d ./ont_assembly.fasta /home/Ill_paired.fastq >ont_mapping.sam

racon -m 8 -x -6 -g -8 -w 500 -t 14 /home/Ill_paired.fastq ./ont_mapping.sam ./ont_assembly.fasta > ./ont_polish_racon.fasta

[racon::Overlap::transmute] error: unequal lengths in sequence and overlap file for sequence ST-E00181:968:HGM7HCCX2:1:1101:20496:2610!

I am looking forward to hearing from you.

@rvaser
Copy link
Collaborator

rvaser commented Mar 11, 2022

Hello,
your Illumina pairs probably have the same name which is not supported in the current implementation. Please use this script to modify them before aligning and polishing.

Best regards,
Robert

@bioramg
Copy link
Author

bioramg commented Mar 11, 2022

Thank you for your prompt response.

I am basic to use scripts and all. Is there any other option to modify using the basic commands?

Thank you.

@rvaser
Copy link
Collaborator

rvaser commented Mar 13, 2022

You can just run it with python racon_preprocess.py Ill_paired.fastq > Ill_paired.preproc.fastq.

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