私はfastaファイルのコンティグ内のヌクレオチドを数えています。私のファイルには、私のコマンドが同じ行のawkプリント出力
awk '/^>/ {if (seqlen){print seqlen}; print ;seqlen=0; } { seqlen += length($0)}END{print seqlen}'
出力は以下のように、同じ行に出力を得るためにどのように
>1
23
>2
21
ようになりましたです
>1
ATACCTACTA
ATTTACGTCA
GTA
>2
ATATTCGTAT
GTCTCGATCT
A
>3
などのように見えます
>1 23
>2 21
感謝。私は正しい例で私の質問に言及しなかった。私は各行だけを数えたいとは思わないが、文字 ">"の間の行は、私の質問を編集している。私が使用しているawkコマンドは正しい動作をしますが、今は入力と同じ行に出力します。あなたのソリューションでは、それは行だけを数えますが、行のすべての行を>で数えなければなりません。 – Gravel
私は私の答えを変えました。今すぐ動作するはずです:) – JFS31
@Gravel:今すぐ確認してください。私の答えが更新されました。 – RavinderSingh13