convert_bim_allele.pl help

Question

convert_bim_allele.pl help

Opened this issue 3 years ago · 6 comments

I am using Illumina data.

rs646776 in dbSNP is C>T. In my data it is G>A
I have created a SNPTable.txt file (the line for rs646776 looks like this):
rs646776 [A/G] TOP TOP

I tried with --outtype dbsnp:
$ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype dbsnp --outfile converted_allele.bim
but I get the same G A genotypes

I'm thinking I need to do --outtype forward:
$ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype forward --outfile converted_allele.bim --strandfile my_strandfile.txt

I made my_strandfile.txt from taking the rsID and strand information from the Illumina manifest.csv
so:
rs646776 -

This seems to run successfully but when I check the now converted_allele.bim
$ grep rs646776 converted_allele.bim
1 rs646776 129.6542 109818530 G A

I'm still getting the same genotype as before?
Do you know how I convert this to the same as dbSNP e.g. C>T

Answer 1 · 2022-01-24T05:06:31.000Z

Did you use the exported data from Illumina genome studio? It seems that you are just using dbSNP data? In that case, there is no need to convert_bim_allele.pl tool. "dbsnp" is not the best out type to use, because it contains a minor fraction of reverse strand SNP information as well. So if you have TOP allele, it is best that you convert to forward strand allele for a specific genome build. Annoar's hg19_avsnp147 has no strand information in. You have to use dbsnp147 directly from UCSC genome browser as you have done. It is likely that this SNP is obsolete, or too new, so it is not included in the snp147 file. Try a more recent version to see if you can find it.

…

On Sun, Jan 23, 2022 at 8:04 PM zoeward-nz ***@***.***> wrote: I am using Illumina data. rs646776 in dbSNP is C>T. In my data it is G>A I have created a SNPTable.txt file (the line for rs646776 looks like this): rs646776 [A/G] TOP TOP I tried with --outtype dbsnp: $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype dbsnp --outfile converted_allele.bim but I get the same G A genotypes I'm thinking I need to do --outtype forward: $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype forward --outfile converted_allele.bim --strandfile my_strandfile.txt I tried to make a my_strandflle.txt by downloading the hg19_avsnp147.txt from ANNOVAR but doesn't have strand information in it?? (the line for rs646776 looks like this): 1 109818530 109818530 C T rs646776 I tried downloading the snp147.txt.gz from here: http://hgdownload.soe.ucsc.edu/goldenPath/galGal5/database/ which does appear to have strand information in it but this doesn't have my rs646776 present. Can you help with the correct strandfile?? — Reply to this email directly, view it on GitHub <#10>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABNG3OFQTOYTWSEEILBVMTDUXSQSLANCNFSM5MUCGWGQ> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

Answer 2 · 2022-01-24T19:48:04.000Z

Hi Kai, Thanks for your speedy reply as I'm totally stuck on this. I want my data to reflect the genotypes that are in dbSNP. If this is also achieved with the forward outtype then I'm happy to this. At the moment my data does NOT reflect what is in dnSNP. The example snp that I used is C>T in dbSNP and G>A in my data

…

On Mon, Jan 24, 2022, 18:06 Kai Wang ***@***.***> wrote: Did you use the exported data from Illumina genome studio? It seems that you are just using dbSNP data? In that case, there is no need to convert_bim_allele.pl tool. "dbsnp" is not the best out type to use, because it contains a minor fraction of reverse strand SNP information as well. So if you have TOP allele, it is best that you convert to forward strand allele for a specific genome build. Annoar's hg19_avsnp147 has no strand information in. You have to use dbsnp147 directly from UCSC genome browser as you have done. It is likely that this SNP is obsolete, or too new, so it is not included in the snp147 file. Try a more recent version to see if you can find it. On Sun, Jan 23, 2022 at 8:04 PM zoeward-nz ***@***.***> wrote: > I am using Illumina data. > > rs646776 in dbSNP is C>T. In my data it is G>A > I have created a SNPTable.txt file (the line for rs646776 looks like this): > rs646776 [A/G] TOP TOP > > I tried with --outtype dbsnp: > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype > dbsnp --outfile converted_allele.bim > but I get the same G A genotypes > > I'm thinking I need to do --outtype forward: > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype > forward --outfile converted_allele.bim --strandfile my_strandfile.txt > > I tried to make a my_strandflle.txt by downloading the hg19_avsnp147.txt > from ANNOVAR but doesn't have strand information in it?? (the line for > rs646776 looks like this): > 1 109818530 109818530 C T rs646776 > > I tried downloading the snp147.txt.gz from here: > http://hgdownload.soe.ucsc.edu/goldenPath/galGal5/database/ > > which does appear to have strand information in it but this doesn't have > my rs646776 present. > > Can you help with the correct strandfile?? > > — > Reply to this email directly, view it on GitHub > <#10>, or unsubscribe > < https://github.com/notifications/unsubscribe-auth/ABNG3OFQTOYTWSEEILBVMTDUXSQSLANCNFSM5MUCGWGQ > > . > You are receiving this because you are subscribed to this thread.Message > ID: ***@***.***> > — Reply to this email directly, view it on GitHub <#10 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AR6LYQ6FYTLVPYJNQKOBXRTUXTM6NANCNFSM5MUCGWGQ> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: ***@***.***>

