zhengxwen / SeqArray

Data management of large-scale whole-genome sequence variant calls (Development version only)
http://www.bioconductor.org/packages/SeqArray
43 stars 12 forks source link

seqGetData genotype with wrong genotype on multiallelic sites #78

Closed no7ren closed 2 years ago

no7ren commented 2 years ago

I am using the latest SeqArray package with the following command: seqGetData(f, 'genotype')

And part of the output look like: sample allele [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] NA NA 0 NA 0 0 0 61 0 0 0 0 0 [2,] NA NA 0 NA 0 0 0 61 0 0 0 0 0 sample allele [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 61 0 0 0 0 0 0 20 20 0 NA NA [2,] 61 0 0 0 0 0 0 20 2 0 NA NA

while the 8th and 14th samples in the gds should contain [13, 13] instead of [61, 61]. And the 20th and 21th samples should be [4, 4] and [4, 18].

I have use the PySeqArray to confirm that the gds file looks fine where I can get the correct genotype from the same gds file. Because of the sensitivity of the data, I can only post the line of vcf file which has the problem.

chr01 10761 . CATATAT CACACACACACACACATATATATAT,CACACACACACACATATATAT,CACACACACACATATATATATATATATAT,CACATATATATAT,CACATATATATATAT,CAT,CATATATAT,CATATATATAT,CATATATATATAT,CATATATATATATAT,CATATATATATATATAT,CATATATATATATATATAT,CATATATATATATATATATAT,CATATATATATATATATATATAT,CATATATATATATATATATATATAT,CATATATATATATATATATATATATAT,CATATATATATATATATATATATATATAT,TATATAT 21797.9 . . GT:DP:FT:GQ:AD ./.:2:PASS:3:2,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:5:DRAGENIndelHardQUAL:12:3,.,.,.,.,.,.,1,.,.,.,.,.,.,.,.,.,.,. 0/0:11:LowGQ:23:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:3:PASS:25:0,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,.,. 0/0:16:LowGQ:23:15,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.0/0:31:PASS:54:30,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:21:PASS:50:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 13/13:11:PASS:29:0,.,.,.,.,.,.,.,.,.,.,.,.,11,.,.,.,.,. 0/0:26:PASS:53:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:13:DRAGENIndelHardQUAL:26:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:16:LowGQ:20:15,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|0:18:DRAGENIndelHardQUAL:22:17,.,.,.,.,.,.,.,.,.,.,.,.,.,1,.,.,.,. 0/0:17:PASS:21:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 13/13:12:PASS:36:4,.,.,.,.,.,.,.,.,.,.,.,.,8,.,.,.,.,. 0/0:25:PASS:42:24,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:30:PASS:48:29,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:13:PASS:27:13,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:12:PASS:27:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:13:LowGQ:6:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:16:PASS:31:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/4:8:PASS:27:0,.,.,.,8,.,.,.,0,.,.,.,.,.,.,.,.,.,. 4/18:18:PASS:49:1,.,.,.,6,.,.,.,0,.,.,.,.,.,.,.,.,.,11 0/0:18:PASS:39:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:14:LowGQ:0:6,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:17:LowGQ:0:7,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:30:LowGQ:0:0,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:33:LowGQ:43:31,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:21:PASS:35:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:11:DRAGENIndelHardQUAL:16:10,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:17:PASS:26:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:21:LowGQ:0:2,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:22:LowGQ:0:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:9:PASS:18:8,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/15:13:PASS:34:7,.,.,.,.,.,.,.,.,.,.,.,.,.,.,3,.,.,. 0/0:8:PASS:22:8,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/15:18:PASS:34:13,.,.,.,.,.,.,.,.,.,.,.,.,.,1,4,.,.,. 0|16:15:PASS:34:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,4,.,. 0|16:13:PASS:11:10,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,3,.,. 5/9:7:PASS:4:1,.,.,.,.,1,.,.,.,5,.,.,.,.,.,.,.,.,. 0|16:16:PASS:16:13,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,3,.,. 9/9:10:PASS:27:0,.,.,.,.,.,.,.,.,10,.,.,.,.,.,.,.,.,. 0/0:63:LowGQ:77:58,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:33:LowGQ:34:31,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:27:LowGQ:39:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:27:PASS:49:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18/18:9:PASS:25:0,.,.,.,.,.,.,.,.,.,0,.,.,.,.,.,.,.,9 0|15:21:PASS:33:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,4,.,.,. 7/7:16:PASS:27:4,.,.,.,.,.,.,12,.,.,.,.,.,.,.,.,.,.,. 7/7:13:PASS:9:0,.,.,.,.,.,.,13,.,.,.,.,.,.,.,.,.,.,0 0|18:28:PASS:24:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,7 0/7:22:PASS:45:7,.,.,.,.,.,.,14,.,.,.,.,.,.,.,.,.,.,1 0|18:24:PASS:41:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,7 0/7:13:PASS:47:8,.,.,.,.,.,.,5,.,.,.,.,.,.,.,.,.,.,0 0/0:20:LowGQ:6:7,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/9:16:PASS:35:0,.,.,.,.,.,.,.,1,15,.,.,.,.,.,.,.,.,. 13/16:9:PASS:26:0,.,.,.,.,.,.,.,.,.,.,.,2,5,.,.,2,.,. 0/0:23:PASS:59:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/9:16:PASS:47:0,.,.,.,.,.,.,.,.,15,.,.,.,.,.,.,.,.,1 18/18:8:PASS:34:0,.,.,.,.,.,.,.,.,.,.,.,.,0,.,.,.,.,8 0/18:25:PASS:48:11,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,.,14 18/18:18:PASS:30:3,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,.,15 0|16:59:PASS:43:36,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,23,.,. 0/0:34:PASS:69:34,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:27:DRAGENIndelHardQUAL:5:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 16/18:29:PASS:31:5,.,.,.,.,.,.,0,0,.,.,.,.,.,.,.,6,.,18 6/18:54:PASS:50:0,.,.,.,.,.,29,.,.,.,.,.,.,.,.,.,.,.,25 0/18:36:PASS:48:18,.,.,.,.,.,.,.,.,.,.,.,1,.,.,.,.,.,17 0/0:11:LowGQ:4:10,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|15:26:PASS:47:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,10,.,.,. ./.:5:PASS:13:0,.,.,.,.,.,.,.,.,.,.,.,0,0,.,.,.,.,5 0/0:11:PASS:26:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:6:PASS:15:0,.,.,.,.,.,.,.,.,6,.,.,.,.,.,.,.,.,. 0/0:32:PASS:62:31,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/16:38:PASS:33:28,.,.,.,.,.,.,2,.,.,.,.,.,.,.,.,8,.,. 0/0:22:PASS:53:22,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/18:45:PASS:48:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,25 0/18:8:PASS:8:2,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,6 0/16:28:PASS:31:10,.,.,.,.,.,.,.,.,.,.,.,.,9,.,1,8,.,. ./.:0:.:.:.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:2:DRAGENSnpHardQUAL:20:1,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,1 0|18:63:PASS:47:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,37 0|18:31:PASS:45:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,11 0|18:82:PASS:35:52,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,30 0/16:28:PASS:30:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,3,7,.,. 0|16:21:PASS:42:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,5,.,. 0/0:15:PASS:41:15,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/9:11:PASS:17:0,.,.,.,.,.,.,.,.,9,.,.,.,.,.,.,.,.,2 9/9:14:PASS:23:1,.,.,.,.,.,.,.,.,13,.,.,.,.,.,.,.,.,. 0/0:21:LowGQ:35:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 1/2:9:PASS:44:0,3,5,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:4:LowGQ:0:2,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/9:7:PASS:17:0,.,.,.,.,.,.,.,.,7,.,.,.,.,.,.,.,.,. ./.:3:PASS:3:0,.,.,.,.,.,.,.,1,2,.,.,.,.,.,.,.,.,. ./.:4:PASS:9:0,.,.,.,.,.,.,.,.,4,.,.,.,.,.,.,.,.,. ./.:4:PASS:9:0,.,.,.,.,.,.,.,.,4,.,.,.,.,.,.,.,.,. ./.:3:PASS:7:0,.,.,.,.,.,.,.,.,3,.,.,.,.,.,.,.,.,. 0/0:7:LowGQ:10:7,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:6:PASS:15:0,.,.,.,.,.,.,.,.,0,3,.,.,.,.,.,.,.,3 ./.:1:PASS:3:0,.,.,.,.,.,.,.,.,.,.,.,.,1,.,.,.,.,. 7/17:28:PASS:50:0,.,.,.,.,.,.,17,.,.,.,.,.,.,.,.,.,11,. 0/0:26:PASS:62:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/4:17:PASS:83:0,.,.,.,17,.,.,.,0,.,.,.,.,.,.,.,.,.,. 13/13:11:PASS:30:0,.,.,.,.,.,.,.,.,.,.,.,.,11,.,.,.,.,. 0/0:28:PASS:68:27,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|15:41:PASS:48:22,.,.,.,.,.,.,.,.,.,.,.,.,.,.,19,.,.,. 0/0:25:LowGQ:51:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 16/18:13:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,0,.,.,.,5,.,8 0|18:37:PASS:45:14,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,23 11/18:22:PASS:44:0,.,.,.,.,.,.,.,.,0,.,4,.,.,.,.,.,.,18 0/0:18:PASS:48:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:11:PASS:30:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:20:LowGQ:39:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/18:31:PASS:47:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,13 3/3:10:PASS:24:0,.,.,9,.,.,.,.,.,.,.,.,0,.,.,.,.,.,. 9/9:16:PASS:22:2,.,.,.,1,1,.,.,0,12,.,.,.,.,.,.,.,.,. 7/18:10:PASS:4:0,.,.,.,.,.,.,6,.,.,.,.,.,1,.,.,.,.,3 7/18:11:PASS:45:0,.,.,.,.,.,.,6,.,.,.,.,.,0,.,.,.,.,5 18/18:12:PASS:38:0,.,.,.,.,.,.,0,0,.,.,.,.,.,.,.,.,.,12 6/11:20:PASS:50:0,.,.,.,.,.,12,.,.,.,.,8,.,.,.,.,.,.,. 0/0:17:LowGQ:36:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/9:16:PASS:45:0,.,.,.,.,.,.,.,.,16,.,.,.,.,.,.,.,.,. 7/18:9:PASS:15:0,.,.,.,.,.,.,5,.,.,.,.,.,0,.,.,.,.,4 18/18:22:PASS:63:0,.,.,.,.,.,.,.,.,.,0,.,.,.,.,.,.,.,22 0/0:21:LowGQ:8:6,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:17:LowGQ:5:8,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/16:27:PASS:39:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,8,.,. 13/13:11:PASS:29:1,.,.,.,.,.,.,.,.,.,.,.,.,10,.,.,.,.,. 13/16:16:PASS:18:0,.,.,.,.,.,.,.,.,.,.,.,.,8,.,.,4,4,. 13/16:18:PASS:47:0,.,.,.,.,.,.,.,.,.,.,.,.,6,.,.,12,.,. 0/0:18:PASS:28:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 16/18:33:PASS:41:0,.,.,.,.,.,.,.,.,.,.,.,0,.,.,.,10,.,23 16/18:32:PASS:32:0,.,.,.,.,.,.,.,.,.,.,.,1,.,.,.,7,.,24 ./.:24:LowGQ:0:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:25:LowGQ:48:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:20:PASS:38:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 13/16:21:PASS:48:0,.,.,.,.,.,.,.,.,.,.,.,1,7,.,2,11,.,. 0/0:29:PASS:53:28,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.0/0:18:PASS:37:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:1:LowDepth;LowGQ:3:1,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:2:PASS:6:2,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:25:PASS:69:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:20:PASS:42:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.0/0:19:PASS:48:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:36:LowGQ:51:32,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:22:PASS:49:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:41:LowGQ:63:37,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18/18:17:PASS:48:0,.,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,17 6/18:24:PASS:47:1,.,.,.,.,.,14,.,.,.,.,.,.,.,.,.,.,.,9 0|18:16:DRAGENSnpHardQUAL:7:8,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,8 11/18:25:PASS:29:0,.,.,.,.,.,.,.,.,.,.,5,.,.,.,.,0,.,20 12/18:35:PASS:41:0,.,.,.,.,.,.,.,.,0,.,.,8,0,.,.,.,.,27 0|15:31:PASS:41:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,10,.,.,. 0/15:22:PASS:40:15,.,.,.,.,.,.,.,.,.,.,.,.,.,.,5,.,.,. 16/18:27:PASS:38:3,.,.,.,.,.,.,1,.,.,.,.,.,.,.,.,5,.,18 0/0:18:PASS:39:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:32:PASS:81:32,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:22:PASS:60:22,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|16:31:PASS:47:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,13,.,. 0|15:27:PASS:47:14,.,.,.,.,.,.,.,.,.,.,.,.,.,.,13,.,.,. 7/7:17:PASS:48:0,.,.,.,.,.,.,17,.,.,.,.,.,.,.,.,.,.,. 0/15:30:PASS:30:21,.,.,.,.,.,.,2,.,.,.,.,.,.,.,7,.,.,. 0/0:33:PASS:81:33,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:31:PASS:65:31,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:24:PASS:60:24,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|16:32:PASS:45:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,12,.,. 0|16:33:PASS:47:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,14,.,. 0|15:24:PASS:48:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,12,.,.,. 0/0:12:PASS:32:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/15:27:PASS:45:17,.,.,.,.,.,.,.,.,.,.,.,.,.,2,8,.,.,. 0/15:33:PASS:48:17,.,.,.,.,.,.,.,.,.,.,.,.,.,3,13,.,.,. 0/16:26:PASS:41:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,5,.,. 0|16:36:PASS:33:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,10,.,. 0/16:25:PASS:46:9,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,14,.,2 0/0:19:PASS:43:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:16:PASS:33:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|15:36:PASS:46:24,.,.,.,.,.,.,.,.,.,.,.,.,.,.,12,.,.,. 0|15:22:PASS:33:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,4,.,.,. 0/0:21:PASS:48:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:17:PASS:40:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/16:35:PASS:37:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,10,.,. 0/0:25:PASS:60:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:24:PASS:54:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/16:35:PASS:48:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,15,.,. 12/16:18:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,8,.,.,.,10,.,. 0|16:26:PASS:48:14,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,12,.,. 0|16:21:PASS:48:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,10,.,. 13/16:13:PASS:40:2,.,.,.,.,.,.,.,.,.,.,.,.,3,.,.,8,.,. 12/12:15:PASS:41:0,.,.,.,.,.,.,.,.,.,.,.,15,.,.,.,.,.,. 12/12:11:PASS:18:0,.,.,.,.,.,.,.,.,.,.,1,10,.,.,.,.,.,. 0|16:22:PASS:46:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,4,.,. 0/0:25:LowGQ:8:8,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 12/12:10:PASS:27:0,.,.,.,.,.,.,.,.,.,.,.,10,.,.,.,.,.,. 12/16:21:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,11,.,.,.,10,.,. 0/0:22:LowGQ:50:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:19:LowGQ:36:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|9:44:PASS:48:23,.,.,.,.,.,.,.,.,21,.,.,.,.,.,.,.,.,. 18/18:21:PASS:57:0,.,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,21 ./.:4:PASS:11:0,.,.,.,.,.,.,.,.,.,.,.,0,.,.,.,.,.,4 0/0:17:LowGQ:6:5,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18/18:16:PASS:53:0,.,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,16 13/15:26:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,.,13,.,13,.,.,. 0|15:35:PASS:47:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,15,.,.,. 18/18:8:PASS:27:0,.,.,.,.,.,.,.,.,.,.,.,1,.,.,.,.,.,7 0/0:25:PASS:60:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 7/7:16:PASS:44:0,.,.,.,.,.,.,16,.,.,.,.,.,.,.,.,.,.,. 0|16:36:PASS:46:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,15,.,. 0/0:18:PASS:48:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 13/16:19:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,.,10,.,.,9,.,. 18/18:15:PASS:40:0,.,.,.,.,.,.,2,.,.,.,.,.,.,.,.,.,.,13 0/0:19:PASS:42:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:19:PASS:44:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|15:22:PASS:43:15,.,.,.,.,.,.,.,.,.,.,.,.,.,.,7,.,.,. 0/16:25:PASS:32:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,2,4,.,. 13/15:16:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,.,8,.,8,.,.,. 16/18:28:PASS:49:0,.,.,.,.,.,.,.,.,.,.,.,0,2,.,.,13,.,13 0|16:28:PASS:31:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,7,.,. 12/16:15:PASS:37:1,.,.,.,.,.,.,.,.,.,.,.,11,.,.,.,3,.,. 13/16:19:PASS:49:2,.,.,.,.,.,.,.,.,.,.,.,.,11,.,.,6,.,. 4/4:11:PASS:41:1,.,.,.,10,.,.,.,0,.,.,.,.,.,.,.,.,.,. 4/16:31:PASS:50:0,.,.,.,16,.,.,.,0,.,.,.,.,.,.,.,14,.,. 4/4:8:PASS:44:0,.,.,.,8,.,.,.,0,.,.,.,.,.,.,.,.,.,. 4/9:27:PASS:49:0,.,.,.,11,0,.,.,0,16,.,.,.,.,.,.,.,.,. 0/0:23:PASS:51:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:26:PASS:49:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:20:PASS:30:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 13/16:13:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,.,7,.,.,6,.,.0/0:22:PASS:59:22,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:17:LowGQ:36:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:19:LowGQ:41:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:26:LowGQ:54:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|15:24:PASS:40:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,7,.,.,.0/0:22:LowGQ:29:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:13:LowGQ:28:13,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:20:PASS:41:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/4:13:PASS:49:0,.,.,.,13,.,.,.,0,.,.,.,.,.,.,.,.,.,. 0/14:30:PASS:43:18,.,.,.,.,.,.,.,.,.,.,.,.,2,10,.,.,.,. 0/0:21:PASS:57:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/18:31:PASS:42:8,.,.,.,.,.,.,.,.,0,.,.,.,.,.,6,.,.,17 0|15:25:PASS:45:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,9,.,.,. 0/0:16:PASS:37:16,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:19:PASS:41:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 17/18:19:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,0,.,.,.,.,5,14 0/0:13:PASS:30:13,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 16/18:47:PASS:45:2,.,.,.,.,.,.,.,.,.,.,.,.,0,.,.,13,.,32 0/0:21:PASS:42:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 8/16:25:PASS:49:2,.,.,.,.,.,.,.,12,.,.,.,.,.,.,.,11,.,. 4/4:8:PASS:25:0,.,.,.,7,.,.,.,1,.,.,.,.,.,.,.,.,.,. 0/0:13:LowGQ:24:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18|18:13:PASS:33:0,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,13 10/18:28:PASS:49:0,.,.,.,.,.,.,.,.,0,11,.,.,.,.,.,.,.,17 10/18:33:PASS:36:2,.,.,.,.,.,.,.,.,.,21,4,.,.,.,.,.,.,6 0/0:18:PASS:40:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:24:PASS:64:24,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:21:PASS:46:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:25:PASS:48:24,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:22:PASS:51:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|16:32:PASS:44:22,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,10,.,. 0/0:22:PASS:46:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:23:PASS:44:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:25:PASS:60:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18/18:12:PASS:36:0,.,.,.,.,.,.,.,0,0,.,.,.,.,.,.,.,.,12 0/0:27:LowGQ:42:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:18:LowGQ:36:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/16:21:PASS:36:0,.,.,.,16,.,.,.,0,.,.,.,.,.,.,.,5,.,. 18/18:11:PASS:53:0,.,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,.,11 0|16:41:PASS:47:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,18,.,. 12/18:21:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,9,.,.,.,.,1,11 18/18:15:PASS:60:0,.,.,.,.,.,.,.,.,.,.,.,0,.,.,.,.,.,15 0|9:45:PASS:48:24,.,.,.,.,.,.,.,.,21,.,.,.,.,.,.,.,.,. 0/0:27:PASS:62:27,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 10/18:30:PASS:49:0,.,.,.,.,.,.,.,.,.,13,.,.,.,.,.,.,.,17 0/0:17:PASS:41:17,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/9:22:PASS:49:0,.,.,.,6,.,.,.,0,16,.,.,.,.,.,.,.,.,.0/0:29:LowGQ:9:12,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/9:28:PASS:50:0,.,.,.,15,.,.,.,0,13,.,.,.,.,.,.,.,.,. 0/15:26:PASS:16:18,.,.,.,.,.,.,.,.,.,.,.,.,.,.,5,.,3,. 0|15:43:PASS:48:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,23,.,.,. 0/15:34:PASS:47:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,5,.,6,.0/0:24:PASS:48:24,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:29:PASS:63:29,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:33:LowGQ:47:32,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|15:38:PASS:44:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,15,.,.,. 0/16:30:PASS:28:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,11,.,. 0/0:36:PASS:69:35,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/15:26:PASS:18:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,4,.,.,. 0/0:23:PASS:42:23,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|18:35:PASS:34:21,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,14 0/15:30:PASS:21:22,.,.,.,.,.,.,.,.,.,.,.,.,.,.,4,.,4,../.:4:LowGQ:8:4,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:0:LowDepth;LowGQ:0:0,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 4/9:27:PASS:46:4,.,.,.,6,.,.,.,0,17,.,.,.,.,.,.,.,.,. 0/0:11:LowGQ:11:10,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 13/13:15:PASS:40:0,.,.,.,.,.,.,.,.,.,.,.,.,15,.,.,.,.,. 7/7:23:PASS:64:0,.,.,.,.,.,.,23,.,.,.,.,.,.,.,.,.,.,. 0/7:39:PASS:47:16,.,.,.,.,.,.,23,.,.,.,.,.,.,.,.,.,.,. 7/7:19:PASS:58:2,.,.,.,.,.,.,17,.,.,.,.,.,.,.,.,.,.,. 0/0:31:PASS:69:30,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:28:LowGQ:57:27,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.0|18:39:PASS:48:20,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,19 7/7:26:PASS:68:6,.,.,.,.,.,.,20,.,.,.,.,.,.,.,.,.,.,. 12/12:8:PASS:19:0,.,.,.,.,.,.,.,.,.,.,.,7,1,.,.,.,.,. 4/9:27:PASS:39:2,.,.,.,10,.,.,.,0,11,.,.,.,.,.,.,.,.,. 0/0:35:PASS:71:35,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18/18:9:PASS:33:0,.,.,.,.,.,.,.,.,.,.,.,.,.,.,2,.,.,7 ./.:2:LowGQ:0:2,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:46:LowGQ:0:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:4:PASS:9:0,.,.,.,.,.,.,.,.,4,.,.,.,.,.,.,.,.,. 0/0:27:LowGQ:59:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:18:LowGQ:0:14,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:26:LowGQ:0:11,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 18/18:40:PASS:38:0,.,.,.,.,.,.,.,.,0,.,1,.,.,.,.,.,.,39 0/0:31:LowGQ:50:29,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:28:LowGQ:54:27,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. ./.:0:.:.:.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0/0:14:DRAGENIndelHardQUAL:12:14,.,.,.,.,.,0,.,.,.,.,.,.,.,.,.,.,.,. 0/0:26:PASS:51:26,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/9:23:PASS:66:0,.,.,.,.,.,.,.,.,23,.,.,.,.,.,.,.,.,. 0|16:15:PASS:44:7,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,8,.,. 0/0:19:PASS:39:19,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 9/12:24:PASS:49:0,.,.,.,.,.,.,.,.,14,.,.,10,.,.,.,.,.,. 0/0:31:LowGQ:51:31,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 16/16:18:PASS:50:0,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,18,.,. 0/0:25:PASS:42:25,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,. 0|16:54:PASS:41:34,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,20,.,.

zhengxwen commented 2 years ago

I will check it as soon as I can.

zhengxwen commented 2 years ago

This bug only occurs when there are more than 15 different alleles at a multi-allelic site (or on a row in the VCF file).

SeqArray_1.36.2 fixes this bug. Please install SeqArray_1.36.2 from Bioconductor.