From 9791287ca19ccc0a2192923a2d4269fe1d7c913b Mon Sep 17 00:00:00 2001 From: Wenfei Xian <2316469036@qq.com> Date: Wed, 18 Oct 2023 10:14:32 +0200 Subject: [PATCH 1/4] Update seqtk.c --- seqtk.c | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/seqtk.c b/seqtk.c index 962b9a8..eebcb84 100644 --- a/seqtk.c +++ b/seqtk.c @@ -2010,7 +2010,8 @@ int stk_telo(int argc, char *argv[]) else if (max - score > max_drop) break; } if (max >= min_score) { - if (!show_profile) printf("%s\t0\t%ld\t%ld\n", seq->name.s, max_i + 1, seq->seq.l); + if (!show_profile) printf(">%s\t0\t%ld\t%ld\n", seq->name.s, max_i + 1, seq->seq.l); + printf("%s\n", seq->seq.s); sum_telo += max_i + 1; st = max_i + 1; } @@ -2028,7 +2029,8 @@ int stk_telo(int argc, char *argv[]) else if (max - score > max_drop) break; } if (max >= min_score) { - if (!show_profile) printf("%s\t%ld\t%ld\t%ld\n", seq->name.s, max_i, seq->seq.l, seq->seq.l); + if (!show_profile) printf(">%s\t%ld\t%ld\t%ld\n", seq->name.s, max_i, seq->seq.l, seq->seq.l); + printf("%s\n", seq->seq.s); sum_telo += seq->seq.l - max_i; } } From 88661cfc77e58abf2c1cf2303279eb3dee56d5e5 Mon Sep 17 00:00:00 2001 From: Wenfei Xian <2316469036@qq.com> Date: Wed, 18 Oct 2023 10:16:30 +0200 Subject: [PATCH 2/4] Update README.md --- README.md | 1 + 1 file changed, 1 insertion(+) diff --git a/README.md b/README.md index 417a529..a7ae563 100644 --- a/README.md +++ b/README.md @@ -62,3 +62,4 @@ Seqtk Examples * Find telomere (TTAGGG)n repeats: seqtk telo seq.fa > telo.bed 2> telo.count + In this folk, the origin sequences will be output From 9c141fb0570d22fb3c5b18f168ded6b96cea1217 Mon Sep 17 00:00:00 2001 From: Wenfei Xian <2316469036@qq.com> Date: Wed, 18 Oct 2023 10:17:02 +0200 Subject: [PATCH 3/4] Update README.md --- README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/README.md b/README.md index a7ae563..7578ffa 100644 --- a/README.md +++ b/README.md @@ -61,5 +61,5 @@ Seqtk Examples * Find telomere (TTAGGG)n repeats: - seqtk telo seq.fa > telo.bed 2> telo.count + seqtk telo seq.fa > telo.fasta In this folk, the origin sequences will be output From 96eb70eabc29bc1a221b2c07c4b265f84a449fb0 Mon Sep 17 00:00:00 2001 From: Wenfei Xian <2316469036@qq.com> Date: Fri, 20 Oct 2023 19:35:20 +0200 Subject: [PATCH 4/4] Add files via upload --- seqtk.c | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/seqtk.c b/seqtk.c index eebcb84..d99af8d 100644 --- a/seqtk.c +++ b/seqtk.c @@ -2012,6 +2012,8 @@ int stk_telo(int argc, char *argv[]) if (max >= min_score) { if (!show_profile) printf(">%s\t0\t%ld\t%ld\n", seq->name.s, max_i + 1, seq->seq.l); printf("%s\n", seq->seq.s); + if (!show_profile) printf(">%s_remove\t%ld\t%ld\t%ld\n", seq->name.s, max_i, seq->seq.l, seq->seq.l); + printf("%s\n", seq->seq.s + max_i + 1); sum_telo += max_i + 1; st = max_i + 1; } @@ -2030,7 +2032,12 @@ int stk_telo(int argc, char *argv[]) } if (max >= min_score) { if (!show_profile) printf(">%s\t%ld\t%ld\t%ld\n", seq->name.s, max_i, seq->seq.l, seq->seq.l); - printf("%s\n", seq->seq.s); + printf("%s\n", seq->seq.s); + if (!show_profile) printf(">%s_remove\t%ld\t%ld\t%ld\n", seq->name.s, max_i, seq->seq.l, seq->seq.l); + char temp_char = seq->seq.s[max_i]; + seq->seq.s[max_i] = '\0'; + printf("%s\n", seq->seq.s); + seq->seq.s[max_i] = temp_char; sum_telo += seq->seq.l - max_i; } }