Answer 3 · 2022-01-25T03:10:32.000Z

dbSNP has numerous errors and by definition, dbSNP's rsid is a locus identifier only (the same dbSNP record can mean a dozen completely different mutations), not a mutation identifier, so it can be problematic to match dbsnp since information in dbsnp is not identifiable. Using forward strand for a specific genome build (such as hg19 and hg38) is really the only way to identify a mutation.

…

On Mon, Jan 24, 2022 at 2:48 PM zoeward-nz ***@***.***> wrote: Hi Kai, Thanks for your speedy reply as I'm totally stuck on this. I want my data to reflect the genotypes that are in dbSNP. Of this is achieved with the forward outtype then I'm happy to this also. At the moment my data does NOT reflect what is in dnSNP. The example snp that I used is C>T in dbSNP and G>A in my data On Mon, Jan 24, 2022, 18:06 Kai Wang ***@***.***> wrote: > Did you use the exported data from Illumina genome studio? It seems that > you are just using dbSNP data? In that case, there is no need to > convert_bim_allele.pl tool. > "dbsnp" is not the best out type to use, because it contains a minor > fraction of reverse strand SNP information as well. So if you have TOP > allele, it is best that you convert to forward strand allele for a specific > genome build. > > Annoar's hg19_avsnp147 has no strand information in. You have to use > dbsnp147 directly from UCSC genome browser as you have done. It is likely > that this SNP is obsolete, or too new, so it is not included in the snp147 > file. Try a more recent version to see if you can find it. > > On Sun, Jan 23, 2022 at 8:04 PM zoeward-nz ***@***.***> wrote: > > > I am using Illumina data. > > > > rs646776 in dbSNP is C>T. In my data it is G>A > > I have created a SNPTable.txt file (the line for rs646776 looks like > this): > > rs646776 [A/G] TOP TOP > > > > I tried with --outtype dbsnp: > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype > > dbsnp --outfile converted_allele.bim > > but I get the same G A genotypes > > > > I'm thinking I need to do --outtype forward: > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top --outtype > > forward --outfile converted_allele.bim --strandfile my_strandfile.txt > > > > I tried to make a my_strandflle.txt by downloading the hg19_avsnp147.txt > > from ANNOVAR but doesn't have strand information in it?? (the line for > > rs646776 looks like this): > > 1 109818530 109818530 C T rs646776 > > > > I tried downloading the snp147.txt.gz from here: > > http://hgdownload.soe.ucsc.edu/goldenPath/galGal5/database/ > > > > which does appear to have strand information in it but this doesn't have > > my rs646776 present. > > > > Can you help with the correct strandfile?? > > > > — > > Reply to this email directly, view it on GitHub > > <#10>, or unsubscribe > > < > https://github.com/notifications/unsubscribe-auth/ABNG3OFQTOYTWSEEILBVMTDUXSQSLANCNFSM5MUCGWGQ > > > > . > > You are receiving this because you are subscribed to this thread.Message > > ID: ***@***.***> > > > > — > Reply to this email directly, view it on GitHub > <#10 (comment)>, or > unsubscribe > < https://github.com/notifications/unsubscribe-auth/AR6LYQ6FYTLVPYJNQKOBXRTUXTM6NANCNFSM5MUCGWGQ > > . > Triage notifications on the go with GitHub Mobile for iOS > < https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 > > or Android > < https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub >. > > You are receiving this because you authored the thread.Message ID: > ***@***.***> > — Reply to this email directly, view it on GitHub <#10 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABNG3OBB7DBB3HU6XILXYALUXWUIBANCNFSM5MUCGWGQ> . You are receiving this because you commented.Message ID: ***@***.***>

