1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
import rtest
datadir = "../samples"
vcf = datadir+"/sample.vcf"
def test01(n):
"""vcf2tsv converts a CSV to a tabulated test file, e.g.
vcf2tsv samples/sample.vcf
outputs
```
#CHROM POS ID REF ALT QUAL FILTER AA AC AF AN DB DP H2NS
19 111 . A C 9.6 . . . . . . . . .
19 112 . A G 10 . . . . . . . . .
20 14370 rs6054257 G A 29 PASS . . 0.5 . . 14.3
```
>>> rtest.run_stdout(f"vcf2tsv {vcf}", ext="tsv")
"""
pass
if __name__ == "__main__":
import doctest
rtest.setup()
doctest.testmod()
|
