mccortex icon indicating copy to clipboard operation
mccortex copied to clipboard

Breakpoint caller should use novel kmers

Open noporpoise opened this issue 9 years ago • 0 comments

The breakpoint caller should mark kmers if they are used in a call. Unused kmers that are in a sample and not the referece ("novel kmers") should be used to seed a second breakpoint caller.

The novel kmer breakpoint caller should perform a breadth-first search (BFS) to find the nearest kmer that occurs in the reference either side. The shortest path between these kmers should be reported as a putative breakpoint call. This caller will have a high false positive rate, but may catch a few reasonable SNPs / indels.

noporpoise avatar Aug 10 '15 17:08 noporpoise