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)
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)
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)
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)
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)
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)
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)
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)
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
re-pcr
RE-PCR(1) General Commands Manual RE-PCR(1)NAME
re-PCR -- Find sequence tagged sites (STS) in DNA sequences
SYNOPSIS
re-PCR [-hV] -p hash-file [-g gaps] [-n mism] [-lq] [primer ...]
re-PCR [-hV] -P hash-file [-g gaps] [-n mism] [-l] [-m margin] [-O+|-] [-C batchcnt] [-o outfile] [-r+|-] [primers-file ...]
re-PCR [-hV] -s hash-file [-g gaps] [-n mism] [-lq] [-m margin] [-o outfile] [-r+|-] [left right lo[-hi] [...]]
re-PCR [-hV] -S hash-file [-g gaps] [-n mism] [-lq] [-m margin] [-O+|-] [-C batchcnt] [-o outfile] [-r+|-] [stsfile ...]
DESCRIPTION
Implements reverse searching (called Reverse e-PCR) to make it feasible to search the human genome sequence and other large genomes by per-
forming STS and primer searches.
OPTIONS -p=hash-file
Perform primer lookup using hash-file
-P=hash-file
Perform primer lookup using hash-file
-s=hash-file
Perform STS lookup using hash-file
-S=hash-file
Perform STS lookup using hash-file
-n=mism Set max allowed mismatches per primer for lookup
-g=gaps Set max allowed indels per primer for lookup
-m=margin Set variability for STS size for lookup
-l Use presize alignments (only if gaps>0)
-G Print alignments in comments
-d=min-max
Set default STS size
-r=+|- Enable/disable reverse STS lookup
-O=+|- Enable/disable syscall optimisation
-C=batchcnt
Set number of STSes per batch
-o=outfile
Set output file name
-q Quiet (no progress indicator)
EXAMPLE
famap -tN -b genome.famap org/chr_*.fa
fahash -b genome.hash -w 12 -f3 ${PWD}/genome.famap
re-PCR -s genome.hash -n1 -g1 ACTATTGATGATGA AGGTAGATGTTTTT 120-200
See famap(1) and fahash(1)SEE ALSO
/usr/share/doc/ncbi-epcr/README.txt
bioperl(1), e-pcr(1), famap(1) and fahash(1)AUTHORS
This manual page was written by Andreas Tille <tille@debian.org> for the Debian system (but may be used by others). Permission is granted
to copy, distribute and/or modify this document under the terms of the GNU General Public License, Version 2 any later version published by
the Free Software Foundation.
On Debian systems, the complete text of the GNU General Public License can be found in /usr/share/common-licenses/GPL.
April 2008 RE-PCR(1)