remove extra output column in pairs.tsv

closes #47
brentp · Apr 6, 2020 · c50a639 · c50a639
1 parent 3dbac36
commit c50a639
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 3 deletions.
diff --git a/CHANGES.md b/CHANGES.md
@@ -1,3 +1,7 @@
+v0.2.10
+=======
++ fix extra output column in pairs.tsv (#47)
+
 v0.2.9
 ======
 + support multi-sample GVCF and fix some GVCF cases (thanks @ameynert for implementing)

diff --git a/src/somalierpkg/relate.nim b/src/somalierpkg/relate.nim
@@ -122,7 +122,7 @@ proc add*(rt:var seq[relations], rel:relation, expected_relatedness:float) =
     rt.add(rel, expected_relatedness)
 
 
-const header = "$sample_a\t$sample_b\t$relatedness\t$ibs0\t$ibs2\t$hom_concordance\t$hets_a\t$hets_b\t$shared_hets\t$hom_alts_a\t$hom_alts_b\t$shared_hom_alts\t$ibs2\t$n\t$x_ibs0\t$x_ibs2\t$expected_relatedness"
+const header = "$sample_a\t$sample_b\t$relatedness\t$ibs0\t$ibs2\t$hom_concordance\t$hets_a\t$hets_b\t$shared_hets\t$hom_alts_a\t$hom_alts_b\t$shared_hom_alts\t$n\t$x_ibs0\t$x_ibs2\t$expected_relatedness"
 
 proc tsv(r:relation, expected_relatedness:float= -1.0): string =
   result = &"{r.sample_a}\t{r.sample_b}\t{r.rel:.3f}\t{r.ibs0}\t{r.ibs2}\t{r.hom_alt_concordance:.3f}\t{r.hets_a}\t{r.hets_b}\t{r.shared_hets}\t{r.hom_alts_a}\t{r.hom_alts_b}\t{r.shared_hom_alts}\t{r.n}\t{r.x_ibs0}\t{r.xibs2}\t{expected_relatedness}"
@@ -302,7 +302,7 @@ proc ab*(c:allele_count, min_depth:int): float {.inline.} =
     return 0
   result = c.nalt.float / (c.nalt + c.nref).float
 
-proc alts*(ab:float, ab_cutoff:float=0.04): int8 {.inline.} =
+proc alts*(ab:float, ab_cutoff:float=0.01): int8 {.inline.} =
   if ab < 0: return -1
   if ab < ab_cutoff: return 0
   if ab > (1 - ab_cutoff): return 2

diff --git a/src/somalierpkg/version.nim b/src/somalierpkg/version.nim
@@ -1 +1 @@
-let somalierVersion* = "0.2.9"
+let somalierVersion* = "0.2.10"
Original file line number	Diff line number	Diff line change
		@@ -1 +1 @@
		let somalierVersion* = "0.2.9"
		let somalierVersion* = "0.2.10"