Skip to content

Commit

Permalink
Apply annotation to all split RS
Browse files Browse the repository at this point in the history
  • Loading branch information
tcezard committed Apr 15, 2024
1 parent ea48d6f commit 7e515c7
Showing 1 changed file with 1 addition and 1 deletion.
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,7 @@ do
if [[ ${INPUT} == *.vcf.gz ]]
then
# There are sometime multiple rs (separated by ;) in one line that needs to be split across multiple lines
zcat "${INPUT}" | grep -v '^#' | awk -v annotation="${ASSEMBLY}-${SC_NAME}-${TYPE}" '{gsub(";","\n",$3); print $3" "annotation}' > ${OUTPUT}
zcat "${INPUT}" | grep -v '^#' | awk '{gsub(";","\n",$3); print $3}' | awk -v annotation="${ASSEMBLY}-${SC_NAME}-${TYPE}" '{print $0" "annotation}' > ${OUTPUT}
elif [[ ${INPUT} == *_unmapped_ids.txt.gz ]]
then
SC_NAME=$(basename $(dirname ${INPUT}));
Expand Down

0 comments on commit 7e515c7

Please sign in to comment.