262
262
done
263
263
```
264
264
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 |
271
270
272
271
If you want to restrict to the truth BED files, use this command:
273
272
@@ -287,12 +286,11 @@ done
287
286
288
287
Which resulted in this table:
289
288
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 |
296
294
297
295
### Rtg vcfstats
298
296
0 commit comments