subset_ont_reads_by_max_sequencing_time.sh

#!/bin/bash
hours=$1;
files=${@:2};
if [[ -z ${hours} || -z ${files} ]]; then
    echo "Usage: reads_for_seqtime <hours> <files>";
    return;
fi;
seqkit fx2tab ${files} \
    | sort -k 8,8 \
    | awk '{ split($8, dt, "[=T:Z-]"); ts=mktime(dt[2]" "dt[3]" "dt[4]" "dt[5]" "dt[6]" "dt[7]); } NR==1 { start=ts } { seqtime=ts-start; print(seqtime"\t"$0) }' \
    | awk '( $1 <= '${hours}'*3600 )' \
    | cut -d$'\t' -f 2- \
    | seqkit tab2fx