Hello,
I am having issues whith ROAST in order to perform Multiz aligment of 37 genomes. I performed pairwise alignments with lastz for all the species and followed the Genome Alignments tool protocol for chains and netting. Once I have my pairwise aligments (all of them against the same reference genome), I tried to run ROAST to perform a multiz aligment. At the beginning the analysis ran smoothly for a week, generating 27 _MZ_* temp files, suggesting me thtat the tree and most of pairwised aligmnents are correct (including the headers), but suddenly stopped with an error of inconsistent row size in N line number from a tmp file.
I checked the line number and it corresponded to p.hastatus sample (see below), I checked and the lenght of the sequence match with the number of the maf format (1118)
This is how the failed batch looks like:
s raegyptiacus.scaffold_m13_p_17 25775659 1112 + 46074814 TTATGTGACAGCAGG-.....
s mblainvillei.PVJD01000009.1 642256 1115 - 1092840 CTATCTGACAGCAGATTG....
s drotundus.KZ672072.1 9922121 1118 - 10751721 CTATCTGACAGCAGATTG....
s mhirsuta.scaffold_7851 99997 1120 + 116090 CTATCTGACAGCAGAGTGGGG....
s acaudifer.scaffold_6171 177889 1124 + 194096 CTATCTGACAGCAGATTG....
s gsoricina.Super-Scaffold_40 33140221 1124 - 52801020 CTATCTGACAGCAGA....
s pdiscolor.scaffold_m19_p_15 4298889 1118 + 53102546 CTATCTGACAGCAG....
Did anyone have seen a similar error?
