"ref2genome": &ref2genome{},
"vcf2fasta": &vcf2fasta{},
"import": &importer{},
+ "export-hgvs": &exportHGVS{},
"export-numpy": &exportNumpy{},
"filter": &filterer{},
"build-docker-image": &buildDockerImage{},
--- /dev/null
+package main
+
+import (
+ "bufio"
+ "bytes"
+ "context"
+ "errors"
+ "flag"
+ "fmt"
+ "io"
+ "net/http"
+ _ "net/http/pprof"
+ "os"
+ "sort"
+ "strings"
+ "sync"
+ "time"
+
+ "git.arvados.org/arvados.git/sdk/go/arvados"
+ "github.com/arvados/lightning/hgvs"
+ log "github.com/sirupsen/logrus"
+)
+
+type exportHGVS struct {
+}
+
+func (cmd *exportHGVS) RunCommand(prog string, args []string, stdin io.Reader, stdout, stderr io.Writer) int {
+ var err error
+ defer func() {
+ if err != nil {
+ fmt.Fprintf(stderr, "%s\n", err)
+ }
+ }()
+ flags := flag.NewFlagSet("", flag.ContinueOnError)
+ flags.SetOutput(stderr)
+ pprof := flags.String("pprof", "", "serve Go profile data at http://`[addr]:port`")
+ runlocal := flags.Bool("local", false, "run on local host (default: run in an arvados container)")
+ projectUUID := flags.String("project", "", "project `UUID` for output data")
+ priority := flags.Int("priority", 500, "container request priority")
+ refname := flags.String("ref", "", "reference genome `name`")
+ inputFilename := flags.String("i", "-", "input `file` (library)")
+ outputFilename := flags.String("o", "-", "fasta output `file`")
+ pick := flags.String("pick", "", "`name` of single genome to export")
+ err = flags.Parse(args)
+ if err == flag.ErrHelp {
+ err = nil
+ return 0
+ } else if err != nil {
+ return 2
+ }
+
+ if *pprof != "" {
+ go func() {
+ log.Println(http.ListenAndServe(*pprof, nil))
+ }()
+ }
+
+ if !*runlocal {
+ if *outputFilename != "-" {
+ err = errors.New("cannot specify output file in container mode: not implemented")
+ return 1
+ }
+ runner := arvadosContainerRunner{
+ Name: "lightning export-hgvs",
+ Client: arvados.NewClientFromEnv(),
+ ProjectUUID: *projectUUID,
+ RAM: 128000000000,
+ VCPUs: 2,
+ Priority: *priority,
+ }
+ err = runner.TranslatePaths(inputFilename)
+ if err != nil {
+ return 1
+ }
+ runner.Args = []string{"export-hgvs", "-local=true", "-pick", *pick, "-ref", *refname, "-i", *inputFilename, "-o", "/mnt/output/export.csv"}
+ var output string
+ output, err = runner.Run()
+ if err != nil {
+ return 1
+ }
+ fmt.Fprintln(stdout, output+"/export.csv")
+ return 0
+ }
+
+ input, err := os.Open(*inputFilename)
+ if err != nil {
+ return 1
+ }
+ defer input.Close()
+
+ // Error out early if seeking doesn't work on the input file.
+ _, err = input.Seek(0, io.SeekEnd)
+ if err != nil {
+ return 1
+ }
+ _, err = input.Seek(0, io.SeekStart)
+ if err != nil {
+ return 1
+ }
+
+ var mtx sync.Mutex
+ var cgs []CompactGenome
+ var tilelib tileLibrary
+ err = tilelib.LoadGob(context.Background(), input, func(cg CompactGenome) {
+ if *pick != "" && *pick != cg.Name {
+ return
+ }
+ log.Debugf("export: pick %q", cg.Name)
+ mtx.Lock()
+ defer mtx.Unlock()
+ cgs = append(cgs, cg)
+ })
+ if err != nil {
+ return 1
+ }
+ sort.Slice(cgs, func(i, j int) bool { return cgs[i].Name < cgs[j].Name })
+ log.Printf("export: pick %q => %d genomes", *pick, len(cgs))
+
+ refseq, ok := tilelib.refseqs[*refname]
+ if !ok {
+ err = fmt.Errorf("reference name %q not found in input; have %v", *refname, func() (names []string) {
+ for name := range tilelib.refseqs {
+ names = append(names, name)
+ }
+ return
+ }())
+ return 1
+ }
+
+ _, err = input.Seek(0, io.SeekStart)
+ if err != nil {
+ return 1
+ }
+
+ var output io.WriteCloser
+ if *outputFilename == "-" {
+ output = nopCloser{stdout}
+ } else {
+ output, err = os.OpenFile(*outputFilename, os.O_CREATE|os.O_WRONLY, 0777)
+ if err != nil {
+ return 1
+ }
+ defer output.Close()
+ }
+ bufw := bufio.NewWriter(output)
+ err = cmd.export(bufw, input, tilelib.taglib.keylen, refseq, cgs)
+ if err != nil {
+ return 1
+ }
+ err = bufw.Flush()
+ if err != nil {
+ return 1
+ }
+ err = output.Close()
+ if err != nil {
+ return 1
+ }
+ err = input.Close()
+ if err != nil {
+ return 1
+ }
+ return 0
+}
+
+func (cmd *exportHGVS) export(out io.Writer, librdr io.Reader, taglen int, refseq map[string][]tileLibRef, cgs []CompactGenome) error {
+ need := map[tileLibRef]bool{}
+ var seqnames []string
+ for seqname, librefs := range refseq {
+ seqnames = append(seqnames, seqname)
+ for _, libref := range librefs {
+ need[libref] = true
+ }
+ }
+ sort.Strings(seqnames)
+
+ for _, cg := range cgs {
+ for i, variant := range cg.Variants {
+ if variant == 0 {
+ continue
+ }
+ libref := tileLibRef{Tag: tagID(i / 2), Variant: variant}
+ need[libref] = true
+ }
+ }
+
+ log.Infof("export: loading %d tile variants", len(need))
+ tileVariant := map[tileLibRef]TileVariant{}
+ err := DecodeLibrary(librdr, func(ent *LibraryEntry) error {
+ for _, tv := range ent.TileVariants {
+ libref := tileLibRef{Tag: tv.Tag, Variant: tv.Variant}
+ if need[libref] {
+ tileVariant[libref] = tv
+ }
+ }
+ return nil
+ })
+ if err != nil {
+ return err
+ }
+
+ log.Infof("export: loaded %d tile variants", len(tileVariant))
+ var missing []tileLibRef
+ for libref := range need {
+ if _, ok := tileVariant[libref]; !ok {
+ missing = append(missing, libref)
+ }
+ }
+ if len(missing) > 0 {
+ if limit := 100; len(missing) > limit {
+ log.Warnf("first %d missing tiles: %v", limit, missing[:limit])
+ } else {
+ log.Warnf("missing tiles: %v", missing)
+ }
+ return fmt.Errorf("%d needed tiles are missing from library", len(missing))
+ }
+
+ refpos := 0
+ for _, seqname := range seqnames {
+ variantAt := map[int][]hgvs.Variant{} // variantAt[chromOffset][genomeIndex*2+phase]
+ for refstep, libref := range refseq[seqname] {
+ reftile := tileVariant[libref]
+ for cgidx, cg := range cgs {
+ for phase := 0; phase < 2; phase++ {
+ if len(cg.Variants) <= int(libref.Tag)*2+phase {
+ continue
+ }
+ variant := cg.Variants[int(libref.Tag)*2+phase]
+ if variant == 0 {
+ continue
+ }
+ genometile := tileVariant[tileLibRef{Tag: libref.Tag, Variant: variant}]
+ if variant == libref.Variant {
+ continue
+ }
+ refSequence := reftile.Sequence
+ // If needed, extend the
+ // reference sequence up to
+ // the tag at the end of the
+ // genometile sequence.
+ refstepend := refstep + 1
+ for refstepend < len(refseq[seqname]) && len(refSequence) >= taglen && !bytes.EqualFold(refSequence[len(refSequence)-taglen:], genometile.Sequence[len(genometile.Sequence)-taglen:]) {
+ if &refSequence[0] == &reftile.Sequence[0] {
+ refSequence = append([]byte(nil), refSequence...)
+ }
+ refSequence = append(refSequence, tileVariant[refseq[seqname][refstepend]].Sequence...)
+ refstepend++
+ }
+ vars, _ := hgvs.Diff(strings.ToUpper(string(refSequence)), strings.ToUpper(string(genometile.Sequence)), time.Second)
+ for _, v := range vars {
+ v.Position += refpos
+ log.Debugf("%s seq %s phase %d tag %d tile diff %s\n", cg.Name, seqname, phase, libref.Tag, v.String())
+ varslice := variantAt[v.Position]
+ if varslice == nil {
+ varslice = make([]hgvs.Variant, len(cgs)*2)
+ }
+ varslice[cgidx*2+phase] = v
+ variantAt[v.Position] = varslice
+ }
+ }
+ }
+ refpos += len(reftile.Sequence) - taglen
+
+ // Flush entries from variantAt that are
+ // behind refpos. Flush all entries if this is
+ // the last reftile of the path/chromosome.
+ var flushpos []int
+ lastrefstep := refstep == len(refseq[seqname])-1
+ for pos := range variantAt {
+ if lastrefstep || pos <= refpos {
+ flushpos = append(flushpos, pos)
+ }
+ }
+ sort.Slice(flushpos, func(i, j int) bool { return flushpos[i] < flushpos[j] })
+ for _, pos := range flushpos {
+ varslice := variantAt[pos]
+ delete(variantAt, pos)
+ for i := range varslice {
+ if varslice[i].Position == 0 {
+ varslice[i].Position = pos
+ }
+ }
+ for i := 0; i < len(cgs); i++ {
+ if i > 0 {
+ out.Write([]byte{'\t'})
+ }
+ var1, var2 := varslice[i*2], varslice[i*2+1]
+ if var1.Position == 0 && var2.Position == 0 {
+ out.Write([]byte{'.'})
+ } else if var1 == var2 {
+ fmt.Fprintf(out, "%s:g.%s", seqname, var1.String())
+ } else {
+ if var1.Position == 0 {
+ var1.Position = pos
+ }
+ if var2.Position == 0 {
+ var2.Position = pos
+ }
+ fmt.Fprintf(out, "%s:g.[%s];[%s]", seqname, var1.String(), var2.String())
+ }
+ }
+ out.Write([]byte{'\n'})
+ }
+ }
+ }
+ return nil
+}
Variants []tileVariantID
}
+type CompactSequence struct {
+ Name string
+ TileSequences map[string][]tileLibRef
+}
+
type TileVariant struct {
Tag tagID
Variant tileVariantID
}
type LibraryEntry struct {
- TagSet [][]byte
- CompactGenomes []CompactGenome
- TileVariants []TileVariant
+ TagSet [][]byte
+ CompactGenomes []CompactGenome
+ CompactSequences []CompactSequence
+ TileVariants []TileVariant
}
func ReadCompactGenomes(rdr io.Reader) ([]CompactGenome, error) {
func (v *Variant) String() string {
switch {
+ case len(v.New) == 0 && len(v.Ref) == 0:
+ return fmt.Sprintf("%d=", v.Position)
case len(v.New) == 0 && len(v.Ref) == 1:
return fmt.Sprintf("%ddel", v.Position)
case len(v.New) == 0:
log.Printf("%s starting", infile)
defer log.Printf("%s done", infile)
tseqs, err := cmd.tileFasta(tilelib, infile)
- var kept, dropped int
- variants[0], kept, dropped = tseqs.Variants()
- variants[1] = variants[0]
- log.Printf("%s found %d unique tags plus %d repeats", infile, kept, dropped)
- return err
+ if err != nil {
+ return err
+ }
+ totlen := 0
+ for _, tseq := range tseqs {
+ totlen += len(tseq)
+ }
+ log.Printf("%s tiled %d seqs, total len %d", infile, len(tseqs), totlen)
+ return cmd.encoder.Encode(LibraryEntry{
+ CompactSequences: []CompactSequence{{Name: infile, TileSequences: tseqs}},
+ })
}
+ // Don't write out a CompactGenomes entry
+ continue
} else if vcfFilenameRe.MatchString(infile) {
for phase := 0; phase < 2; phase++ {
phase := phase
if len(errs) > 0 {
return
}
- ntags := len(variants[0])
- if ntags < len(variants[1]) {
- ntags = len(variants[1])
- }
- flat := make([]tileVariantID, ntags*2)
- for i := 0; i < ntags; i++ {
- for hap := 0; hap < 2; hap++ {
- if i < len(variants[hap]) {
- flat[i*2+hap] = variants[hap][i]
- }
- }
- }
err := cmd.encoder.Encode(LibraryEntry{
- CompactGenomes: []CompactGenome{{Name: infile, Variants: flat}},
+ CompactGenomes: []CompactGenome{{Name: infile, Variants: flatten(variants)}},
})
if err != nil {
select {
}
return
}
+
+func flatten(variants [][]tileVariantID) []tileVariantID {
+ ntags := 0
+ for _, v := range variants {
+ if ntags < len(v) {
+ ntags = len(v)
+ }
+ }
+ flat := make([]tileVariantID, ntags*2)
+ for i := 0; i < ntags; i++ {
+ for hap := 0; hap < 2; hap++ {
+ if i < len(variants[hap]) {
+ flat[i*2+hap] = variants[hap][i]
+ }
+ }
+ }
+ return flat
+}
cmd.tilelib = &tileLibrary{
encoder: encoder,
includeNoCalls: true,
- onLoadGenome: func(cg CompactGenome) {
- err := encoder.Encode(LibraryEntry{CompactGenomes: []CompactGenome{cg}})
- if err != nil {
- cmd.setError(err)
- cancel()
- }
- },
}
cmd.mapped = map[string]map[tileLibRef]tileVariantID{}
"bytes"
"fmt"
"io"
+ "io/ioutil"
"os"
"sync"
c.Logf("len(merged) %d", merged.Len())
statsout := &bytes.Buffer{}
- code = (&stats{}).RunCommand("lightning stats", []string{"-local"}, merged, statsout, os.Stderr)
+ code = (&stats{}).RunCommand("lightning stats", []string{"-local"}, bytes.NewReader(merged.Bytes()), statsout, os.Stderr)
c.Check(code, check.Equals, 0)
c.Check(statsout.Len() > 0, check.Equals, true)
c.Logf("%s", statsout.String())
+
+ c.Check(ioutil.WriteFile(tmpdir+"/merged.gob", merged.Bytes(), 0666), check.IsNil)
+
+ hgvsout := &bytes.Buffer{}
+ code = (&exportHGVS{}).RunCommand("lightning export-hgvs", []string{"-local", "-ref", "testdata/ref.fasta", "-i", tmpdir + "/merged.gob"}, bytes.NewReader(nil), hgvsout, os.Stderr)
+ c.Check(code, check.Equals, 0)
+ c.Check(hgvsout.Len() > 0, check.Equals, true)
+ c.Logf("%s", hgvsout.String())
+ c.Check(hgvsout.String(), check.Equals, `chr1:g.[41_42delinsAA];[41=]
+chr1:g.[161=];[161A>T]
+chr1:g.[178=];[178A>T]
+chr1:g.222_224del
+chr1:g.[302=];[302_305delinsAAAA]
+chr2:g.[813_826del];[813=]
+chr2:g.[830_841delinsAA];[830=]
+chr2:g.[887C>A];[887=]
+chr2:g.[1042_1044del];[1042=]
+chr2:g.[1043=];[1043_1044delinsAA]
+`)
}
>chr1
-ggcgtctacctcgagaagccccgacctctgaataagatctaagaacatctcaagggattgtgtatcttgttgggtgtacgcgcgccagcccgcagcatta
+ggcgtctacctcgagaagccccgacctctgaataagatctAAgaacatctcaagggattgtgtatcttgttgggtgtacgcgcgccagcccgcagcatta
ggagaactgtgctccgccttcaga
ccccttgggtaaaatgccgcgcaatatgttgattacacttgctgcccatctgaaaggtcgccttatcaatcctatgctgaatgccctctaaggagtt
acacatgctagcgcgtcggggtgg
tttttatagcagtagcggttgcgataatgcgcactaaggtggccataacttagccacacagactgcgacctcggtgtcaatctttaggcgatgactagtg
gttattaataataacttatcatca
cttccggggaaactagcgtaaaaaccgccgtcgcagtataccaccttatacctgtgccactcaatacaggagttgctcagcccaagacaccaacgactaa
-gctctcaaaccttgtatttttctt
-atgtcctcgcttggattatggggactcgcagtaaatgacaaccgtgcccgctggccctggccggaccgcgcgtccgtaagtagcgatcgagtcctgcagt
+gctctcaaaccttgtaAAAAActt
+atgtcctcgcttggattatggggactcgcagtaaatgacaacAgtgcccgctggccctggccggaccgcgcgtccgtaagtagcgatcgagtcctgcagt
aaaactgatccaaaaaaaatacaa
acaccacgtttataccctgttgagtcgccacgtaacgtattaacgtatgaacggcccggttttcttatctcgcaccgctagtttgccctggcggtcg
cctatgagtcaatcctattttcaa
>chr1
ggcgtctacctcgagaagccccgacctctgaataagatctttgaacatctcaagggattgtgtatcttgttgggtgtacgcgcgccagcccgcagcatta
ggagaactgtgctccgccttcaga
-ccccttgggtaaaatgccgcgcaatatgttgattacacttgctgcccatctgaaaggtcgccttatcaatcctatgctgaatgccctctaaggagtt
+ccccttgggtaaaatgccgcgcaatatgttgattacTcttgctgcccatctgaTaggtcgccttatcaatcctatgctgaatgccctctaaggagtt
acacatgctagcgcgtcggggtgg
-tcgccgcatgggacatgttggggtcccgtagctccggtcgatgtaggcacgcgttttgccgaggagatagatcatcagctgacctatagattcgtctgtc
+tcgccgcatgggacatgttggggtcccgtagctccggtcgatgtaggcacgcgAAAAgccgaggagatagatcatcagctgacctatagattcgtctgtc
gactctagcagagtggccagccac
aaggtttatattcagtcgaaatggacgggtcccgaacttgcacggacctaacgggactcgcctttcgtaataaccgaacctctaggccgccgcgagatca
cctcccgagccgagccacccgtca
type tileVariantID uint16 // 1-based
type tileLibRef struct {
- tag tagID
- variant tileVariantID
+ Tag tagID
+ Variant tileVariantID
}
type tileSeq map[string][]tileLibRef
maxtag := 0
for _, refs := range tseq {
for _, ref := range refs {
- if maxtag < int(ref.tag) {
- maxtag = int(ref.tag)
+ if maxtag < int(ref.Tag) {
+ maxtag = int(ref.Tag)
}
}
}
var kept, dropped int
for _, refs := range tseq {
for _, ref := range refs {
- if vars[int(ref.tag)] != 0 {
+ if vars[int(ref.Tag)] != 0 {
dropped++
} else {
kept++
}
- vars[int(ref.tag)] = ref.variant
+ vars[int(ref.Tag)] = ref.Variant
}
}
return vars, kept, dropped
skipOOO bool
taglib *tagLibrary
variant [][][blake2b.Size256]byte
+ refseqs map[string]map[string][]tileLibRef
// count [][]int
// seq map[[blake2b.Size]byte][]byte
variants int
// if non-nil, write out any tile variants added while tiling
encoder *gob.Encoder
- // if non-nil, call this func upon loading a genome
- onLoadGenome func(CompactGenome)
mtx sync.Mutex
}
for _, tv := range tvs {
// Assign a new variant ID (unique across all inputs)
// for each input variant.
- variantmap[tileLibRef{tag: tv.Tag, variant: tv.Variant}] = tilelib.getRef(tv.Tag, tv.Sequence).variant
+ variantmap[tileLibRef{Tag: tv.Tag, Variant: tv.Variant}] = tilelib.getRef(tv.Tag, tv.Sequence).Variant
}
return nil
}
-func (tilelib *tileLibrary) loadGenomes(genomes map[string][]tileVariantID, variantmap map[tileLibRef]tileVariantID, onLoadGenome func(CompactGenome)) error {
+func (tilelib *tileLibrary) loadCompactGenomes(cgs []CompactGenome, variantmap map[tileLibRef]tileVariantID, onLoadGenome func(CompactGenome)) error {
+ log.Debugf("loadCompactGenomes: %d", len(cgs))
var wg sync.WaitGroup
errs := make(chan error, 1)
- for name, variants := range genomes {
- name, variants := name, variants
+ for _, cg := range cgs {
+ name, variants := cg.Name, cg.Variants
wg.Add(1)
go func() {
defer wg.Done()
continue
}
tag := tagID(i / 2)
- newvariant, ok := variantmap[tileLibRef{tag: tag, variant: variant}]
+ newvariant, ok := variantmap[tileLibRef{Tag: tag, Variant: variant}]
if !ok {
err := fmt.Errorf("oops: genome %q has variant %d for tag %d, but that variant was not in its library", name, variant, tag)
select {
}
return
}
+ log.Tracef("loadCompactGenomes: cg %s tag %d variant %d => %d", name, tag, variant, newvariant)
variants[i] = newvariant
}
+ if onLoadGenome != nil {
+ onLoadGenome(cg)
+ }
if tilelib.encoder != nil {
- for name, variants := range genomes {
- cg := CompactGenome{
- Name: name,
- Variants: variants,
- }
- if onLoadGenome != nil {
- onLoadGenome(cg)
- }
- err := tilelib.encoder.Encode(LibraryEntry{
- CompactGenomes: []CompactGenome{cg},
- })
- if err != nil {
- select {
- case errs <- err:
- default:
- }
- return
+ err := tilelib.encoder.Encode(LibraryEntry{
+ CompactGenomes: []CompactGenome{cg},
+ })
+ if err != nil {
+ select {
+ case errs <- err:
+ default:
}
+ return
}
}
}()
return <-errs
}
+func (tilelib *tileLibrary) loadCompactSequences(cseqs []CompactSequence, variantmap map[tileLibRef]tileVariantID) error {
+ log.Debugf("loadCompactSequences: %d", len(cseqs))
+ for _, cseq := range cseqs {
+ for _, tseq := range cseq.TileSequences {
+ for i, libref := range tseq {
+ v, ok := variantmap[libref]
+ if !ok {
+ return fmt.Errorf("oops: CompactSequence %q has variant %d for tag %d, but that variant was not in its library", cseq.Name, libref.Variant, libref.Tag)
+ }
+ tseq[i].Variant = v
+ }
+ }
+ if tilelib.encoder != nil {
+ if err := tilelib.encoder.Encode(LibraryEntry{
+ CompactSequences: []CompactSequence{cseq},
+ }); err != nil {
+ return err
+ }
+ }
+ }
+ tilelib.mtx.Lock()
+ defer tilelib.mtx.Unlock()
+ if tilelib.refseqs == nil {
+ tilelib.refseqs = map[string]map[string][]tileLibRef{}
+ }
+ for _, cseq := range cseqs {
+ tilelib.refseqs[cseq.Name] = cseq.TileSequences
+ }
+ return nil
+}
+
+// Load library data from rdr. Tile variants might be renumbered in
+// the process; in that case, genomes variants will be renumbered to
+// match.
+//
+// If onLoadGenome is non-nil, call it on each CompactGenome entry.
func (tilelib *tileLibrary) LoadGob(ctx context.Context, rdr io.Reader, onLoadGenome func(CompactGenome)) error {
- genomes := map[string][]tileVariantID{}
+ cgs := []CompactGenome{}
+ cseqs := []CompactSequence{}
variantmap := map[tileLibRef]tileVariantID{}
err := DecodeLibrary(rdr, func(ent *LibraryEntry) error {
if ctx.Err() != nil {
}
if err := tilelib.loadTagSet(ent.TagSet); err != nil {
return err
- } else if err = tilelib.loadTileVariants(ent.TileVariants, variantmap); err != nil {
+ }
+ if err := tilelib.loadTileVariants(ent.TileVariants, variantmap); err != nil {
return err
- } else {
- for _, cg := range ent.CompactGenomes {
- genomes[cg.Name] = cg.Variants
- }
- return nil
}
+ cgs = append(cgs, ent.CompactGenomes...)
+ cseqs = append(cseqs, ent.CompactSequences...)
+ return nil
})
if err != nil {
return err
if ctx.Err() != nil {
return ctx.Err()
}
- err = tilelib.loadGenomes(genomes, variantmap, onLoadGenome)
+ err = tilelib.loadCompactGenomes(cgs, variantmap, onLoadGenome)
+ if err != nil {
+ return err
+ }
+ err = tilelib.loadCompactSequences(cseqs, variantmap)
if err != nil {
return err
}
log.Tracef("%s %s found[%d] == %#v", filelabel, job.label, i, f)
if last.taglen > 0 {
path = append(path, tilelib.getRef(last.tagid, job.fasta[last.pos:f.pos+f.taglen]))
+ } else {
+ f.pos = 0
}
last = f
}
if last.taglen > 0 {
path = append(path, tilelib.getRef(last.tagid, job.fasta[last.pos:]))
+ } else {
+ log.Warnf("%s %s no tags found", filelabel, job.label)
}
pathcopy := make([]tileLibRef, len(path))
if b != 'a' && b != 'c' && b != 'g' && b != 't' {
// return "tile not found" if seq has any
// no-calls
- return tileLibRef{tag: tag}
+ return tileLibRef{Tag: tag}
}
}
}
for i, varhash := range tilelib.variant[tag] {
if varhash == seqhash {
tilelib.mtx.Unlock()
- return tileLibRef{tag: tag, variant: tileVariantID(i + 1)}
+ return tileLibRef{Tag: tag, Variant: tileVariantID(i + 1)}
}
}
tilelib.variants++
}},
})
}
- return tileLibRef{tag: tag, variant: variant}
+ return tileLibRef{Tag: tag, Variant: variant}
}