This is a test version of Biostars. For the public version, visit https://www.biostars.org.
prostate cancer detection question

SNP: rs10993994 linked to prostate cancer, info here says: at chromosome 10, position 46046326:

if "CC" ==> Mag=0 (i.e. good)

else if "CT" ==> Mag=1.5 & increased prostate cancer risk (odds ratio 1.2)

else if "TT" ==> Mag=2 & increased prostate cancer risk (odds ratio 1.6)

http://www.snpedia.com/index.php/Rs10993994

i downloaded chromosome 10 from this link:

http://hgdownload.cse.ucsc.edu/goldenpath/hg19/chromosomes/chr10.fa.gz

and removed the first line from the file (i.e. ">chr10") and removed every space occurs in the end of the lines. then using command line(bash):

dd bs=1 skip=46046326 count=2 if=InputFastaFile of=OutputFastaFile

this line extracts the two nucleotides in the position 46046326. but when i open the OutputFastaFile i find: "CA" not CC or CT or TT !!! am i making something wrong ?!

snp dna phenotype genotype genome

0 answers

No answers yet.

Log in to answer this question.