Question: Runs of homozygosity calculation
1
gravatar for BAGeno
2.2 years ago by
BAGeno170
BAGeno170 wrote:

Hi,

I am using plink to calculate runs of homozygosity.I prepared plink format file with vcftools I ran following command

plink --file Plink.input --homozyg --out Plink.result

I got empty Plink.result.hom, there are headings in this file but no results.

I repeated above mentioned procedure with 1000 genomes file and I got results. So I concluded there is problem with my file.

I have searched bcftools roh command but problem is that it needs genetic map file as input. My files are of hg38 build and I have not found genetic map of hg38 build.

Can anyone please tell me how should I calculate roh from my file.

Here are some lines from my file.

  1 11796321    rs1801133   G   A   2199.01 .   BaseCounts=22,0,19,0;BaseQRankSum=-0.851;DB;Dels=0;FS=1.257;GC=59.35;HaplotypeScore=0.9997;MQ=60;MQ0=0;MQRankSum=0.499;QD=13.85;ReadPosRankSum=1.732;DP=4533;AF=0.5;MLEAC=1;MLEAF=0.5;AN=58;AC=32   GT:AD:DP:GQ:PL:FL   0/1:20,14:34:99:295,0,588:7 ./.:.:.:.:.:.   1/1:0,46:46:99:1354,108,0:9 0/1:45,34:76:99:860,0,1218:7    ./.:.:.:.:.:.   ./.:.:.:.:.:.   1/1:0,61:60:99:1864,156,0:6 0/1:24,16:38:99:444,0,647:9 0/1:14,10:24:99:220,0,360:0 0/1:24,15:38:99:423,0,702:9 ./.:.:.:.:.:.   0/1:26,22:46:99:559,0,698:9 0/1:20,16:36:99:377,0,589:7 ./.:.:.:.:.:.   0/1:20,32:50:99:776,0,630:9 ./.:.:.:.:.:.   0/1:26,13:38:99:301,0,857:9 ./.:.:.:.:.:.   0/1:15,13:28:99:345,0,357:2 0/1:17,14:30:99:326,0,444:4 1/1:0,35:35:93:1136,93,0:8  ./.:.:.:.:.:.   0/1:12,13:26:99:340,0,392:1 ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   0/1:31,21:50:99:496,0,854:9 0/1:18,18:36:99:492,0,479:0 ./.:.:.:.:.:.   0/1:22,22:44:99:453,0,554:4 0/1:18,24:40:99:603,0,427:7 ./.:.:.:.:.:.   0/1:17,28:43:99:648,0,523:4 0/1:15,24:38:99:605,0,376:5 0/1:27,31:57:99:842,0,807:9 0/1:18,16:33:99:387,0,565:3 0/1:18,18:36:99:410,0,480:1 0/1:23,29:52:99:662,0,694:9 ./.:.:.:.:.:.   0/1:22,19:39:99:539,0,564:2 ./.:.:.:.:.:.   0/1:23,33:54:99:812,0,658:9 0/1:19,23:41:99:639,0,568:5 0/1:28,20:46:99:571,0,753:9 ./.:.:.:.:.:.   ./.:.:.:.:.:.
1   11846011    rs5065  A   G   1949.01 .   BaseCounts=19,0,33,0;BaseQRankSum=-0.227;DB;Dels=0;FS=13.341;GC=49.38;HaplotypeScore=0;MQ=60;MQ0=0;MQRankSum=-0.207;QD=19.33;ReadPosRankSum=-1.547;DP=2886;AF=0.5;MLEAC=1;MLEAF=0.5;AN=0;AC=0   GT:AD:DP:GQ:PL:FL   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.
1   16379580    rs6603859   G   T   2787.01 .   BaseCounts=0,0,19,15;BaseQRankSum=-0.173;DB;Dels=0;FS=0;GC=41.15;HaplotypeScore=0;MQ=60;MQ0=0;MQRankSum=0.277;QD=14.62;ReadPosRankSum=-1.387;DP=8264;AF=0.5;MLEAC=1;MLEAF=0.5;AN=80;AC=58   GT:AD:DP:GQ:PL:FL   0/1:15,20:34:99:573,0,500:7 0/1:14,16:29:99:508,0,474:1 1/1:0,36:36:99:1390,105,0:4 1/1:0,76:76:99:2787,214,0:5 1/1:0,30:30:84:1102,84,0:5  ./.:.:.:.:.:.   0/1:35,29:61:99:998,0,961:7 0/1:13,17:29:99:537,0,396:1 1/1:0,29:29:84:1072,84,0:3  1/1:0,44:44:99:1683,129,0:9 1/1:0,23:23:66:862,66,0:0   0/1:23,26:47:99:827,0,769:9 0/1:21,22:42:99:675,0,710:9 0/1:18,19:36:99:618,0,567:9 ./.:.:.:.:.:.   ./.:.:.:.:.:.   0/1:16,8:24:99:202,0,558:0  1/1:0,32:32:96:1289,96,0:2  0/1:18,18:36:99:613,0,551:8 0/1:18,21:38:99:678,0,596:9 1/1:0,37:37:99:1451,111,0:9 0/1:14,21:34:99:595,0,433:9 1/1:0,36:36:99:1371,102,0:9 1/1:0,30:30:81:1078,81,0:3  0/1:13,15:27:99:476,0,433:0 1/1:1,46:46:99:1766,138,0:9 1/1:0,51:51:99:2005,153,0:9 ./.:.:.:.:.:.   ./.:.:.:.:.:.   0/1:27,20:45:99:665,0,807:7 0/1:18,18:36:99:601,0,582:2 1/1:0,36:36:99:1415,108,0:0 1/1:0,40:40:99:1502,114,0:1 ./.:.:.:.:.:.   0/1:28,23:49:99:776,0,841:4 0/1:27,19:44:99:620,0,827:0 1/1:0,48:48:99:1830,141,0:9 0/1:16,19:34:99:556,0,541:1 1/1:0,44:44:99:1717,129,0:8 1/1:1,44:44:99:1728,132,0:5 0/1:35,29:61:99:917,0,1027:9    0/1:19,28:45:99:808,0,591:3 0/1:15,22:36:99:722,0,498:1 0/1:27,19:44:99:608,0,768:7 1/1:0,49:49:99:1895,144,0:6 0/1:20,11:30:99:348,0,663:9
1   97082391    rs67376798  T   A   854.01  .   BaseCounts=27,1,0,16;BaseQRankSum=-0.138;DB;Dels=0;FS=4.669;GC=35.66;HaplotypeScore=0.9789;MQ=60;MQ0=0;MQRankSum=-1.193;QD=19.41;ReadPosRankSum=-0.415;DP=44;AF=0.5;MLEAC=1;MLEAF=0.5;AN=0;AC=0 GT:AD:DP:GQ:PL:FL   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.   ./.:.:.:.:.:.
plink bcftools roh • 748 views
ADD COMMENTlink modified 2.2 years ago by Wietje210 • written 2.2 years ago by BAGeno170
1
gravatar for Wietje
2.2 years ago by
Wietje210
Germany
Wietje210 wrote:

PLINK assumes that you have a ped and map file when you use --file (same prefix, different file ending); did you try to specify for VCF format with --vcf? PLINK would otherwise also allow you to convert your VCF to ped and map with a --recode option. Does this help you along?

ADD COMMENTlink written 2.2 years ago by Wietje210

I have used this command

vcftools --vcf Input.vcf --plink --out Plink.input

It gave me both ped and map file.

ADD REPLYlink written 2.2 years ago by BAGeno170
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1648 users visited in the last hour