DaehwanKimLab / hisat-genotype

GNU General Public License v3.0
23 stars 15 forks source link

Error when running tutorials code with --assembly #77

Open Runda-Xu opened 1 year ago

Runda-Xu commented 1 year ago

When running "hisatgenotype --base hla --locus-list A -1 ILMN/NA12892.extracted.1.fq.gz -2 ILMN/NA12892.extracted.2.fq.gz", everything is alright.

However, when I added --assembly, the program terminated after outputting the result of DMA. The log file shows this error:

Base and Files:hla NA12892.extracted.1.fq.gz-hla-extracted-1.fq.gz NA12892.extracted.1.fq.gz-hla-extracted-2.fq.gz multiprocessing.pool.RemoteTraceback: """ Traceback (most recent call last): File "/home/.conda/envs/neo/lib/python3.10/multiprocessing/pool.py", line 125, in worker result = (True, func(*args, **kwds)) File "/home/Software/hisatgenotype/hisatgenotype_modules/hisatgenotype_typing_core.py", line 2655, in genotyping_locus typing(simulation, File "/home/Software/hisatgenotype/hisatgenotype_modules/hisatgenotype_typing_core.py", line 1532, in typing assembly_graph.Node(node_read_id, File "/home/Software/hisatgenotype/hisatgenotype_modules/hisatgenotype_assembly_graph.py", line 197, in init assert len(seq) == len(var) AssertionError """

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

Traceback (most recent call last): File "/home/Software/hisatgenotype/./hisatgenotype", line 675, in typing_process ofnlog.write(str(x.get())) File "/home/.conda/envs/neo/lib/python3.10/multiprocessing/pool.py", line 774, in get raise self._value AssertionError

Besides, in the output fasta file, the allele names for HLA-B are "B15:01:01:04" and "Novel". While in the report, they are "B15:01:01:04", "B56:01:01:02", and "B56:01:01:03"