Question: determine unique snp
0
gravatar for slimane.khayi
2.3 years ago by
France
slimane.khayi40 wrote:

Dear colleagues, I have tab-limited file in this format

#CHROM                  POS REF  a1  a2    b1   b2   c1  c2
NW_008246507.1  16  T   C/C C/C T/T C/C C/C T/C
NW_008246507.1  1624    A   C/C C/C C/C C/C C/C C/C
NW_008246507.1  1656    C   T/T T/T T/T T/T T/T T/T
NW_008246507.1  1666    C   T/T T/T T/T T/T T/T T/T
NW_008246507.1  1679    C   T/T T/T T/T T/T T/T T/T
NW_008246507.1  1681    G   A/A A/A A/A A/A A/A A/A
NW_008246507.1  1682    T   A/A A/A A/A A/A A/A A/A
NW_008246507.1  1695    T   C/C C/C C/C C/C C/C C/C

I want to identify the unique SNPs for each species a, b, c (not strain a1, a2, b1..),have you any python script or any idea to do this job, I am not familiar within scripting languages. Thank you in advance for your help. Sincerely.

ADD COMMENTlink modified 2.3 years ago by Pierre Lindenbaum120k • written 2.3 years ago by slimane.khayi40

Could you please clarify what is a unique SNP in your example data, I find it difficult to see (and I think you are missing a couple of line breaks, as there are currently two positions per line...) -- Would T/C for species c at NW_008246507.1:16 be what you are looking for?

ADD REPLYlink written 2.3 years ago by cschu1811.7k

Here is a capture from my file

ADD REPLYlink written 2.3 years ago by slimane.khayi40

commercial photography locations

ADD REPLYlink written 2.3 years ago by slimane.khayi40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1936 users visited in the last hour