Sponsored Content
Top Forums UNIX for Dummies Questions & Answers Redirecting stdout inside a loop Post 302912750 by hmortens on Monday 11th of August 2014 04:34:22 PM
Old 08-11-2014
Quote:
You have to organize it before it gets dumped into the same giant pile, not after, so we need to alter their output, not edit the great giant file, I think.
is it possible to organize it real time? like as it is being written to file?

Here is the code for CutAdapt, where the /Processed reads/ and /Trimmed reads/ come from:
Code:
#run CutAdapt to trim Wafergen adapter 
#/usr/local/bin/python2.7 /illumina/runs/Runs/cutadapt-1.2.1/bin/cutadapt -a AGATCGGAAGAGCACACGTCT -o $i/${y#*/Sample_*}_cutadapt_AdapterRemoved.fastq.gz $infile

Here is the call to sRNAbench where the other bits come from.
Code:
#run sRNAbench w/o adapter trimming
#java -jar /illumina/runs/Runs/miRanalyzer/sRNAbenchDB/sRNAbench.jar input=$infile2 output=${y#*/Sample_*}.sRNAbench.$date_formatted dbPath=/illumina/runs/Runs/miRanalyzer/sRNAbenchDB/ species=$speciesBuild microRNA=$speciesAbbrev p=14 noMM=0 alignType=v minRC=2

I was not specific about this before because I thought it would be the same pattern search, but I need to pull lines containing "No. input reads:", "No. reads in analysis:", "mapped...reads to genomes(s)", and "Detected:" The stdout for that portion looks like this:
Code:
           START WITH THE PRE-PROCESSING OF THE READS             

               Reading file: /illumina/runs/Runs/140513_H207_0249_AD2B9LACXX/Unaligned_Lane1/Project_DefaultProject/Sample_GRC270_DEHP2_67C/GRC270_DEHP2_67C_cutadapt_AdapterRemoved.fastq.gz

             Result of pre-processing

               No. input reads: 2139267
               No. cleaned input reads (adapter found and trimmed): 0
               No. input reads where the adapter was not found: 2139267
               No. length filtered input reads (min. Length): 192482
               No. length filtered input reads (max. Length): 0
               No. reads in analysis: 1570818
               No. unique reads in analysis: 78295
               Max. read length in analysis: 51
               Filtered reads (low quality or low read count): 375967
               Max. read length in input file: 51

           FINISHED PRE-PROCESSING

I have already shown you the summary.txt file.
The call to that program is
#run QC with FastQC
#/illumina/runs/Runs/FastQC/fastqc --outdir=$LOCATION --quiet $infile2
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

redirecting STDOUT & STDERR

In bash, I need to send the STDOUT and STDERR from a command to one file, and then just STDERR to another file. Doing one or the other using redirects is easy, but trying to do both at once is a bit tricky. Anyone have any ideas? (9 Replies)
Discussion started by: jshinaman
9 Replies

2. Shell Programming and Scripting

implicitly redirecting stdout to a file

