vdjtools
vdjtools copied to clipboard
A bug of OverlapPair: div1 and div2 are wrong when overlap type isn't nt
Hi, I find a bug of OverlapPair: div1 and div2 are wrong when 'overlap type' isn't 'nt'.
Here is an example: java -jar vdjtools.jar OverlapPair -i aa sample1.txt sample2.txt aa/result java -jar vdjtools.jar OverlapPair -i nt sample1.txt sample2.txt nt/result java -jar vdjtools.jar OverlapPair -i strict sample1.txt sample2.txt strict/result
$ cat nt/result.paired.nt.summary.txt | cut -f1-6 1_sample_id 2_sample_id div1 div2 div12 div21 sample1 sample2 62917 76574 4236 4236
$ cat aa/result.paired.aa.summary.txt | cut -f1-6 1_sample_id 2_sample_id div1 div2 div12 div21 sample1 sample2 62917 76574 5392 5392
$ cat strict/result.paired.strict.summary.txt | cut -f1-6 1_sample_id 2_sample_id div1 div2 div12 div21 sample1 sample2 62917 76574 4096 4096