hmm, my program (
http://magictour.free.fr/seq/seq1.exe)
with switch t=1 , still ignores positions where no full
nucleotide (A,C,G,T) mutation occurs.
It makes some sense to list those positions in caes like this.
(but not so much in Alberta-swine or Singapore/2009 or such)
So I changed it for you, honey.
I hope it works correctly and I didn't include another bug
Code:
00000000000000000000000000011111
00011223334455666677788889901224
16934080296813056811204680585130
64152723462966783156619553708538
-codon-position-----------------------------111 1 211 12 1 1
---Index------------------------------------GAGCAGCTCGTAAGGTACGAGACGTTAGGATC
1 >A/Mexflu/index/2009/02/01 ................................
2 >A/New Jersey/AF21791,2011/01/13 ....G..C.A...AAATT.....ACGGA.---
3 >A/Florida/AF21771,2011/01/12 A.A.G....A...AAA.T.....A.GGA.---
4 >A/North Carolina/AF21796,2011/01/18 AG..G...TA...AAA.T....TA.GGA.---
5 >A/South Carolina/AF21803,2011/01/14 AG..G...TA...AAA.T....TA.GGAA---
6 >A/Arizona/AF21768,2011/01/07 AG..G...TA..GAAA.T..A.TA.GGA.---
7 >A/Mexico/InDRE1945,2011/03/ ...TGAT..AAC.AAA.TA....A.GGA.GCT
8 >A/Mexico/InDRE1946,2011/03/22 ...TGAT..AAC.AAA.T.....A.GGA.GCT
9 >A/Mexico/InDRE1947,2011/03/23 ...TGA...AAC.AAA.T...G.A.GGA.GCT
10 >A/Mexico/InDRE2192,2011/03/28 ---------AAC.AAA.T...-----------
11 >A/Mexico/InDRE2195,2011/03/28 ---------AAC.AAA.TRR.-----------
12 >A/Mexico/InDRE2197,2011/03/28 ---------AAC.AAA.T...-----------
13 >A/Mexico/InDRE2200,2011/03/29 ---------AAC.AAA.T...-----------
14 >A/Mexico/InDRE2222,2011/03/31 ---------AAC.AAA.T...-----------
---Index------------------------------------GAGCAGCTCGTAAGGTACGAGACGTTAGGATC
-codon-position-----------------------------111 1 211 12 1 1
00000000000000000000000000011111
00011223334455666677788889901224
16934080296813056811204680585130
64152723462966783156619553708538
.
.