Is there a way to redirect all stdout to a file implicitly - like defining stdout=/home/me/process.log - so that all "echo" commands in several scripts/subscripts are written to that file; instead of having to edit all scripts to redirect the "echo" (e.g. echo 'This is a test ' >>... (1 Reply)
Discussion started by: ALTRUNVRSOFLN
1 Replies

3. Shell Programming and Scripting

Redirecting part of output to stdout

Hi, I am trying to execute a command like this: find ./ -name "*.gz" -exec sh -c 'zcat {} | awk -f parse.awk' \; >> output If I want to print the filename, i generally use the -print argument to the find command but when I am redirecting the output to a file, how can I print just the... (2 Replies)
Discussion started by: Legend986
2 Replies

4. UNIX for Dummies Questions & Answers

Redirecting several outputs to /dev/stdout

I have an executable that, depending on its input, outputs to either one file or several. It usually prints nothing on screen. The usual way to call this program is to specify an input and output filenames, like this: ./executable.exe -i inputfile -o outputfileIt will then try to use the output... (1 Reply)
Discussion started by: aplaydoc
1 Replies

5. Shell Programming and Scripting

redirecting to stdout in betwen command

can anyone help me in making singleline command for Capital Letters are folders ,small letter are files X,Y,Z are subfolders of A as shown below A - X,Y,Z Folder X has three files a.txt,b.txt,c.txt similarly Y,Z. as shown below X- a.txt,b.txt,c.txt Y- a.txt,b.txt,c.txt Z-... (4 Replies)
Discussion started by: phoenix_nebula
4 Replies

6. Shell Programming and Scripting

Redirecting stdout to variable while printing it

Hi everybody, I am trying to do the thing you see in the title, and I can't simply do a=$(svn up) echo $a because the program (svn) gives output on lots of lines and in the variable the output is stored on only one line (resulting in a horribly formatted text). Any tips? Thanks,... (2 Replies)
Discussion started by: ocirne94
2 Replies

7. Shell Programming and Scripting

Redirecting stdout on background task

Hello, I have a script (videostream.sh) which invokes the GStreamer command-line tool gst-launch with all the correct command line parameters. When I invoke this program, I add the '&' character at the end to make it a background task, so that my script can complete and exit, i.e. gst-launch... (1 Reply)
Discussion started by: salukibob
1 Replies

8. Shell Programming and Scripting

Redirecting stdout problem

I have a simple bash script that prints sth every 5 seconds. What I do is the following. I redirect the output of the script to a file, tail the file and see that it works and then from another console I delete the file where the output is redirected to. Even though I have deleted the file, the... (2 Replies)
Discussion started by: igurov
2 Replies

9. Shell Programming and Scripting

Redirecting stdout continously to a file

I have a C program that continously outputs info to stdout. The problem is that I am redirecting the stdout and stderr to a file and stdout is written at the end of the problem rather than continously to the file. This could be a problem if for example the program is killed and the stdout output is... (3 Replies)
Discussion started by: igurov
3 Replies

10. Shell Programming and Scripting

Redirecting stdion, stdout within an AT command

Hello, I'm strugling with some redirecting and all help is apreciated. The following program is working as expected, but the result of the AT command doesn't go to any file. Thanks in advance for the help. #!/bin/bash modem=/dev/ttyUSB1 file=/root/imsi.txt # print error to stderr and exit... (4 Replies)
Discussion started by: cleitao
4 Replies
fst-infl(1)							     fst-infl							       fst-infl(1)

NAME
fst-infl - morphological analysers SYNOPSIS
fst-infl [ options ] file [ input-file [ output-file ] ] fst-infl2 [ options ] file [ input-file [ output-file ] ] fst-infl3 [ options ] file [ input-file [ output-file ] ] OPTIONS
-t file Read an alternative transducer from file and use it if the main transducer fails to find an analysis. By iterating this option, a cascade of transducers may be tried to find an analysis. -b Print surface and analysis symbols. (fst-infl2 only) -n Print multi-character symbols without the enclosing angle brackets. (fst-infl only) -d The analyses are symbolically disambiguated by returning only analyses with a minimal number of morphemes. This option requires that morpheme boundaries are marked with the tag <X>. If no <X> tag is found in the analysis string, then the program (basically) counts the number of multi-character symbols consisting entirely of upper-case characters and uses this count for disambiguation. The lat- ter heuristic was developed for the German SMOR morphology. (This option is only available with fst-infl2 and fst-infl3.) -e n If no regular analysis is found, do robust matching and print analyses with up to n edit errors. The set of edit operations cur- rently includes replacement, insertion and deletion. Each operation has currently a fixed error weight of 1. (fst-infl2 only) -% f Disambiguates the analyses statistically and prints the most likely analyses with at least f % of the total probability mass of the analyses. The transducer weights are read from a file obtained by appending .prob to the name of the transducer file. The weight files are created with fst-train. (fst-infl2 only) -p Print the probability of each analysis. (fst-infl2 only) -c use this option if the transducer was compiled on a computer with a different endianness. If you have a transducer which was com- piled on a Sparc computer and you want to use it on a Pentium, you need to use this option. (fst-infl2 only) -q Suppress status messages. -h Print usage information. DESCRIPTION
fst-infl is a morphological analyser. The first argument is the name of a file which was generated by fst-compiler. The second argument is the name of the input file. The third argument is the output file. If the third argument is missing, output is directed to stdout. If the second argument is missing, as well, input is read from stdin. fst-infl2 is similar to fst-infl but needs a transducer in compact format (see the man pages for fst-compiler and fst-compact). fst-infl2 is implemented differently from fst-infl and usually much faster. fst-infl3 is also similar to fst-infl but needs a transducer in lowmem format (see the man pages for fst-compiler and fst-lowmem). fst- infl3 accesses the transducer on disc rather than reading it into memory. It starts very fast and needs very little memory, but is slower than fst-infl2. fst-infl reads the transducer which is stored in the argument file. Then it reads the input file line by line. Each line is analysed with the transducer and all resulting analyses are printed (see also the man pages for fst-mor). BUGS
No bugs are known so far. SEE ALSO
fst-compiler, fst-mor AUTHOR
Helmut Schmid, Institute for Computational Linguistics, University of Stuttgart, Email: schmid@ims.uni-stuttgart.de, This software is available under the GNU Public License. November 2004 fst-infl(1)
All times are GMT -4. The time now is 04:24 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy