schneebergerlab / syri

Synteny and Rearrangement Identifier
https://schneebergerlab.github.io/syri/
MIT License
303 stars 36 forks source link

i was using syri with two genome files its give me a error #208

Closed Hemaprasanth24 closed 7 months ago

Hemaprasanth24 commented 10 months ago

(py3.8) hemaprasanth@hemaprasanth-Aspire-A715-51G:~/main_syri/radiodurans_protealyticus$ syri -c r_pr.sam -r radiodurans.fasta -q protealyticus.fasta -k -F S --prefix r_pr Reading Coords - WARNING - Chromosomes IDs do not match. Reading Coords - WARNING - Matching them automatically. For each reference genome, most similar query genome will be selected. Check mapids.txt for mapping used. multiprocessing.pool.RemoteTraceback: """ Traceback (most recent call last): File "/home/hemaprasanth/anaconda/envs/py3.8/lib/python3.8/multiprocessing/pool.py", line 125, in worker result = (True, func(*args, *kwds)) File "/home/hemaprasanth/anaconda/envs/py3.8/lib/python3.8/multiprocessing/pool.py", line 48, in mapstar return list(map(args)) File "syri/pyxFiles/synsearchFunctions.pyx", line 803, in syri.synsearchFunctions.syri IndexError: list index out of range """

The above exception was the direct cause of the following exception:

Traceback (most recent call last): File "/home/hemaprasanth/anaconda/envs/py3.8/bin/syri", line 6, in main(sys.argv[1:]) File "/home/hemaprasanth/anaconda/envs/py3.8/lib/python3.8/site-packages/syri/scripts/syri.py", line 326, in main syri(args) File "/home/hemaprasanth/anaconda/envs/py3.8/lib/python3.8/site-packages/syri/scripts/syri.py", line 214, in syri startSyri(args, coords[["aStart", "aEnd", "bStart", "bEnd", "aLen", "bLen", "iden", "aDir", "bDir", "aChr", "bChr"]]) File "syri/pyxFiles/synsearchFunctions.pyx", line 505, in syri.synsearchFunctions.startSyri File "syri/pyxFiles/synsearchFunctions.pyx", line 506, in syri.synsearchFunctions.startSyri File "/home/hemaprasanth/anaconda/envs/py3.8/lib/python3.8/multiprocessing/pool.py", line 364, in map return self._map_async(func, iterable, mapstar, chunksize).get() File "/home/hemaprasanth/anaconda/envs/py3.8/lib/python3.8/multiprocessing/pool.py", line 771, in get raise self._value IndexError: list index out of range

mnshgl0110 commented 10 months ago

Duplicate https://github.com/schneebergerlab/syri/issues/154

Hemaprasanth24 commented 10 months ago

Hii, gene_presence_absence_phylogeny

as you mentioned above i just did the reverse complement of my query genome i got the out fie but i don't now is it correct because i have i phylogenetic tree based on gene presence and absence that shows the reference r_p_c and query genome are not closely related but the synteny plot from plotsr with this out file shows syntenic regions

mnshgl0110 commented 10 months ago

Please start a separate issue for new questions.