Skip to content

Commit 703d303

Browse files
remove redundant filter for minimal hit length
1 parent ef33eb0 commit 703d303

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

polishTE

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -317,7 +317,7 @@ function get_blast_hit() {
317317
blastn -query $2 -db $3 -outfmt "6 qseqid sseqid pident length mismatch qstart qend sstart send sstrand" -evalue $4 | awk -v "ml=$5" '{OFS="\t"; if ($4 > ml) {print $0}}' >> $1.blast.o
318318

319319
# filter blastn results
320-
awk -v cons_len="$8" '{OFS="\t"; if ($1~/^#/) {} else {if ($3 > 80 && $7-$6 > 30 && $7-$6<3*cons_len) {print}}}' < $1.blast.o > $1.blast.tmp.o
320+
awk -v cons_len="$8" '{OFS="\t"; if ($1~/^#/) {} else {if ($3 > 80 && $7-$6<3*cons_len) {print}}}' < $1.blast.o > $1.blast.tmp.o
321321

322322
# parse blast result into a bed file
323323
awk '{OFS="\t"; if ($10~/plus/) {print $2, $8, $9, $1, $3, "+"} else {print $2, $9, $8, $1, $3, "-"}}' < $1.blast.tmp.o > $1.blast.bed

0 commit comments

Comments
 (0)