09-19-2012
Is 59 zero based as in start with the 60th character, and is 140 1-based, as in stop with the 140th character (I guss it could be zero-based stop before the cahracter at offset 140)?
The values are indented and divided into groups -- is that what the file looks like? You want the same indentation and division in the output?
ksh/bash can read the 6 fields into an array of variables, save the pairs of numbers in a secondary array, mathematically decompose the numbers into line, field and offset, and compose the output with the same white space. Some line's field values will set a state variable, and the right state begins the field cutting. Initial state says you capture and print from the ID line, second state is look for and print from the FT COILED lines or find SQ SEQUENCE, and SQ SEQUENCE says start looking for the subset of output fields. I assume the COILED lines are forward sequenced.
Last edited by DGPickett; 09-19-2012 at 02:23 PM..
10 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
Hi there,
I'm wanting to produce a shell script that will check through some file names and identify a skip in sequence (four digit seq num in file name).
I have played on the idea of havng a file that has a sorted list of file names which I can read line at a time and cut out the sequence... (1 Reply)
Discussion started by: nhatch
1 Replies
2. UNIX for Dummies Questions & Answers
A developer of mine has this requirement - I couldn't tell her quickly how to do it with UNIX commands or a quick script so she's writing a quick program to do it - but that got my curiousity up and thought I'd ask here for advice.
In a text file, there are some records (about half of them)... (4 Replies)
Discussion started by: LisaS
4 Replies
3. Shell Programming and Scripting
Hi,
I have a data as follow:
1 400
2 239
3 871
4 219
5 543
6 ...
7 ...
.. ...
.. ...
99 818
100 991
I want to replace the sequence number (column 1) that start from 150. The output should like this:
150 400
151 239 (3 Replies)
Discussion started by: nica
3 Replies
4. Shell Programming and Scripting
Hi,
I have a string like:
DBMS stats (Number Used | Percentage of total): 10 | 1.00%
I have a sed command to extract numbers from this string:
sed "s///g;s/^$/-1/;"
Output: 10100
However what I want the sed command to return is only the first number(regardless of its size) i.e.... (3 Replies)
Discussion started by: mccartj5
3 Replies
5. Shell Programming and Scripting
Dear Perl users,
I need your help to solve my problem below.
I want to print the sequence number without missing number within the range.
E.g. my sequence number :
1 2 3 4 5 6 7 8 11 12 13 14
my desired output:
1 -8 , 11-14
my code below but still problem with the result:
1 - 14
1 -... (2 Replies)
Discussion started by: mandai
2 Replies
6. Shell Programming and Scripting
Am using unix aix KSH...
I have the files called
MMRR0106.DAT
MMRR0206.DAT
MMRR0406.DAT
MMRR0506.DAT
MMRR0806.DAT
....
...
MMRR3006.DAT
MMRR0207.DAT
These files are in one dircetory /venky ?
I want the output like this ?
Missing files are :
MMRR0306.DAT
MMRR0606.DAT... (7 Replies)
Discussion started by: Venkatesh1
7 Replies
7. Shell Programming and Scripting
Hi Forum:
I have struggle with it and decide to use my eye ball to accomplish this.
Basically I am looking for sequence of date inside a file.
If one of the sequence repeat 2-3 time or skip once; it's still consider a match.
input text file:
Sep 6 A
Sep 6 A
Sep 10 A
Sep 7 B
Sep 8... (7 Replies)
Discussion started by: chirish
7 Replies
8. Shell Programming and Scripting
Need some help doing this ... with awk maybe
Input
0DF6
0DF7
0DF8
0DF9
0DFA
0DFB
0DFC
0DFD
0DFF
0E00
0E01
0E02
0E03
0E04
0E05
0E06 (11 Replies)
Discussion started by: greycells
11 Replies
9. Shell Programming and Scripting
Hi,
Need to add sequnce number to one of the csv file and please find below
actual requirement.
Input file
ABC,500
XXQ,700
ADF,400,
ART,200
Out put file should be
1,ABC,500
2,XXQ,700
3,ADF,400,
4,ART,200 (3 Replies)
Discussion started by: siva83
3 Replies
10. Shell Programming and Scripting
Hello, here I am posting my query again with modified data input files.
see my query is :
i have two input files file1 and file2.
file1 is smalldata.fasta
>gi|546671471|gb|AWWX01449637.1| Bubalus bubalis breed Mediterranean WGS:AWWX01:contig449636, whole genome shotgun sequence... (20 Replies)
Discussion started by: harpreetmanku04
20 Replies
LEARN ABOUT DEBIAN
bup-margin
bup-margin(1) General Commands Manual bup-margin(1)
NAME
bup-margin - figure out your deduplication safety margin
SYNOPSIS
bup margin [options...]
DESCRIPTION
bup margin iterates through all objects in your bup repository, calculating the largest number of prefix bits shared between any two
entries. This number, n, identifies the longest subset of SHA-1 you could use and still encounter a collision between your object ids.
For example, one system that was tested had a collection of 11 million objects (70 GB), and bup margin returned 45. That means a 46-bit
hash would be sufficient to avoid all collisions among that set of objects; each object in that repository could be uniquely identified by
its first 46 bits.
The number of bits needed seems to increase by about 1 or 2 for every doubling of the number of objects. Since SHA-1 hashes have 160 bits,
that leaves 115 bits of margin. Of course, because SHA-1 hashes are essentially random, it's theoretically possible to use many more bits
with far fewer objects.
If you're paranoid about the possibility of SHA-1 collisions, you can monitor your repository by running bup margin occasionally to see if
you're getting dangerously close to 160 bits.
OPTIONS
--predict
Guess the offset into each index file where a particular object will appear, and report the maximum deviation of the correct answer
from the guess. This is potentially useful for tuning an interpolation search algorithm.
--ignore-midx
don't use .midx files, use only .idx files. This is only really useful when used with --predict.
EXAMPLE
$ bup margin
Reading indexes: 100.00% (1612581/1612581), done.
40
40 matching prefix bits
1.94 bits per doubling
120 bits (61.86 doublings) remaining
4.19338e+18 times larger is possible
Everyone on earth could have 625878182 data sets
like yours, all in one repository, and we would
expect 1 object collision.
$ bup margin --predict
PackIdxList: using 1 index.
Reading indexes: 100.00% (1612581/1612581), done.
915 of 1612581 (0.057%)
SEE ALSO
bup-midx(1), bup-save(1)
BUP
Part of the bup(1) suite.
AUTHORS
Avery Pennarun <apenwarr@gmail.com>.
Bup unknown- bup-margin(1)