Synonmous And Non Synonmous Snps
2
0
Entering edit mode
13.8 years ago
Payal ▴ 150

Dear friends,

I've compared the two sequences and i ve got the locations of SNPs between the two sequences. for eg:

Sequence 1: ---------GCCTGCTGCTGGCCGGGCGGGGGACGGGG------------CGGGACCGGAGCCGGAGCTGCGGGGCGCACCGGCTAGA

sequence 2: CAGAGCTGCGCCTGCTGCTGGCCGGGCGGGGGACGGGG------------

o/p- 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109

now how can i distinguish which are synonymous and which are non synonymous SNPs from the given output? (NOTE: only for example i ve given two sequences with less nucleotides, but there are multiple sequences with thousands of nucleotide)

snp perl non • 4.0k views
ADD COMMENT
1
Entering edit mode

Answers from this previous question might be helpful: Synonymous And Non-Synonymous Snps

ADD REPLY
2
Entering edit mode
13.8 years ago
Chris ★ 1.6k

synonymous/non-syn. SNPs refer to protein coding regions. The latter alters the aminoacid while the first does not. So first of all you should be sure that the genomic sequences above are in exonic regions. Apart from that, a SNP is defined as such when it happens with a certain frequency in a population (not sure about the threshold though).

btw, I don't understand the dashes in your sequence. Does this represent an alignment?

ADD COMMENT
0
Entering edit mode

ya -- represents the alignment

ADD REPLY
1
Entering edit mode
13.8 years ago

Simply comparing two sequences of unknown quality makes it difficult to identify SNPs. If the 2 seqs are single reads from a high-throughput device, I would be hesitant to call any polymorphisms. If the above are high-quality and/or aligned to a reference genome, then you can begin to call SNPs and map them to a genome and its encoded genes.

So, with what you present above, you can perform a translated BLAST search (BLASTX) with seq as query against a protein database. The result should indicate if your seq and its SNPs fall into protein-coding region(s) and if so, if the SNP alters that protein sequence (non-synonymous) or not (synonymous).

ADD COMMENT

Login before adding your answer.

Traffic: 2597 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6