Skip diffs on long ref seqs.
authorTom Clegg <tom@tomclegg.ca>
Thu, 15 Jul 2021 15:07:08 +0000 (11:07 -0400)
committerTom Clegg <tom@tomclegg.ca>
Thu, 15 Jul 2021 15:07:08 +0000 (11:07 -0400)
Arvados-DCO-1.1-Signed-off-by: Tom Clegg <tom@curii.com>

export.go

index 2ecf984ba1e094cff0be1b62ba6df73fa15cd3b5..47ff751afa4f6cc02a317ac1cfa0e34b51e4bb70 100644 (file)
--- a/export.go
+++ b/export.go
@@ -412,9 +412,11 @@ func (cmd *exporter) exportSeq(outw, bedw io.Writer, taglen int, seqname string,
                                                refstepend++
                                        }
                                        // (TODO: handle no-calls)
-                                       refstr := strings.ToUpper(string(refSequence))
-                                       genomestr := strings.ToUpper(string(genomeseq))
-                                       vars, _ = hgvs.Diff(refstr, genomestr, time.Second)
+                                       if len(refSequence) <= cmd.maxTileSize {
+                                               refstr := strings.ToUpper(string(refSequence))
+                                               genomestr := strings.ToUpper(string(genomeseq))
+                                               vars, _ = hgvs.Diff(refstr, genomestr, time.Second)
+                                       }
                                        diffs[glibref] = vars
                                }
                                for _, v := range vars {