新しい出力csvファイルに "samtools"を使用して複数の入力ファイルの結果を見つけるbashスクリプト

新しい出力csvファイルに "samtools"を使用して複数の入力ファイルの結果を見つけるbashスクリプト

生物情報学に関する質問があります。複数のbamファイルがあります。各bamファイルで、IDを持つ各bamファイルに対してマップされた読み取り数を見つける必要があります。これを行うには、次のコマンドがあります。

$ samtools view -c -f 1 -F 12 HG00173.chrom11.ILLUMINA.bwa.FIN.low_coverage.20111114.bam

これを行うには、1222456などの数字が出力されます。

私の要件は、2つの列と名前列1)Bam-idと列2)No_of_mapped_readsを含むcsvファイルからすべての入力bamファイルの出力を取得する必要があることです。

答え1

たぶん、次のようなものがあります。

parallel --tag samtools view -c -f 1 -F 12 ::: *.bam

これが望ましくない場合は、入力例と出力例を使用して質問を更新してください。

関連情報