Answer 4 · 2022-01-25T04:05:49.000Z

Thanks Kai, So if I set convert_bim_allele.pl to --outype forward All of the genotypes in my data will be referencing the forward strand of the sequence? <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Virus-free. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

…

On Tue, Jan 25, 2022 at 4:10 PM Kai Wang ***@***.***> wrote: dbSNP has numerous errors and by definition, dbSNP's rsid is a locus identifier only (the same dbSNP record can mean a dozen completely different mutations), not a mutation identifier, so it can be problematic to match dbsnp since information in dbsnp is not identifiable. Using forward strand for a specific genome build (such as hg19 and hg38) is really the only way to identify a mutation. On Mon, Jan 24, 2022 at 2:48 PM zoeward-nz ***@***.***> wrote: > Hi Kai, > > Thanks for your speedy reply as I'm totally stuck on this. > > I want my data to reflect the genotypes that are in dbSNP. Of this is > achieved with the forward outtype then I'm happy to this also. > > At the moment my data does NOT reflect what is in dnSNP. The example snp > that I used is C>T in dbSNP and G>A in my data > > On Mon, Jan 24, 2022, 18:06 Kai Wang ***@***.***> wrote: > > > Did you use the exported data from Illumina genome studio? It seems that > > you are just using dbSNP data? In that case, there is no need to > > convert_bim_allele.pl tool. > > "dbsnp" is not the best out type to use, because it contains a minor > > fraction of reverse strand SNP information as well. So if you have TOP > > allele, it is best that you convert to forward strand allele for a > specific > > genome build. > > > > Annoar's hg19_avsnp147 has no strand information in. You have to use > > dbsnp147 directly from UCSC genome browser as you have done. It is likely > > that this SNP is obsolete, or too new, so it is not included in the > snp147 > > file. Try a more recent version to see if you can find it. > > > > On Sun, Jan 23, 2022 at 8:04 PM zoeward-nz ***@***.***> wrote: > > > > > I am using Illumina data. > > > > > > rs646776 in dbSNP is C>T. In my data it is G>A > > > I have created a SNPTable.txt file (the line for rs646776 looks like > > this): > > > rs646776 [A/G] TOP TOP > > > > > > I tried with --outtype dbsnp: > > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top > --outtype > > > dbsnp --outfile converted_allele.bim > > > but I get the same G A genotypes > > > > > > I'm thinking I need to do --outtype forward: > > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top > --outtype > > > forward --outfile converted_allele.bim --strandfile my_strandfile.txt > > > > > > I tried to make a my_strandflle.txt by downloading the > hg19_avsnp147.txt > > > from ANNOVAR but doesn't have strand information in it?? (the line for > > > rs646776 looks like this): > > > 1 109818530 109818530 C T rs646776 > > > > > > I tried downloading the snp147.txt.gz from here: > > > http://hgdownload.soe.ucsc.edu/goldenPath/galGal5/database/ > > > > > > which does appear to have strand information in it but this doesn't > have > > > my rs646776 present. > > > > > > Can you help with the correct strandfile?? > > > > > > — > > > Reply to this email directly, view it on GitHub > > > <#10>, or unsubscribe > > > < > > > https://github.com/notifications/unsubscribe-auth/ABNG3OFQTOYTWSEEILBVMTDUXSQSLANCNFSM5MUCGWGQ > > > > > > . > > > You are receiving this because you are subscribed to this > thread.Message > > > ID: ***@***.***> > > > > > > > — > > Reply to this email directly, view it on GitHub > > <#10 (comment)>, or > > unsubscribe > > < > https://github.com/notifications/unsubscribe-auth/AR6LYQ6FYTLVPYJNQKOBXRTUXTM6NANCNFSM5MUCGWGQ > > > > . > > Triage notifications on the go with GitHub Mobile for iOS > > < > https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 > > > > or Android > > < > https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub > >. > > > > You are receiving this because you authored the thread.Message ID: > > ***@***.***> > > > > — > Reply to this email directly, view it on GitHub > <#10 (comment)>, or > unsubscribe > < https://github.com/notifications/unsubscribe-auth/ABNG3OBB7DBB3HU6XILXYALUXWUIBANCNFSM5MUCGWGQ > > . > You are receiving this because you commented.Message ID: > ***@***.***> > — Reply to this email directly, view it on GitHub <#10 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AR6LYQ7CKZ7XUXS6SWWQ6CDUXYIDTANCNFSM5MUCGWGQ> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: ***@***.***>

