Originally posted by KaiYe
View Post
Thanks,
Boel
You are currently viewing the SEQanswers forums as a guest, which limits your access. Click here to register now, and join the discussion
cat XXXX_D.vcf | grep DEL | awk '{print $1,"\t",$2,"\t&&&"$4}'| sed -e 's/&&&./-/'> XXXX_D.dindel.del cat XXXX_SI.vcf | grep INS | awk '{print $1,"\t",$2,"\t&&&"$5}'| sed -e 's/&&&./+/' > XXXX_SI.dindel.ins cat XXXX_D.dindel.del XXXX_SI.dindel.ins > XXXX.dindel.pindel
cat XXXX_D.vcf XXXX_INV.vcf XXXX_SI.vcf XXXX_TD.vcf | awk '/DEL/ {pos=$2+1; print $1,"\t",pos,"\t-"substr($4,2)}; /RPL|INV/ {pos=$2+1; print $1,"\t",pos,"\t-"substr($4,2)"\n"$1,"\t",pos,"\t+"substr($5,2)}; /INS/ {pos=$2+1; print $1"\t"pos"\t+"substr($5,2)}' > XXXX-all.pindel.dindel
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 05-10-2024, 06:35 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
05-10-2024, 06:35 AM
|
||
Started by seqadmin, 05-09-2024, 02:46 PM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
05-09-2024, 02:46 PM
|
||
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
||
Started by seqadmin, 05-06-2024, 07:17 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
05-06-2024, 07:17 AM
|
Comment