cat d1              |   # 入力データ
tr -dc 'atgcATGC\n' |   # 塩基配列文字と改行以外を除く
tr 'atgc' 'ATGC'    |   # 大文字に統一
while read seq ; do
  revcomp=$(
    echo $seq |         # 元の塩基配列
    rev       |         # 逆順変換
    tr 'ATGC' 'TACG' )  # 相補関係に変換

  N_total=$(echo -n $seq              | wc -m)  # 全塩基長
  N_gc=$(   echo -n $seq | tr -d 'AT' | wc -m)  # GC配列長

  echo $seq $revcomp $N_total $N_gc
done