Answer 5 · 2022-01-25T04:37:36.000Z

Yes, if the snptable file is correct for the specific genome coordinate that you are using. On Mon, Jan 24, 2022 at 11:06 PM zoeward-nz ***@***.***> wrote:

…

Thanks Kai, So if I set convert_bim_allele.pl to --outype forward All of the genotypes in my data will be referencing the forward strand of the sequence? < https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail > Virus-free. www.avast.com < https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> On Tue, Jan 25, 2022 at 4:10 PM Kai Wang ***@***.***> wrote: > dbSNP has numerous errors and by definition, dbSNP's rsid is a locus > identifier only (the same dbSNP record can mean a dozen completely > different mutations), not a mutation identifier, so it can be problematic > to match dbsnp since information in dbsnp is not identifiable. Using > forward strand for a specific genome build (such as hg19 and hg38) is > really the only way to identify a mutation. > > On Mon, Jan 24, 2022 at 2:48 PM zoeward-nz ***@***.***> wrote: > > > Hi Kai, > > > > Thanks for your speedy reply as I'm totally stuck on this. > > > > I want my data to reflect the genotypes that are in dbSNP. Of this is > > achieved with the forward outtype then I'm happy to this also. > > > > At the moment my data does NOT reflect what is in dnSNP. The example snp > > that I used is C>T in dbSNP and G>A in my data > > > > On Mon, Jan 24, 2022, 18:06 Kai Wang ***@***.***> wrote: > > > > > Did you use the exported data from Illumina genome studio? It seems > that > > > you are just using dbSNP data? In that case, there is no need to > > > convert_bim_allele.pl tool. > > > "dbsnp" is not the best out type to use, because it contains a minor > > > fraction of reverse strand SNP information as well. So if you have TOP > > > allele, it is best that you convert to forward strand allele for a > > specific > > > genome build. > > > > > > Annoar's hg19_avsnp147 has no strand information in. You have to use > > > dbsnp147 directly from UCSC genome browser as you have done. It is > likely > > > that this SNP is obsolete, or too new, so it is not included in the > > snp147 > > > file. Try a more recent version to see if you can find it. > > > > > > On Sun, Jan 23, 2022 at 8:04 PM zoeward-nz ***@***.***> wrote: > > > > > > > I am using Illumina data. > > > > > > > > rs646776 in dbSNP is C>T. In my data it is G>A > > > > I have created a SNPTable.txt file (the line for rs646776 looks like > > > this): > > > > rs646776 [A/G] TOP TOP > > > > > > > > I tried with --outtype dbsnp: > > > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top > > --outtype > > > > dbsnp --outfile converted_allele.bim > > > > but I get the same G A genotypes > > > > > > > > I'm thinking I need to do --outtype forward: > > > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top > > --outtype > > > > forward --outfile converted_allele.bim --strandfile my_strandfile.txt > > > > > > > > I tried to make a my_strandflle.txt by downloading the > > hg19_avsnp147.txt > > > > from ANNOVAR but doesn't have strand information in it?? (the line > for > > > > rs646776 looks like this): > > > > 1 109818530 109818530 C T rs646776 > > > > > > > > I tried downloading the snp147.txt.gz from here: > > > > http://hgdownload.soe.ucsc.edu/goldenPath/galGal5/database/ > > > > > > > > which does appear to have strand information in it but this doesn't > > have > > > > my rs646776 present. > > > > > > > > Can you help with the correct strandfile?? > > > > > > > > — > > > > Reply to this email directly, view it on GitHub > > > > <#10>, or unsubscribe > > > > < > > > > > > https://github.com/notifications/unsubscribe-auth/ABNG3OFQTOYTWSEEILBVMTDUXSQSLANCNFSM5MUCGWGQ > > > > > > > > . > > > > You are receiving this because you are subscribed to this > > thread.Message > > > > ID: ***@***.***> > > > > > > > > > > — > > > Reply to this email directly, view it on GitHub > > > <#10 (comment)>, > or > > > unsubscribe > > > < > > > https://github.com/notifications/unsubscribe-auth/AR6LYQ6FYTLVPYJNQKOBXRTUXTM6NANCNFSM5MUCGWGQ > > > > > > . > > > Triage notifications on the go with GitHub Mobile for iOS > > > < > > > https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 > > > > > > or Android > > > < > > > https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub > > >. > > > > > > You are receiving this because you authored the thread.Message ID: > > > ***@***.***> > > > > > > > — > > Reply to this email directly, view it on GitHub > > <#10 (comment)>, or > > unsubscribe > > < > https://github.com/notifications/unsubscribe-auth/ABNG3OBB7DBB3HU6XILXYALUXWUIBANCNFSM5MUCGWGQ > > > > . > > You are receiving this because you commented.Message ID: > > ***@***.***> > > > > — > Reply to this email directly, view it on GitHub > <#10 (comment)>, or > unsubscribe > < https://github.com/notifications/unsubscribe-auth/AR6LYQ7CKZ7XUXS6SWWQ6CDUXYIDTANCNFSM5MUCGWGQ > > . > Triage notifications on the go with GitHub Mobile for iOS > < https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 > > or Android > < https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub >. > > You are receiving this because you authored the thread.Message ID: > ***@***.***> > — Reply to this email directly, view it on GitHub <#10 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABNG3OAQ5GUMWPQA5UXFHELUXYOSZANCNFSM5MUCGWGQ> . You are receiving this because you commented.Message ID: ***@***.***>

