Skip to content
Snippets Groups Projects
Commit a0ca3a4e authored by Monika Komorowska's avatar Monika Komorowska
Browse files

don't limit seq_region to chromosomes, so the script can be used to fix...

don't limit seq_region to chromosomes, so the script can be used to fix overlaps in mappings between other coord_systems
parent f48b9801
No related branches found
No related tags found
No related merge requests found
#!/usr/local/ensembl/bin/perl #!/usr/bin/env perl
=head1 NAME =head1 NAME
fix_overlaps.pl - remove overlapping mappings between two closely related fix_overlaps.pl - remove overlapping mappings in the assembly coordinates (asm_start - asm_end) between two closely related assemblies
assemblies
=head1 SYNOPSIS =head1 SYNOPSIS
...@@ -19,11 +18,10 @@ Required arguments: ...@@ -19,11 +18,10 @@ Required arguments:
--assembly=ASSEMBLY assembly version ASSEMBLY --assembly=ASSEMBLY assembly version ASSEMBLY
--altassembly=ASSEMBLY alternative assembly version ASSEMBLY --altassembly=ASSEMBLY alternative assembly version ASSEMBLY
--chromosomes, --chr=LIST only process LIST seq_regions
Optional arguments: Optional arguments:
--chromosomes, --chr=LIST only process LIST toplevel seq_regions
--conffile, --conf=FILE read parameters from FILE --conffile, --conf=FILE read parameters from FILE
(default: conf/Conversion.ini) (default: conf/Conversion.ini)
...@@ -38,10 +36,9 @@ Optional arguments: ...@@ -38,10 +36,9 @@ Optional arguments:
=head1 DESCRIPTION =head1 DESCRIPTION
This script removes overlapping mappings that were generated by the code in This script removes overlapping mappings from the assembly coordinates.
align_nonident_regions.pl. Mappings are trimmed so that no overlaps are present Mappings are trimmed so that such overlaps don't break the AssemblyMapper when
in the assembly table, because such overlaps may break the AssemblyMapper when projecting between the reference and non-reference assemblies.
projecting between the two assemblies.
It also merges adjacent assembly segments which can result from alternating It also merges adjacent assembly segments which can result from alternating
alignments from clone identity and blastz alignment. alignments from clone identity and blastz alignment.
...@@ -114,7 +111,8 @@ $support->init_log; ...@@ -114,7 +111,8 @@ $support->init_log;
$support->check_required_params( $support->check_required_params(
'assembly', 'assembly',
'altassembly' 'altassembly',
'chromosomes',
); );
# database connection # database connection
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment