SubPhaser
SubPhaser copied to clipboard
IndexError: index -1 is out of bounds for axis 0 with size 0
Hi, Thanks for developing the tool. I tried the example of ginger and successfully procressed. But When I used my own triploid genome (3n=63), I met an error. My config file is as follow: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 21 22 23 24 25 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 45 47 48 49 50 51 52 53 54 59 60 64 65 67 68 71 72 73 74 75 76 77
The command was 'subphaser -i ref.fa -c config.txt -pre out', The I get the error like this:
22-06-02 16:24:49 [INFO] Summary of overall LTR insertion age (million years):
/home/wangyue/software/miniconda2/envs/SubPhaser/lib/python3.8/site-packages/numpy/core/fromnumeric.py:3440: RuntimeWarning: Mean of empty slice.
return _methods._mean(a, axis=axis, dtype=dtype,
/home/wangyue/software/miniconda2/envs/SubPhaser/lib/python3.8/site-packages/numpy/core/_methods.py:189: RuntimeWarning: invalid value encountered in double_scalars
ret = ret.dtype.type(ret / rcount)
Traceback (most recent call last):
File "/home/wangyue/software/miniconda2/envs/SubPhaser/bin/subphaser", line 33, in
And the results in the file "outk15_q200_f2.chrom-subgenome.tsv" showed different number of chromosomes for each genotype.
I don't know where is my problem. Can you give me any advises? Thanks a lot
Can you show me the whole log and the plots produced by subphaser? I suspect that the genome is not well phased.
As clearly shown by the results, the subgenomes are not well assigned, because there are too few differential kmers between subgenomes. It may be artifactual or natural. My suggestions (refer to the Supplementary Material):
- Check the HiC contact map whether there are too many switch errors between subgenomes. If not, it is natural.
- If natural, it may be because 1) there are too many recombinations between subgenomes after hybridization, or 2) their progenitors are too close (possible autopolyploid) to distinguish. In these cases, subgenome phasing is meaningless and subphaser will not work.
So many chromosomes may also be an issue. You can have a try to reduce the number of homeologous chromosome sets by setting the config like:
1 2 3
...
13 14 15
16
17
18
...
75
76
77
Remove check point files by rm tmp/*ok
before re-runing to avoid confusing.
Thank you for your quick reply. I will take a try