Skip to content

Commit a5bc0ea

Browse files
committed
Add test for breakend detection
Depends on samtools/htslib#1858
1 parent f602523 commit a5bc0ea

File tree

4 files changed

+29
-0
lines changed

4 files changed

+29
-0
lines changed

test/norm.breakend.1.1.out

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,8 @@
1+
##fileformat=VCFv4.5
2+
##FILTER=<ID=PASS,Description="All filters passed">
3+
##INFO=<ID=MATEID,Number=A,Type=String,Description="ID of mate breakend">
4+
##contig=<ID=chrA,length=500>
5+
##contig=<ID=chrB,length=500>
6+
#CHROM POS ID REF ALT QUAL FILTER INFO
7+
chrA 321 bnd_V T ]chrB:123]AGTNNNNNCAT 6 PASS MATEID=bnd_U
8+
chrB 123 bnd_U C CAGTNNNNNCA[chrA:321[ 6 PASS MATEID=bnd_V

test/norm.breakend.1.fa

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
>chrA
2+
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
3+
ACTGACTGACTGACTGACTGCCCCCACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGA
4+
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
5+
ACTGACTGACTGACTGACTGTTTTTACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGA
6+
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
7+
8+
>chrB
9+
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
10+
ACTGACTGACTGACTGACTGCCCCCACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGA
11+
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
12+
ACTGACTGACTGACTGACTGTTTTTACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGAACTGA
13+
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN

test/norm.breakend.1.vcf

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
##fileformat=VCFv4.5
2+
##INFO=<ID=MATEID,Number=A,Type=String,Description="ID of mate breakend">
3+
##contig=<ID=chrA,length=500>
4+
##contig=<ID=chrB,length=500>
5+
#CHROM POS ID REF ALT QUAL FILTER INFO
6+
chrA 321 bnd_V T ]chrB:123]AGTNNNNNCAT 6 PASS MATEID=bnd_U
7+
chrB 123 bnd_U C CAGTNNNNNCA[chrA:321[ 6 PASS MATEID=bnd_V

test/test.pl

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -271,6 +271,7 @@
271271
run_test(\&test_vcf_query,$opts,in=>'query.header',out=>'query.98.2.out',args=>q[-HH -f'%CHROM %POS[ %SAMPLE][ %DP][ %GT]']);
272272
run_test(\&test_vcf_query,$opts,in=>'query.filter-or',out=>'query.filter-or.1.out',args=>q[-f'[%SAMPLE %DP\\n]' -i'DP=1 || DP=2']);
273273
run_test(\&test_vcf_query,$opts,in=>'query.filter-or',out=>'query.filter-or.2.out',args=>q[-f'[%SAMPLE %DP\\n]' -i'DP=1 | DP=2']);
274+
run_test(\&test_vcf_norm,$opts,in=>'norm.breakend.1',fai=>'norm.breakend.1',out=>'norm.breakend.1.1.out',args=>'-m -');
274275
run_test(\&test_vcf_norm,$opts,in=>'norm.sort',out=>'norm.sort.1.out',args=>'-m -');
275276
run_test(\&test_vcf_norm,$opts,in=>'norm.sort',out=>'norm.sort.2.out',args=>'-m - -S lex');
276277
run_test(\&test_vcf_norm,$opts,in=>'norm.join-missing-ploidy',out=>'norm.join-missing-ploidy.1.out',args=>'-m +both');

0 commit comments

Comments
 (0)