Investigate collapse_ancestry1.py
Closed this issue · 2 comments
The paint_karyograms.py script is grabbing the incorrect values in the bed files. This needs to be investigated. Bed files are skipping sections of the chromosome and may be leading to mispainting.
This bug is rescued by a bug in the karyogram painting script but still impacts the global lai script
A bug was found in both collapse_ancestry1.py and plot_karyogram1.py. The plot_karyogram1.py bug rescues the collapse ancestry. Collapse_ancestry1.py has a fence post problem where the last window in a file is not being added for the ancestry in the bed file. This creates gaps on the right side of the ancestry windows in the karyograms. Plot ancestry rescues this by grabbing the incorrect field, bp, and plotting the rectangles from a cM position to the bp position instead of cM to cM. This fills in the gap as the anc rectangle extend off the plotted areas.
This has not been fixed as it is not urgent. Closing and opening an issue to fix plotting in icebox.