262262done
263263```
264264
265- | Sample | [ 3] ts | [ 4] tv | [ 5] ts/tv | [ 6] ts (1st | [ 7] tv (1st | [ 8] ts/tv (1st |
266- : : : : : ALT) : ALT) : ALT) :
267- | ------ | ----- | ----- | -------- | ---------- | ---------- | ------------- |
268- | HG002 | 29955 | 11693 | 2.56 | 29942 | 11673 | 2.57 |
269- | HG003 | 29852 | 11769 | 2.54 | 29842 | 11746 | 2.54 |
270- | HG004 | 30048 | 11838 | 2.54 | 30037 | 11821 | 2.54 |
265+ | Sample | [ 3] ts | [ 4] tv | [ 5] ts/tv | [ 6] ts (1st ALT) | [ 7] tv (1st ALT) | [ 8] ts/tv (1st ALT) |
266+ | ------ | ----- | ----- | -------- | --------------- | --------------- | ------------------ |
267+ | HG002 | 29955 | 11693 | 2.56 | 29942 | 11673 | 2.57 |
268+ | HG003 | 29852 | 11769 | 2.54 | 29842 | 1174 | 2.54 |
269+ | HG004 | 30048 | 11838 | 2.54 | 3003 | 11821 | 2.54 |
271270
272271If you want to restrict to the truth BED files, use this command:
273272
@@ -287,12 +286,11 @@ done
287286
288287Which resulted in this table:
289288
290- | Sample | [ 3] ts | [ 4] tv | [ 5] ts/tv | [ 6] ts (1st | [ 7] tv (1st | [ 8] ts/tv (1st |
291- : : : : : ALT) : ALT) : ALT) :
292- | ------ | ----- | ----- | -------- | ---------- | ---------- | ------------- |
293- | HG002 | 27716 | 10549 | 2.63 | 27708 | 10536 | 2.63 |
294- | HG003 | 27382 | 10527 | 2.60 | 27378 | 10515 | 2.60 |
295- | HG004 | 27503 | 10607 | 2.59 | 27496 | 10596 | 2.59 |
289+ | Sample | [ 3] ts | [ 4] tv | [ 5] ts/tv | [ 6] ts (1st ALT) | [ 7] tv (1st ALT) | [ 8] ts/tv (1st ALT) |
290+ | ------ | ----- | ----- | -------- | --------------- | --------------- | ------------------ |
291+ | HG002 | 27716 | 10549 | 2.63 | 27708 | 10536 | 2.63 |
292+ | HG003 | 27382 | 10527 | 2.60 | 27378 | 10515 | 2.60 |
293+ | HG004 | 27503 | 10607 | 2.59 | 27496 | 10596 | 2.59 |
296294
297295### Rtg vcfstats
298296
0 commit comments