Answer 6 · 2022-01-25T04:44:24.000Z

Then I don't think it has worked? After running convert_bim_allele.pl with --outtype forward I get (using one snp as an example): 19 rs1122608 31.58171 11163601 A C On hg19 position 19:11163601 the forward strand is a G <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Virus-free. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

…

On Tue, Jan 25, 2022 at 5:37 PM Kai Wang ***@***.***> wrote: Yes, if the snptable file is correct for the specific genome coordinate that you are using. On Mon, Jan 24, 2022 at 11:06 PM zoeward-nz ***@***.***> wrote: > Thanks Kai, > > So if I set convert_bim_allele.pl to > --outype forward > All of the genotypes in my data will be referencing the forward strand of > the sequence? > > > > > > > > > > > < > https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail > > > Virus-free. > www.avast.com > < > https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail > > > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> > > On Tue, Jan 25, 2022 at 4:10 PM Kai Wang ***@***.***> wrote: > > > dbSNP has numerous errors and by definition, dbSNP's rsid is a locus > > identifier only (the same dbSNP record can mean a dozen completely > > different mutations), not a mutation identifier, so it can be problematic > > to match dbsnp since information in dbsnp is not identifiable. Using > > forward strand for a specific genome build (such as hg19 and hg38) is > > really the only way to identify a mutation. > > > > On Mon, Jan 24, 2022 at 2:48 PM zoeward-nz ***@***.***> wrote: > > > > > Hi Kai, > > > > > > Thanks for your speedy reply as I'm totally stuck on this. > > > > > > I want my data to reflect the genotypes that are in dbSNP. Of this is > > > achieved with the forward outtype then I'm happy to this also. > > > > > > At the moment my data does NOT reflect what is in dnSNP. The example > snp > > > that I used is C>T in dbSNP and G>A in my data > > > > > > On Mon, Jan 24, 2022, 18:06 Kai Wang ***@***.***> wrote: > > > > > > > Did you use the exported data from Illumina genome studio? It seems > > that > > > > you are just using dbSNP data? In that case, there is no need to > > > > convert_bim_allele.pl tool. > > > > "dbsnp" is not the best out type to use, because it contains a minor > > > > fraction of reverse strand SNP information as well. So if you have > TOP > > > > allele, it is best that you convert to forward strand allele for a > > > specific > > > > genome build. > > > > > > > > Annoar's hg19_avsnp147 has no strand information in. You have to use > > > > dbsnp147 directly from UCSC genome browser as you have done. It is > > likely > > > > that this SNP is obsolete, or too new, so it is not included in the > > > snp147 > > > > file. Try a more recent version to see if you can find it. > > > > > > > > On Sun, Jan 23, 2022 at 8:04 PM zoeward-nz ***@***.***> wrote: > > > > > > > > > I am using Illumina data. > > > > > > > > > > rs646776 in dbSNP is C>T. In my data it is G>A > > > > > I have created a SNPTable.txt file (the line for rs646776 looks > like > > > > this): > > > > > rs646776 [A/G] TOP TOP > > > > > > > > > > I tried with --outtype dbsnp: > > > > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top > > > --outtype > > > > > dbsnp --outfile converted_allele.bim > > > > > but I get the same G A genotypes > > > > > > > > > > I'm thinking I need to do --outtype forward: > > > > > $ perl convert_bim_allele.pl my.bim SNPTable.txt --intype top > > > --outtype > > > > > forward --outfile converted_allele.bim --strandfile > my_strandfile.txt > > > > > > > > > > I tried to make a my_strandflle.txt by downloading the > > > hg19_avsnp147.txt > > > > > from ANNOVAR but doesn't have strand information in it?? (the line > > for > > > > > rs646776 looks like this): > > > > > 1 109818530 109818530 C T rs646776 > > > > > > > > > > I tried downloading the snp147.txt.gz from here: > > > > > http://hgdownload.soe.ucsc.edu/goldenPath/galGal5/database/ > > > > > > > > > > which does appear to have strand information in it but this doesn't > > > have > > > > > my rs646776 present. > > > > > > > > > > Can you help with the correct strandfile?? > > > > > > > > > > — > > > > > Reply to this email directly, view it on GitHub > > > > > <#10>, or unsubscribe > > > > > < > > > > > > > > > > https://github.com/notifications/unsubscribe-auth/ABNG3OFQTOYTWSEEILBVMTDUXSQSLANCNFSM5MUCGWGQ > > > > > > > > > > . > > > > > You are receiving this because you are subscribed to this > > > thread.Message > > > > > ID: ***@***.***> > > > > > > > > > > > > > — > > > > Reply to this email directly, view it on GitHub > > > > <#10 (comment) >, > > or > > > > unsubscribe > > > > < > > > > > > https://github.com/notifications/unsubscribe-auth/AR6LYQ6FYTLVPYJNQKOBXRTUXTM6NANCNFSM5MUCGWGQ > > > > > > > > . > > > > Triage notifications on the go with GitHub Mobile for iOS > > > > < > > > > > > https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 > > > > > > > > or Android > > > > < > > > > > > https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub > > > >. > > > > > > > > You are receiving this because you authored the thread.Message ID: > > > > ***@***.***> > > > > > > > > > > — > > > Reply to this email directly, view it on GitHub > > > <#10 (comment)>, > or > > > unsubscribe > > > < > > > https://github.com/notifications/unsubscribe-auth/ABNG3OBB7DBB3HU6XILXYALUXWUIBANCNFSM5MUCGWGQ > > > > > > . > > > You are receiving this because you commented.Message ID: > > > ***@***.***> > > > > > > > — > > Reply to this email directly, view it on GitHub > > <#10 (comment)>, or > > unsubscribe > > < > https://github.com/notifications/unsubscribe-auth/AR6LYQ7CKZ7XUXS6SWWQ6CDUXYIDTANCNFSM5MUCGWGQ > > > > . > > Triage notifications on the go with GitHub Mobile for iOS > > < > https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 > > > > or Android > > < > https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub > >. > > > > You are receiving this because you authored the thread.Message ID: > > ***@***.***> > > > > — > Reply to this email directly, view it on GitHub > <#10 (comment)>, or > unsubscribe > < https://github.com/notifications/unsubscribe-auth/ABNG3OAQ5GUMWPQA5UXFHELUXYOSZANCNFSM5MUCGWGQ > > . > You are receiving this because you commented.Message ID: > ***@***.***> > — Reply to this email directly, view it on GitHub <#10 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AR6LYQ4FKPZM7VSZP7QMYXDUXYSJ5ANCNFSM5MUCGWGQ> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: ***@***.***>