Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Evidence for joins in contigs #764

Open
donkirkby opened this issue Oct 28, 2021 · 0 comments
Open

Evidence for joins in contigs #764

donkirkby opened this issue Oct 28, 2021 · 0 comments
Milestone

Comments

@donkirkby
Copy link
Member

Sometimes, assembled contigs look strange, and you wonder whether they are real or a misassembly. Could we add a display to show evidence for some of the strange features in contigs? The first two features I can think of are large deletions and contig sections in a different order from the reference. I found examples of each in the 22 Oct 2021 run: Y04280A-HIVB_S10 has both. Contig 1 has a large deletion of most of gag that is covered by contig 2. Contig 1 also has a section out of order with 1.3 aligning in the 5' LTR. (I wonder if 1.3 should actually align in the 3' LTR, and then it would just be a large deletion of GP41 covered by contig 3.) Y04279A-HIVB_S9 has almost the same deletion in gag, but without the other problem.
In the remapped version, we get clues from the amount of soft clipping when the reads don't agree with the reference. Seeing that displayed across the whole contig might be helpful. I also wonder if it would be helpful to show the bias toward right or left at each point in the contig. If reads are random fragments of the genome, then we'd expect any point away from the ends to be covered by roughly the same number of reads extending farther left than right as those that extend farther right than left. As you come close to the ends of the genome or close to a misassembly, you'd expect more reads to extend into the correct assembly than would extend across the misassembly or off the end.
Maybe plotting the bias with a diverging colour map like bwr would be helpful. Blue could show left bias and red could show right bias. Alternatively, just use a line plot of right bias from 0.0 to 1.0 across each contig.

@CBeelen CBeelen added this to the 7.17 milestone May 25, 2023
@Donaim Donaim modified the milestones: 7.18, near future Jul 11, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants