その行がない行セットを印刷します。

これはawkのやや複雑なバージョンです。 Steeldriverのsedバージョンとのいくつかの違いは次のとおりです。

1. mm10記録が記録される順序についてはrn5何もしません。
2. rn5不足しているレコードを処理できます。
3. 一致しないレコードをstderr。
4. 追加コードは次のようになります:-)

以下から実行できます。

awk -f my_program.awk infile

パスワード：

# find and store a header
/^NM.*/ { header = $0; next }

# we found an mm10 line
header ~ /_rn5/ {

    # get the mm10 line that matches this rn5
    mm_match = header
    sub("_rn5", "_mm10", mm_match)

    # if we have a previous mm10, then print the pair
    if (mm_match in headers) {
        print header
        print
        print mm_match
        print headers[mm_match]

        delete headers[mm_match]
    } else {
        headers[header] = $0
    }
    next
}

# we found an mm10 line
header ~ /_mm10/ {

    # get the rn5 line that matches this mm10
    mm_match = header
    sub("_mm10", "_rn5", mm_match)

    # if we have a previous rn5, then print the pair
    if (mm_match in headers) {
        print mm_match
        print headers[mm_match]
        print header
        print

        delete headers[mm_match]
    } else {
        headers[header] = $0
    }
    next
}

また、このコードをファイルの末尾に追加して、一致しない行を次に出力することもできますstandard error。

# The END block is here just to output anything that was unmatched
END {
    # dump the unmatched to stderr
    for (header in headers) {
        print header > "/dev/stderr"
        print headers[header] > "/dev/stderr"
    }
}

以下から実行できます。

awk -f my_program.awk infile > outfile 2> unmatched

outfile要求された出力（標準出力を介して）をに出力し、残りの入力（標準エラーを介して）をに出力しますunmatched。さまざまな入出力リダイレクトの詳細については、次の章を参照してください。リダイレクトBashリファレンスマニュアルにあります。