Login or Register to Ask a Question and Join Our Community


Shell Programming and Scripting

Match and store numerical prefix to update files

Login to Discuss or Reply to this Discussion in Our Community

 
Thread Tools Search this Thread
Top Forums Shell Programming and Scripting Match and store numerical prefix to update files
# 1  
Old 02-02-2017
Match and store numerical prefix to update files
In the bash below the unique headers of each vcf.gz are stored in a text file with the same name. That is if 16-0000-file.vcf.gz was used the header text file would be 16-0000-file_header.txt.
There can be multiple vcf.gz in a directory, usually 3, that I need to fix the header in each file before further processing it. My question is how can I match each text file with its vcf.gz and pass the stored variables of each to the reheader code ?

In the below I strip off the unique numerical prefix 16-0000 from both the vcf.gz and text file, but am not sure how to match the two files

Code:
IAm=${0##*/}

InDir1='/home/cmccabe/Desktop/NGS/test'
InDir2='/home/cmccabe/Desktop/NGS/test'
OutDir='/home/cmccabe/Desktop/NGS/test'

cd "$InDir1"
for file1 in *.txt
do	# Grab file prefix.
	p=${file1%%_*}

	# Find matching file2.
	file2=$(printf '%s' "$InDir2/$p"_*.vcf.gz)
	if [ ! -f "$file2" ]
	then	printf '%s: No single file matching %s found.\n' "$IAm" \
		    "$file1" >&2
		continue
	fi

# store matches
    out=${file1##*/} && ${file2##*/}

vcf.gz in directory (file2)
Code:
16-0000-file1.vcf.gz
16-0001-file2.vcf.gz
16-0002-file3.vcf.gz

matching text file in directory (file1)
Code:
16-0000-file1_header.txt
16-0001-file2_header.txt
16-0002-file3_header.txt

So the contents of 16-0000.txt would be used to update 16-0000.vcf.gz using the code below.

reheader code
Code:
# edit the header
logfile=/home/cmccabe/Desktop/NGS/test/process.log
for f in /home/cmccabe/Desktop/NGS/test/*.vcf.gz ; do
     echo "Start vcf add header creation: $(date) - file: $f"
     bname=`basename $f`
     pref=${bname%%.vcf.gz}
      bcftools reheader -h $file1 $file2 > ${pref}_fixed.vcf.gz
     echo "End add header creation: $(date) - file: $f"
done >> "$logfile"

Last edited by cmccabe; 02-02-2017 at 01:40 PM.. Reason: fixed format, added details to bash
# 2  
Old 02-02-2017
From your example I don't understand what your script is doing that you don't want it to do or is not doing that you do want it to do.

Please show us:
  1. the output you are getting,
  2. the output you are hoping to get,
  3. the arguments that are currently being passed to reheader,
  4. the arguments you want to be passed to reheader, and
  5. ls -l output from the three directories involved.
This User Gave Thanks to Don Cragun For This Post:
cmccabe
# 3  
Old 02-02-2017
Code:
#!/bin/bash

# gunzip files
logfile=/home/cmccabe/Desktop/NGS/test/process.log
for f in /home/cmccabe/Desktop/NGS/test/*.vcf ; do
     echo "Start vcf.gz creation: $(date) - file: $f"
     bname=`basename $f`
     gzip $f
     echo "End vcf.gz creation: $(date) - file: $f"
done >> "$logfile"

# find undefined annotations
logfile=/home/cmccabe/Desktop/NGS/test/process.log
for f in /home/cmccabe/Desktop/NGS/test/*.vcf.gz ; do
     echo "Start vcf missing header creation: $(date) - file: $f"
     bname=`basename $f`
     pref=${bname%%.vcf.gz}
    bcftools view -h $f > /home/cmccabe/Desktop/NGS/test/${pref}_header.txt
     echo "End missing header creation: $(date) - file: $f"
done >> "$logfile"

# match files
IAm=${0##*/}

InDir1='/home/cmccabe/Desktop/NGS/test'
InDir2='/home/cmccabe/Desktop/NGS/test'
OutDir='/home/cmccabe/Desktop/NGS/test'

cd "$InDir1"
for file1 in *.txt
do	# Grab file prefix.
	p=${file1%%_*}

	# Find matching file2.
	file2=$(printf '%s' "$InDir2/$p"_*.vcf)
	if [ ! -f "$file2" ]
	then	printf '%s: No single file matching %s found.\n' "$IAm" \
		    "$file1" >&2
		continue
	fi
# store matches
    out=${file1##*/} && ${file2##*/}

# edit the header
logfile=/home/cmccabe/Desktop/NGS/test/process.log
for f in /home/cmccabe/Desktop/NGS/test/*.vcf.gz ; do
     echo "Start vcf edit header creation: $(date) - file: $f"
     bname=`basename $f`
     pref=${bname%%.vcf.gz}
      bcftools reheader -h $file1 $file2 > ${pref}_fixed.vcf.gz
     echo "End edit header creation: $(date) - file: $f"
done >> "$logfile"

Current output
syntax error: unexpected end of file, though the portion in bold creates the files that are need in the directory.

desired output
Code:
file1_header.txt is stored as $file1 matched with file1.vcf.gz stored as $file2
file2_header.txt is stored as $file1 matched with file2.vcf.gz stored as $file2f
file3_header.txt is stored as $file1 matched with file3.vcf.gz stored as $file2

Currently no arguments are being passed to reheader (in italics in the command)

desired arguments passed to reheader
Code:
$file1
$file2

The loop would use each argument in the below to create a new file updated with the fixed header:
Code:
bcftools reheader -h $file1 $file2 > ${pref}_fixed.vcf.gz

Code:
ls -l
total 12288
-rwxrwx---+ 1 cmccabe Domain Users 2031823 Jan 31 11:07 file1.vcf.gz
-rwxrwx---+ 1 cmccabe Domain Users    9312 Feb  2 13:17 file1_header.txt
-rwxrwx---+ 1 cmccabe Domain Users 2361873 Jan 31 11:07 file2.vcf.gz
-rwxrwx---+ 1 cmccabe Domain Users    9315 Feb  2 13:17 file2_header.txt
-rwxrwx---+ 1 cmccabe Domain Users 1816662 Jan 31 11:07 file3.vcf.gz
-rwxrwx---+ 1 cmccabe Domain Users    9313 Feb  2 13:17 file3_header.txt
-rwxrwx---+ 1 cmccabe Domain Users    1356 Feb  2 13:22 process.log
-rwxrwx---+ 1 cmccabe Domain Users    1278 Feb  2 13:17 process.log~

I hope this helps and thank you very much Smilie.
# 4  
Old 02-02-2017
Please explain in English what the following code is intended to do:
Code:
# store matches
    out=${file1##*/} && ${file2##*/}

Given that the variable out is never used after being defined by the above statement, are these two lines of code needed?

What happens if you insert the missing done in your script before or after the code discussed above (which should resolve the syntax error in your script)?
This User Gave Thanks to Don Cragun For This Post:
cmccabe
# 5  
Old 02-02-2017
The following code:

Code:
# store matches
    out=${file1##*/} && ${file2##*/}

is intended to match and store file1_header.txt as $file1 and match and store file1_header.vcf.gz as $file2.
These two variables are then passed to the reheader command to be processed. After the command executes the variables are reset using the remaining two files.
I will rerun the code omitting the out and adding the done. Thank you Smilie.

---------- Post updated at 04:53 PM ---------- Previous update was at 02:25 PM ----------

The portion of code in bold in the previous post is unchanged, but I updated the match portion of the code to:
adding the missing done did allow the script to execute.

Code:
# match files
IAm=${0##*/}

InDir1='/home/cmccabe/variants'
InDir2='/home/cmccabe/variants'

cd "$InDir1"
for file1 in *.txt
do    # Grab file prefix.
    p=${file1%%_*}

    # Find matching file2.
    file2=$(printf '%s' "$InDir2/$p"_*.vcf)
    if [ ! -f "$file2" ]
    then    printf '%s: No single file matching %s found.\n' "$IAm" \
            "$file1" >&2
        continue
    fi
        echo "file1is:"$file1
        echo "file2 is:"$file2
done

The output that I get in the terminal is:
Code:
No single file matching 16_0000_file-a_header.txt found.
No single file matching 16_0001_file-b_header.txt found.

Code:
ls -l /home/cmccabe/variants
total 92
-rw-rw-r-- 1 cmccabe cmccabe  8895 Feb  2 16:37 16_0000_file-a_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 29066 Jan 24 12:06 16_0000_file-a.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe  8895 Feb  2 16:37 16_0001_file-b_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 29066 Jan 24 12:06 16_0001_file-b.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe   860 Feb  2 16:37 process.log

desired output
Code:
$file1=16-0000.txt
$file2=16-0000.vcf.gz

$file1=16-0001.txt
$file2=16-0001.vcf.gz

Basically, on the first pass the $file1 variable is the .txt and $file2 variable is the matching vcf.gz. Those variable are passed to the reheader command it is executed. After it executes the variables are reset using the other files in the directory. Since the numerical prefix is always unique that is used to perform the match between files. There will also also be a match vcf.gz and .txt. Thank you Smilie.
Last edited by cmccabe; 02-02-2017 at 06:56 PM.. Reason: added details
# 6  
Old 02-03-2017
Quote:
Originally Posted by cmccabe
The following code:

Code:
# store matches
    out=${file1##*/} && ${file2##*/}

is intended to match and store file1_header.txt as $file1 and match and store file1_header.vcf.gz as $file2.
These two variables are then passed to the reheader command to be processed. After the command executes the variables are reset using the remaining two files.
I will rerun the code omitting the out and adding the done. Thank you Smilie.

---------- Post updated at 04:53 PM ---------- Previous update was at 02:25 PM ----------

The portion of code in bold in the previous post is unchanged, but I updated the match portion of the code to:
adding the missing done did allow the script to execute.

Code:
# match files
IAm=${0##*/}

InDir1='/home/cmccabe/variants'
InDir2='/home/cmccabe/variants'

cd "$InDir1"
for file1 in *.txt
do    # Grab file prefix.
    p=${file1%%_*}

    # Find matching file2.
    file2=$(printf '%s' "$InDir2/$p"_*.vcf)
    if [ ! -f "$file2" ]
    then    printf '%s: No single file matching %s found.\n' "$IAm" \
            "$file1" >&2
        continue
    fi
        echo "file1is:"$file1
        echo "file2 is:"$file2
done

The output that I get in the terminal is:
Code:
No single file matching 16_0000_file-a_header.txt found.
No single file matching 16_0001_file-b_header.txt found.

Code:
ls -l /home/cmccabe/variants
total 92
-rw-rw-r-- 1 cmccabe cmccabe  8895 Feb  2 16:37 16_0000_file-a_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 29066 Jan 24 12:06 16_0000_file-a.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe  8895 Feb  2 16:37 16_0001_file-b_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 29066 Jan 24 12:06 16_0001_file-b.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe   860 Feb  2 16:37 process.log

desired output
Code:
$file1=16-0000.txt
$file2=16-0000.vcf.gz

$file1=16-0001.txt
$file2=16-0001.vcf.gz

Basically, on the first pass the $file1 variable is the .txt and $file2 variable is the matching vcf.gz. Those variable are passed to the reheader command it is executed. After it executes the variables are reset using the other files in the directory. Since the numerical prefix is always unique that is used to perform the match between files. There will also also be a match vcf.gz and .txt. Thank you Smilie.
Given that there is nothing in your script that attempts to print the string $file1= nor the string $file2=, I have no idea how you would expect your script to produce that output???
Furthermore, since your ls output doesn't show any files with any of the filenames shown on the right hand side of those strings, I can't guess at how you would expect your script to create those names and find those names as existing files. (But of course I did ask for ls -l output from all three of your directories and you only showed output from one of those directories and didn't indicate which one of those three directories was the current working directory for the listing you provided, so maybe your names would make sense.)

The confusion between the the punctuation you use in the names of your files completely confuses any attempt to guess at what you want. Above we have .txt file filenames:
Code:
16_0000_file-a_header.txt
16_0001_file-b_header.txt
16-0000.txt
16-0001.txt

quoted from post #5 (with two showing underscores between the leading sequences of digits and two showing hyphens between the leading sequences of digits, while the names shown in post #1 were:
Code:
16-0000-file1_header.txt
16-0001-file2_header.txt
16-0002-file3_header.txt

using a completely different filename format.

Is it really that hard to give us an accurate description of the files you're working with and what you're trying to do?
This User Gave Thanks to Don Cragun For This Post:
cmccabe
# 7  
Old 02-03-2017
Some of the file names and locations were different because I was at the office then at home and that was a problem leading to typos. I apologize for that and am back in the office. The details below are the exact files in the directory.

There are always 3 vcf.gz files in the following format below in each directory:

Files in /home/cmccabe/Desktop/NGS/test
Code:
16-0000_File-A_variant_strandbias_readcount.vcf.gz
16-0002_File-A_variant_strandbias_readcount.vcf.gz
16-0005_File-A_variant_strandbias_readcount.vcf.gz

The numeric prefix of each is always 7 characters xx-xxxx and unique. The portion of code in the previous post uses each vcf.gz file to create a corresponding .txt file with only the unique numeric prefix followed by _header.

Files in /home/cmccabe/Desktop/NGS/test (if code works as expected)
Code:
16-0000_File-A_variant_strandbias_readcount.vcf.gz
16-0002_File-B_variant_strandbias_readcount.vcf.gz
16-0005_File-A_variant_strandbias_readcount.vcf.gz
16-0000_header.txt
16-0002_header.txt
16-0005_header.txt

I am trying to match the 16-0000.vcf.gz with the 16-0000.txt and read 16-0000.vcf.gz into $file1 and 16-0000.txt into $file2.

These two variables $file1 and $file2 would be passed to the reheader command, in bold,below which is part of a loop.

After the $file1 and $file2 are processed for 16-0000, the process is repeated for the remaining two files, 16-0002 and 16-0005.

reheader command
Code:
logfile=/home/cmccabe/Desktop/NGS/test/process.log
for f in /home/cmccabe/Desktop/NGS/test/*.vcf.gz ; do
     echo "Start vcf edit header creation: $(date) - file: $f"
     bname=`basename $f`
     pref=${bname%%.vcf.gz}
     bcftools reheader -h $file1 $file2 > ${pref}_fixed.vcf.gz
     echo "End edit header creation: $(date) - file: $f"
done >> "$logfile"

There is only one directory that contains the files and that is /home/cmccabe/Desktop/NGS/test. The numerical prefixes are not being stripped off and used to perform the match between files nor are the variables being passed to reheader. Thank you Smilie.

Code:
ls -l /home/cmccabe/Desktop/NGS/test
total 6116
-rw-rw-r-- 1 cmccabe cmccabe    9350 Feb  3 07:08 16-0000_File-A_variant_strandbias_readcount_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 2031861 Jan 31 11:07 16-0000_File-A_variant_strandbias_readcount.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe    9353 Feb  3 07:08 16-0002_File-B_variant_strandbias_readcount_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 2361911 Jan 31 11:07 16-0002_File-B_variant_strandbias_readcount.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe    9351 Feb  3 07:08 16-0005_File-C_variant_strandbias_readcount_header.txt
-rw-rw-r-- 1 cmccabe cmccabe 1816700 Jan 31 11:07 16-0005_File-C_variant_strandbias_readcount.vcf.gz
-rw-rw-r-- 1 cmccabe cmccabe    2622 Feb  3 07:08 process.log
-rw-rw-r-- 1 cmccabe cmccabe    1278 Feb  2 13:17 process.log~

Does this help? Thank you very much Smilie.
Last edited by cmccabe; 02-07-2017 at 05:19 PM..
Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Bash to update file on prefix match in two directories

I am trying to use bash to loop through a directory /path/to/data using a prefix match from /path/to/file. That match is obtained and works using the code below (in green)... what I can not seem to do is populate or update the corresponding prefix_file.txt in /path/to/data with the values in each... (3 Replies)
Discussion started by: cmccabe
3 Replies

2. Shell Programming and Scripting

awk move select fields to match file prefix in two directories

In the awk below I am trying to use the file1 as a match to file2. In file2 the contents of $5,&6,and $7 (always tab-delimited) and are copied to the output under the header Quality metrics. The below executes but the output is empty. I have added comments to help and show my thinking. Thank you... (0 Replies)
Discussion started by: cmccabe
0 Replies

3. Shell Programming and Scripting

awk to update file with numerical difference if condition is met

In the file1 below if $9 and $12 are . (dot) then the value in $8 of file1 is used as a key (exact match) to lookup in each $2 of file2, when a match is found then the value of $4 in file1 is used to look for a range match within +/- 50 using the values in $4 and after in file2. The number of... (9 Replies)
Discussion started by: cmccabe
9 Replies

4. Shell Programming and Scripting

Bash to add portion of text to files in directory using numerical match

In the below bash I am trying to rename eachof the 3 text files in /home/cmccabe/Desktop/percent by matching the numerical portion of each file to lines 3,4, or 5 in /home/cmccabe/Desktop/analysis.txt. There will always be a match between the files. When a match is found each text file in... (2 Replies)
Discussion started by: cmccabe
2 Replies

5. Shell Programming and Scripting

awk to update file based on partial match in field1 and exact match in field2

I am trying to create a cronjob that will run on startup that will look at a list.txt file to see if there is a later version of a database using database.txt as the source. The matching lines are written to output. $1 in database.txt will be in list.txt as a partial match. $2 of database.txt... (2 Replies)
Discussion started by: cmccabe
2 Replies

6. Shell Programming and Scripting

Bash to match and store line as variable

The bash below loops through a specific directory dir and finds and writes the oldest folder to a variable called $filename. #!/bin/bash # oldest folder stored as variable for analysis, version log created, and quality indicators matched to run dir=/home/cmccabe/Desktop/NGS/test find... (2 Replies)
Discussion started by: cmccabe
2 Replies

7. Shell Programming and Scripting

Extract Uniq prefix from a start and end prefix

Dear All, assume i have a file with content: <Start>6000</Start> <Stop>7599</Stop> the output is: 6000 7000 7100 7200 7300 7400 7599 how should we use any awk, sed, perl can do this task, means to extract the uniq prefixes from the start and stop prefix. Thanks Jimmy (3 Replies)
Discussion started by: jimmy_y
3 Replies

8. Shell Programming and Scripting

Match columns from two csv files and update field in one of the csv file

Hi, I have a file of csv data, which looks like this: file1: 1AA,LGV_PONCEY_LES_ATHEE,1,\N,1,00020460E1,0,\N,\N,\N,\N,2,00.22335321,0.00466628 2BB,LES_POUGES_ASF,\N,200,200,00006298G1,0,\N,\N,\N,\N,1,00.30887539,0.00050312... (10 Replies)
Discussion started by: djoseph
10 Replies

9. Shell Programming and Scripting

Match Pattern and store next value into array

Hi, I am trying to write a script which parses a log file and will eventually put the values in an array so that I can perform some math on it. In this file I am only interested in the last 200 lines so here is the command I use to display the contents in a manageable manner. tail -200... (3 Replies)
Discussion started by: Filter500
3 Replies

10. Programming

Fuzzy Match Logic for Numerical Values

I have searched the internet (including these forums) and perhaps I'm not using the right wording. What I'm looking for is a function (preferably C) that analyzes the similitude of two numerical or near-numerical values, and returns either a true/false (match/nomatch) or a return code that... (4 Replies)
Discussion started by: marcus121
4 Replies
Login or Register to Ask a Question