0

我有两个文件： SCR_location - 包含按升序排列的 SNP 位置信息。

 19687

 36075

 n...

modi_VCF - 一个包含每个 SNP 信息的 vcf 表。

  19687  G A     xxx:255,0,195 xxx:255,0,206

  20398  G C     0/0:0,255,255 0/0:0,208,255

  n...

我只想将具有匹配 SNP 位置的行保存到一个新文件中我编写了以下脚本，但它不起作用

cat SCR_location |while read SCR_l; do
    cat modi_VCF |while read line; do

            if  [ "$SCR_l" -eq "$line" ] ;
            then echo "$line" >> file
            else :
            fi

    done

done
4

1 回答 1

1

请您尝试 bash 解决方案：

declare -A seen
while read -r line; do
    seen[$line]=1
done < SCR_location

while read -r line; do
    read -ra ary <<< "$line"
    if [[ ${seen[${ary[0]}]} ]]; then
        echo "$line"
    fi
done < modi_VCF > file
  • 它首先遍历 SCR_location 并将 SNP 位置存储在关联数组中seen。
  • 接下来它扫描 modi_VCF 并且如果在关联数组中找到第一列值，则打印该行。

如果awk是你的选择，你也可以说：

awk 'NR==FNR {seen[$1]++; next} {if (seen[$1]) print}' SCR_location modi_VCF > file

[编辑]为了过滤掉未加工的行，只需将逻辑否定为：

awk 'NR==FNR {seen[$1]++; next} {if (!seen[$1]) print}' SCR_location modi_VCF > file_unmatched

上面的代码只输出不匹配的行。如果您想一次对匹配的行和不匹配的行进行排序，请尝试：

awk 'NR==FNR {seen[$1]++; next} {if (seen[$1]) {print >> "file_matched"} else {print >> "file_unmatched"} }' SCR_location modi_VCF

希望这可以帮助。

于 2019-11-19T23:30:28.223 回答