Difference between revisions of "Fabry:Mapping point mutations/Journal"
Rackersederj (talk | contribs) (→SNPdbe) |
Rackersederj (talk | contribs) (→SNPdbe) |
||
Line 15: | Line 15: | ||
== SNPdbe == |
== SNPdbe == |
||
− | The textfile containing all informations gathered in SNPdbe ([https://www.dropbox.com/s/eeboc9nratgeuj5/snps.txt snps.txt]) was downloaded after searching for the gene identifier "NP_000160" in the organism "human" ([http://www.rostlab.org/services/snpdbe/dosearch.php?id=name&val=NP_000160&organism2=human&organism1= see]). The textfile was parsed with the perl scripts [https://www.dropbox.com/s/qp7gr2nxpj91yr9/readSNPdbe.pl.txt readSNPdbe.pl] and [https://www.dropbox.com/s/ygrds0su20ojafi/SNPdbe2table.pl.txt SNPdbe2table.pl] and vizualised with the R script [https://www.dropbox.com/s/4vw4elo84ywkk3m/Hotspots_SNPdbe.R.txt Hotspots_SNPdbe.R]. For now, a SNP was considered disease causing if there was a disease listed, if it was a "N/A" it was assigned "non-disease". Later we added those SNPs that were also found in HGMD to the disease causing SNPs. |
+ | The textfile containing all informations gathered in SNPdbe ([https://www.dropbox.com/s/eeboc9nratgeuj5/snps.txt snps.txt]) was downloaded after searching for the gene identifier "NP_000160" in the organism "human" ([http://www.rostlab.org/services/snpdbe/dosearch.php?id=name&val=NP_000160&organism2=human&organism1= see]). The textfile was parsed with the perl scripts [https://www.dropbox.com/s/qp7gr2nxpj91yr9/readSNPdbe.pl.txt readSNPdbe.pl] and [https://www.dropbox.com/s/ygrds0su20ojafi/SNPdbe2table.pl.txt SNPdbe2table.pl] and vizualised with the R script [https://www.dropbox.com/s/4vw4elo84ywkk3m/Hotspots_SNPdbe.R.txt Hotspots_SNPdbe.R]. For now, a SNP was considered disease causing if there was a disease listed, if it was a "N/A" it was assigned "non-disease". Later we added those SNPs that were also found in HGMD to the disease causing SNPs (see [https://www.dropbox.com/s/ca1hkcr1ighozx3/SNPdbe_manuallyCur.Rtable SNPdbe_manuallyCur.Rtable]). |
perl [https://www.dropbox.com/s/qp7gr2nxpj91yr9/readSNPdbe.pl.txt readSNPdbe.pl] |
perl [https://www.dropbox.com/s/qp7gr2nxpj91yr9/readSNPdbe.pl.txt readSNPdbe.pl] |
Revision as of 12:14, 10 June 2012
Fabry Disease » Mapping point mutations » Journal
OMIM
In this database we searched for the gene GLA and downloaded the Table View table of the Allelic Variants (Allelic_variants.txt), which can be found in the Table of Contents. This table was parsed with the perl scripts readOMIM.pl and Omim2table.pl, which created all the output we needed for creation of the tables and statistics.
perl readOMIM.pl perl Omim2table.pl > omim.wiki
SNPedia
Since there is no special site for the gene GLA, we performed a query with the search term "Gene=GLA". These dbSNP identifiers were downloaded and mapped onto a list of all informations on all identifiers (201) found, when searching for ""snp"[SNP_CLASS] AND GLA[GENE] AND "human"[ORGN]" in the dbSNP (Flat File display). This was done with the perl script parse_dbSNP.pl, which needed the input files snp_result.txt and rsnumber.txt. We realised that the amino acid position in this representation ignores the initial Methionine in the sequence, so we had to add 1 to the read aa positions.
A wiki table was created with SNPedia2table.pl.
perl parse_dbSNP.pl perl SNPedia2table.pl > SNPedia.wiki
SNPdbe
The textfile containing all informations gathered in SNPdbe (snps.txt) was downloaded after searching for the gene identifier "NP_000160" in the organism "human" (see). The textfile was parsed with the perl scripts readSNPdbe.pl and SNPdbe2table.pl and vizualised with the R script Hotspots_SNPdbe.R. For now, a SNP was considered disease causing if there was a disease listed, if it was a "N/A" it was assigned "non-disease". Later we added those SNPs that were also found in HGMD to the disease causing SNPs (see SNPdbe_manuallyCur.Rtable).
perl readSNPdbe.pl perl SNPdbe2table.pl > SNPdbe.wiki R CMD BATCH Hotspots_SNPdbe.R
Mapping
First of all, we created the plots showing the distribution of SNPs along the sequence of the gene GLA with the R Script snp_distr.R. For this we prepared the data with map.pl
perl map.pl R CMD BATCH snp_distr.R
Afterwards we tried to vizualise disease and non-disease causing point mutations in one plot with Hotspots_all.R. We needed R tables of all databases for this, so some preparation had to be done with readHGMDforR.pl (input hgmd_snps.txt) and readdbSNPforR.pl (input dbSNP_snps.txt)
perl readHGMDforR.pl perl readdbSNPforR.pl R CMD BATCH Hotspots_all.R