mikolmogorov / Ragout

Chromosome-level scaffolding using multiple references
Other
146 stars 27 forks source link

ERROR: permutation ids #71

Closed ardy20 closed 3 years ago

ardy20 commented 3 years ago

Hi

I am trying to run Ragout using a .maf alignment file created with SibeliaZ. But I get the following error:

ERROR: An error occured while running Ragout: ERROR: permutation ids in /ragout/results/maf-workdir/100/blocks_coords.txt do not follow naming convention: 'genome.chromosome'

Could you please explain which file should I change and how? Is there any script for that?

Regards Ardy

sajjadasaf commented 3 years ago

I have faced the same issue by using maf alignment created with SibeliaZ.

mikolmogorov commented 3 years ago

Hi,

For the maf workflow, you need to have your fasta sequence names in a form genome.chromosome, where genome is the same name you are using in the recipe file, and chromosome is any unique identifier.

Mikhail

ardy20 commented 3 years ago

Thanks Mikhaeil My fasta file has 250 contigs. Is there any tool or script that can do this? Regards Ardy

On Tue, 15 Dec 2020 at 3:41 am, Mikhail Kolmogorov notifications@github.com wrote:

Hi,

For the maf workflow, you need to have your fasta sequence names in a form genome.chromosome, where genome is the same name you are using in the recipe file, and chromosome is any unique identifier.

Mikhail

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/fenderglass/Ragout/issues/71#issuecomment-744599131, or unsubscribe https://github.com/notifications/unsubscribe-auth/APA35LY4RZSOED3A4EV2HE3SUZE3FANCNFSM4UV6NWUA .

--

Regards

Ardy Kharabian-Masouleh (PhD)

Email: kharabian1@gmail.com

a.kharabianmasouleh@uq.edu.au

Mobile: +61405134478

mikolmogorov commented 3 years ago

I'm not aware of any scripts, unfortunately.

sajjadasaf commented 3 years ago

Dear Mikhaeil, You can do it by using awk with following command awk '/^>/{print ">genome.chromosome" ++i; next}{print}' < file.fasta > output.fasta