- fmt.Fprintf(annow, "%s,%d,%d,%s,%s,%d,%s,%s\n", fields[0], incol+startcol/2, tileVariant, hgvsID, seqname, pos, refseq, fields[7])
- hgvscols := hgvs[hgvsID]
- if hgvscols[0] == nil {
- hgvscols = [2][]int16{make([]int16, len(cgnames)), make([]int16, len(cgnames))}
- hgvs[hgvsID] = hgvscols
+ hgvsColPair := hgvsCols[hgvsID]
+ if hgvsColPair[0] == nil {
+ // values in new columns start
+ // out as -1 ("no data yet")
+ // or 0 ("=ref") here, may
+ // change to 1 ("hgvs variant
+ // present") below, either on
+ // this line or a future line.
+ hgvsColPair = [2][]int16{make([]int16, len(cgnames)), make([]int16, len(cgnames))}
+ rt, ok := reftile[tagID(tag)]
+ if !ok {
+ err = fmt.Errorf("bug: seeing annotations for tag %d, but it has no reftile entry", tag)
+ return 1
+ }
+ for ph := 0; ph < 2; ph++ {
+ for row := 0; row < rows; row++ {
+ v := chunk[row*chunkcols+incol*2+ph]
+ if tileVariantID(v) == rt.variant {
+ hgvsColPair[ph][row] = 0
+ } else {
+ hgvsColPair[ph][row] = -1
+ }
+ }
+ }
+ hgvsCols[hgvsID] = hgvsColPair
+ hgvsref := hgvs.Variant{
+ Position: pos,
+ Ref: string(refseq),
+ New: string(refseq),
+ }
+ fmt.Fprintf(annow, "%d,%d,%d,%s:g.%s,%s,%d,%s,%s,%s\n", tag, incol+startcol/2, rt.variant, seqname, hgvsref.String(), seqname, pos, refseq, refseq, fields[8])