git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/gmapR@101133 bc3139a8-67e5-0310-9ffc-ced21a209358
... | ... |
@@ -18,10 +18,14 @@ enum{ CODON_MINUS = -1, |
18 | 18 |
|
19 | 19 |
|
20 | 20 |
static char *codon_table[64] = |
21 |
- {"AAA", "AAC", "AAG", "AAT", "ACA", "ACC", "ACG", "ACT", "AGA", "AGC", "AGG", "AGT", "ATA", "ATC", "ATG", "ATT", |
|
22 |
- "CAA", "CAC", "CAG", "CAT", "CCA", "CCC", "CCG", "CCT", "CGA", "CGC", "CGG", "CGT", "CTA", "CTC", "CTG", "CTT", |
|
23 |
- "GAA", "GAC", "GAG", "GAT", "GCA", "GCC", "GCG", "GCT", "GGA", "GGC", "GGG", "GGT", "GTA", "GTC", "GTG", "GTT", |
|
24 |
- "TAA", "TAC", "TAG", "TAT", "TCA", "TCC", "TCG", "TCT", "TGA", "TGC", "TGG", "TGT", "TTA", "TTC", "TTG", "TTT"}; |
|
21 |
+ {"AAA", "AAC", "AAG", "AAT", "ACA", "ACC", "ACG", "ACT", |
|
22 |
+ "AGA", "AGC", "AGG", "AGT", "ATA", "ATC", "ATG", "ATT", |
|
23 |
+ "CAA", "CAC", "CAG", "CAT", "CCA", "CCC", "CCG", "CCT", |
|
24 |
+ "CGA", "CGC", "CGG", "CGT", "CTA", "CTC", "CTG", "CTT", |
|
25 |
+ "GAA", "GAC", "GAG", "GAT", "GCA", "GCC", "GCG", "GCT", |
|
26 |
+ "GGA", "GGC", "GGG", "GGT", "GTA", "GTC", "GTG", "GTT", |
|
27 |
+ "TAA", "TAC", "TAG", "TAT", "TCA", "TCC", "TCG", "TCT", |
|
28 |
+ "TGA", "TGC", "TGG", "TGT", "TTA", "TTC", "TTG", "TTT"}; |
|
25 | 29 |
|
26 | 30 |
|
27 | 31 |
|
... | ... |
@@ -288,7 +292,6 @@ read_allele_counts(unsigned char **bytes, int row, SEXP read_R, |
288 | 292 |
#ifdef DEBUG2 |
289 | 293 |
printf("Starting read_allele_counts (codon strand %d) at row %d. Total length of read_R is %d\n", strand, row, LENGTH(read_R)); |
290 | 294 |
#endif |
291 |
-// while((allele = read_char(bytes)) != '\0' && (strand == 0 || (int) allele != 255)) { |
|
292 | 295 |
while((allele = read_char(bytes)) != stop) { |
293 | 296 |
#ifdef DEBUG2 |
294 | 297 |
printf("Parsing counts for allele: %d (%s) row %d \n", allele, &allele, row ); |