ANGSD: Analysis of next generation Sequencing Data

Latest tar.gz version is (0.938/0.939 on github), see Change_log for changes, and download it here.

Vcf: Difference between revisions

From angsd
Jump to navigation Jump to search
No edit summary
Line 29: Line 29:
</div>
</div>
Notice that the sampleDI simply an ind followed by and integer. These relate to the samples form the -b filelist
Notice that the sampleDI simply an ind followed by and integer. These relate to the samples form the -b filelist
==NB==
This is still very beta, please let us know if it doesn't work. I'm especially concerned about the scaling of gls...
==References==
==References==

Revision as of 11:32, 23 May 2022

Newer versions of angsd (master 27april2015) supports basic vcf output. This will only include gl and gp tags which can be useful for certain external programs

Supply

-doVcf 1

Which is simply a wrapper around -gl -dopost -domajorminor -domaf.


Example

A full example commandline is given below:

 ./angsd -b list.list -dovcf 1 -gl 1 -dopost 1 -domajorminor 1 -domaf 1 -snp_pval 1e-6

Output files

This will generate a vcf file called angsdput.vcf.gz

angsdput.vcf.gz

##FORMAT=<ID=GL,Number=G,Type=Float,Description="scaled Genotype Likelihoods (these are really llh eventhough they sum to one)">
#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	ind0	ind1	ind2	ind3	ind4	ind5	ind6	ind7	ind8	ind9	ind10	ind11	ind12	ind13	ind14	ind15	ind16	ind17	ind18	ind19	ind20	ind21	ind22	ind23	ind24	ind25	ind26	ind27	ind28	ind29	ind30	ind31	ind32
1	14000202	.	G	A	.	PASS	.	GL:GP	0.013409,0.986591,0.000001:0.009959,0.990038,0.000003	0.532745,0.394297,0.072957:0.333333,0.333333,0.333333	0.532745,0.394297,0.072957:0.333333,0.333333,0.333333	0.729804,0.270070,0.000126:0.666110,0.333052,0.000839	0.532745,0.394297,0.072957:0.333333,0.333333,0.333333	0.013409,0.986589,0.000003:0.009959,0.990026,0.000015	0.843814,0.156129,0.000057:0.799685,0.199918,0.000397	0.003405,0.996582,0.000013:0.002523,0.997410,0.000068	0.915318,0.084679,0.000003:0.888870,0.111106,0.000025	0.843862,0.156138,0.000000:0.800003,0.199997,0.000000	0.001243,0.728985,0.269772:0.000420,0.333191,0.666388	0.955789,0.044211,0.000000:0.941178,0.058822,0.000000	0.843860,0.156137,0.000003:0.799987,0.199993,0.000020	0.000001,0.999999,0.000000:0.000001,0.999999,0.000000	0.824856,0.152621,0.022523:0.689947,0.172484,0.137569	0.701963,0.259767,0.038270:0.526842,0.263419,0.209740	0.532745,0.394297,0.072957:0.333333,0.333333,0.333333	0.004272,0.995611,0.000117:0.003164,0.996204,0.000632	0.532745,0.394297,0.072957:0.333333,0.333333,0.333333	0.021392,0.791754,0.186854:0.008712,0.435643,0.555645	0.843681,0.156104,0.000215:0.798816,0.199700,0.001484	0.729804,0.270070,0.000126:0.666110,0.333052,0.000839	0.002698,0.996374,0.000928:0.001990,0.993012,0.004997	0.977395,0.022605,0.000000:0.969698,0.030302,0.000000	0.532745,0.394297,0.072957:0.333333,0.333333,0.333333	0.002152,0.997847,0.000002:0.001593,0.998397,0.000010	0.915213,0.084669,0.000118:0.888149,0.111016,0.000835	0.032903,0.965645,0.001452:0.024405,0.967730,0.007865	0.000006,0.999994,0.000000:0.000005,0.999995,0.000000	0.701963,0.259767,0.038270:0.526842,0.263419,0.209740	0.843846,0.156135,0.000019:0.799896,0.199971,0.000133	0.907127,0.083921,0.008952:0.835380,0.104420,0.060201	0.332066,0.644735,0.023200:0.241926,0.634652,0.123421
1	14000873	.	G	A	.	PASS	.	GL:GP	0.000000,0.124151,0.875849:0.000000,0.030302,0.969698	0.692531,0.305335,0.002134:0.659698,0.329846,0.010456	0.900720,0.099279,0.000000:0.888891,0.111109,0.000000	0.993158,0.006842,0.000000:0.992249,0.007751,0.000000	0.999140,0.000860,0.000000:0.999024,0.000976,0.000000	0.000091,0.994381,0.005528:0.000078,0.975325,0.024596	0.900720,0.099279,0.000000:0.888891,0.111109,0.000000	0.819369,0.180627,0.000004:0.799988,0.199993,0.000018	0.000000,0.124151,0.875849:0.000000,0.030302,0.969698	0.693938,0.305955,0.000107:0.666316,0.333155,0.000529	0.986410,0.013590,0.000000:0.984616,0.015384,0.000000	0.693919,0.305947,0.000135:0.666224,0.333109,0.000666	0.007101,0.988541,0.004358:0.006172,0.974344,0.019484	0.000000,1.000000,0.000000:0.000000,1.000000,0.000000	0.030451,0.672875,0.296674:0.013127,0.328956,0.657917	0.973184,0.026816,0.000000:0.969698,0.030302,0.000000	0.973184,0.026816,0.000000:0.969698,0.030302,0.000000	0.993158,0.006842,0.000000:0.992249,0.007751,0.000000	0.004535,0.995465,0.000000:0.004001,0.995998,0.000000	0.000395,0.693734,0.305871:0.000167,0.333276,0.666557	0.900720,0.099279,0.000000:0.888891,0.111109,0.000000	0.693965,0.305967,0.000068:0.666446,0.333220,0.000334	0.947768,0.052232,0.000000:0.941178,0.058822,0.000000	0.000000,0.220881,0.779119:0.000000,0.058822,0.941178	0.900720,0.099279,0.000000:0.888891,0.111109,0.000000	0.000000,0.017410,0.982590:0.000000,0.003891,0.996109	0.000046,0.999954,0.000000:0.000040,0.999960,0.000000	0.000001,0.999999,0.000000:0.000001,0.999999,0.000000	0.947768,0.052232,0.000000:0.941178,0.058822,0.000000	0.310721,0.689279,0.000000:0.284441,0.715559,0.000000	0.900720,0.099279,0.000000:0.888891,0.111109,0.000000	0.973184,0.026816,0.000000:0.969698,0.030302,0.000000	0.000992,0.693320,0.305688:0.000420,0.333191,0.666388
1	14001018	.	T	C	.	PASS	.	GL:GP	0.000000,0.069163,0.930837:0.000000,0.015384,0.984616	0.826258,0.173742,0.000000:0.800002,0.199997,0.000002	0.159620,0.840380,0.000000:0.137752,0.862248,0.000000	0.950058,0.049942,0.000000:0.941178,0.058822,0.000000	0.904865,0.095135,0.000000:0.888891,0.111109,0.000000	0.826258,0.173742,0.000000:0.800003,0.199997,0.000000	0.495544,0.416810,0.087646:0.333333,0.333333,0.333333	0.993472,0.006528,0.000000:0.992249,0.007751,0.000000	0.001055,0.703204,0.295741:0.000420,0.333191,0.666388	0.703933,0.296042,0.000025:0.666580,0.333287,0.000133	0.826258,0.173742,0.000000:0.800003,0.199997,0.000000	0.703933,0.296042,0.000025:0.666580,0.333287,0.000133	0.000420,0.703651,0.295929:0.000167,0.333276,0.666557	0.000008,0.999992,0.000000:0.000006,0.999994,0.000000	0.085860,0.913986,0.000153:0.073175,0.926086,0.000739	0.904865,0.095135,0.000000:0.888891,0.111109,0.000000	0.495544,0.416810,0.087646:0.333333,0.333333,0.333333	0.993472,0.006528,0.000000:0.992249,0.007751,0.000000	0.703731,0.295957,0.000313:0.665554,0.332774,0.001672	0.000007,0.543141,0.456852:0.000002,0.199997,0.800001	0.904865,0.095135,0.000000:0.888891,0.111109,0.000000	0.904865,0.095135,0.000000:0.888891,0.111109,0.000000	0.904865,0.095135,0.000000:0.888891,0.111109,0.000000	0.000007,0.543141,0.456852:0.000002,0.199997,0.800001	0.904865,0.095135,0.000000:0.888892,0.111108,0.000000	0.000000,0.229117,0.770883:0.000000,0.058822,0.941178	0.045152,0.954211,0.000637:0.038160,0.958794,0.003045	0.000599,0.999394,0.000006:0.000504,0.999465,0.000031	0.974389,0.025611,0.000000:0.969698,0.030302,0.000000	0.000005,0.543142,0.456853:0.000002,0.199997,0.800002	0.974389,0.025611,0.000000:0.969698,0.030302,0.000000	0.703916,0.296035,0.000050:0.666492,0.333243,0.000265	0.495544,0.416810,0.087646:0.333333,0.333333,0.333333
1	14001867	.	A	G	.	PASS	.	GL:GP	0.000407,0.698622,0.300971:0.000167,0.333276,0.666557	0.902773,0.097227,0.000000:0.888891,0.111109,0.000000	0.000645,0.698456,0.300899:0.000265,0.333243,0.666492	0.902773,0.097227,0.000000:0.888891,0.111109,0.000000	0.902773,0.097227,0.000000:0.888891,0.111109,0.000000	0.000006,0.998276,0.001719:0.000005,0.992066,0.007929	0.986717,0.013283,0.000000:0.984616,0.015384,0.000000	0.822776,0.177224,0.000000:0.800002,0.199997,0.000001	0.000003,0.537165,0.462832:0.000001,0.199997,0.800002	0.948903,0.051097,0.000000:0.941178,0.058822,0.000000	0.948903,0.051097,0.000000:0.941178,0.058822,0.000000	0.986717,0.013283,0.000000:0.984616,0.015384,0.000000	0.000000,0.537166,0.462833:0.000000,0.199997,0.800003	0.000060,0.999940,0.000000:0.000051,0.999949,0.000000	0.000079,0.999921,0.000000:0.000068,0.999932,0.000000	0.902773,0.097227,0.000000:0.888891,0.111109,0.000000	0.902773,0.097227,0.000000:0.888891,0.111109,0.000000	0.999159,0.000841,0.000000:0.999024,0.000976,0.000000	0.004641,0.995359,0.000000:0.004001,0.995999,0.000000	0.000645,0.698456,0.300899:0.000265,0.333243,0.666492	0.948903,0.051097,0.000000:0.941178,0.058822,0.000000	0.822774,0.177223,0.000002:0.799994,0.199995,0.000011	0.996646,0.003354,0.000000:0.996109,0.003891,0.000000	0.000000,0.008985,0.991015:0.000000,0.001949,0.998051	0.822776,0.177224,0.000000:0.800002,0.199997,0.000001	0.000000,0.367208,0.632792:0.000000,0.111109,0.888891	0.000000,0.999135,0.000865:0.000000,0.995998,0.004001	0.005821,0.994177,0.000002:0.005019,0.994970,0.000011	0.993314,0.006686,0.000000:0.992249,0.007751,0.000000	0.000000,0.999998,0.000002:0.000000,0.999992,0.000008	0.948903,0.051097,0.000000:0.941178,0.058822,0.000000	0.822776,0.177224,0.000001:0.800000,0.199997,0.000003	0.007316,0.992667,0.000017:0.006309,0.993611,0.000079
[capped]

Notice that the sampleDI simply an ind followed by and integer. These relate to the samples form the -b filelist

References