You do not have permission to delete messages in this group
Copy link
Report message
Sign in to report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to gen...@soe.ucsc.edu
Hi,
I’m trying to align 2 whole genomes and am following the procedure described in this wiki page http://genomewiki.ucsc.edu/index.php/Whole_genome_alignment_howto.
Basically, I’ve done the axtChain
| chainMergeSort
| chainPreNet
| chainNet
| netSyntenic
pipeline. I would now like to filter out, i.e. exclude from the net file, overlapping alignments in the target sequence, keeping only the best one in terms of alignment score, and also to remove other alignments that for example are shorter than 30% for that
aligned query.
I looked at `netFilter` and `netChainSubset` but don’t think they would do what I want. I could try to write some parser myself but was wondering if this is something that is already available and that I have missed.
Thanks,
Pedro
Christopher Lee
unread,
Apr 18, 2018, 7:58:37 PM4/18/18
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Sign in to report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to Coelho De Almeida, Pedro, gen...@soe.ucsc.edu
Hello Pedro,
Thank you for your question about whole genome alignments and
filtering your net file. Can you provide an example of these
overlapping alignments on the target sequence that you would like to
filter? Nets are by definition single coverage on the target, although
lower level chains may be used to "fill in" gaps in the top level
chain.