Diff command on two Fastq.gz files Post: 302876927

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

diff 2 files; output diff's to 3rd file

Hello, I want to compare two files. All records in file 2 that are not in file 1 should be output to file 3. For example: file 1 123 1234 123456 file 2 123 2345 23456 file 3 should have 2345 23456 I have looked at diff, bdiff, cmp, comm, diff3 without any luck!

2. Shell Programming and Scripting

Find duplicates from multuple files with 2 diff types of files

I need to compare 2 diff type of files and find out the duplicate after comparing each types of files: Type 1 file name is like: file1.abc (the extension abc could any 3 characters but I can narrow it down or hardcode for 10/15 combinations). The other file is file1.bcd01abc (the extension...

3. AIX

diff command

hello i've two files. how i get the diff between the two files to new file. thanks best regards ariec

4. Shell Programming and Scripting

diff command help

Hi all diff file1 file 2 command will give us op of diff between two file. But it aslo give its position and sign "<" or ">". I dont want position and sign in op. Only diff of content should be come as op. Kindly help me for this. Regards Jaydeep

5. UNIX for Dummies Questions & Answers

Diff command of two files

Hi, I use the diff command to compare two files and append this output to a file. I would like to now not only produce the differences but be able to output the total number of changes made, the number of new files added and the number of files deleted, is there I can do this using the diff...

6. Shell Programming and Scripting

Combining 3 fastq files

Hello, I am working with next-gen short-read sequence data, which we receive in 3 fastq files. These are arranged in 4-line groups for each read: line1: read identifier, beginning, e.g., "@HWI-ST1342..." line2: DNA sequence, for files 1 and 2, 101 characters, for file 3, 7 chars. line3: "+"...

7. Shell Programming and Scripting

Extract length wise sequences from fastq file

I have a fastq file from small RNA sequencing with sequence lengths between 15 - 30. I wanted to filter sequence lengths between 21-25 and write to another fastq file. how can i do that?

8. Shell Programming and Scripting

One-way diff command?

Hello, I am trying to find the different files between multiple directories in Linux, here is a small assumption of what is inside the directories dir1 dir2 dir3 1.txt 1.txt 1.txt 2.txt 3.txt 3.txt 5.txt 4.txt 5.txt 6.txt 7.txt 8.txt I am using the following...

9. Shell Programming and Scripting

Diff 3 files, but diff only their 2nd column

Guys i have 3 files, but i want to compare and diff only the 2nd column path=`/home/whois/doms` for i in `cat domain.tx` do whois $i| sed -n '/Registry Registrant ID:/,/Registrant Email:/p' > $path/$i.registrant whois $i| sed -n '/Registry Admin ID:/,/Admin Email:/p' > $path/$i.admin...

10. UNIX for Beginners Questions & Answers

Comparing fastq files and outputting common records

I have two files: File_1: @M04961:22:000000000-B5VGJ:1:1101:9280:7106 1:N:0:86 GGCATGAAAACATACAAACCGTCTTTCCAGAAATTGTTCCAAGTATCGGCAACAGCTTTATCAATACCATGAAAAATATCAACCACACCAGAAGCAGCAT + GGGGGGGGGGGGGGGGGCCGGGGGF,EDFFGEDFG,@DGGCGGEGGG7DCGGGF68CGFFFGGGG@CGDGFFDFEFEFF:30CGAFFDFEFF8CAF;;8F ...

LEARN ABOUT DEBIAN

lastal

LASTAL(1)							   User Commands							 LASTAL(1)

NAME

       lastal - genome-scale comparison of biological sequences

SYNOPSIS

       lastal [options] lastdb-name fasta-sequence-file(s)

DESCRIPTION

       Find local sequence alignments.

       Score  options  (default  settings): -r: match score   (DNA: 1, protein: blosum62, 0<Q<5:  6) -q: mismatch cost (DNA: 1, protein: blosum62,
       0<Q<5: 18) -p: file for residue pair scores -a: gap existence cost (DNA: 7, protein: 11, 0<Q<5: 21) -b: gap extension cost  (DNA:  1,  pro-
       tein:   2,  0<Q<5:   9)	-c:  unaligned	residue  pair cost (100000) -F: frameshift cost (off) -x: maximum score drop for gapped alignments
       (max[y, a+b*20]) -y: maximum score drop for gapless alignments (t*10) -z: maximum score drop for final gapped alignments  (x)  -d:  minimum
       score for gapless alignments (e*3/5) -e: minimum score for gapped alignments (DNA: 40, protein: 100, 0<Q<5: 180)

       Cosmetic  options  (default  settings): -h: show all options and their default settings -v: be verbose: write messages about what lastal is
       doing -o: output file -f: output format: 0=tabular, 1=maf (1)

       Miscellaneous options (default settings): -s: strand: 0=reverse, 1=forward, 2=both (2 for DNA, 1 for protein) -m: maximum multiplicity  for
       initial	matches (10) -l: minimum length for initial matches (1) -n: maximum number of gapless alignments per query position (infinity) -k:
       step-size along the query sequence (1) -i: query batch size (1 MiB if Q>0, else 16 MiB if j=0, else 128	MiB)  -u:  mask  lowercase  during
       extensions: 0=never, 1=gapless,

	      2=gapless+gapped but not final, 3=always (2 if lastdb -c and Q<5, else 0)

       -w:  supress  repeats  inside exact matches, offset by this distance or less (1000) -G: genetic code file -t: 'temperature' for calculating
       probabilities (1/lambda) -g: 'gamma' parameter for gamma-centroid and LAMA (1) -j: output  type:  0=match  counts,  1=gapless,  2=redundant
       gapped, 3=gapped,

	      4=column ambiguity estimates, 5=gamma-centroid, 6=LAMA (3)

       -Q: input format: 0=fasta, 1=fastq-sanger, 2=fastq-solexa, 3=fastq-illumina,

	      4=prb, 5=PSSM (0)

REPORTING BUGS

       Report bugs to: last (ATmark) cbrc (dot) jp
       LAST home page: http://last.cbrc.jp/

lastal 199							     May 2012								 LASTAL(1)

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

diff 2 files; output diff's to 3rd file

Discussion started by: blt123

2. Shell Programming and Scripting

Find duplicates from multuple files with 2 diff types of files

Discussion started by: ricky007

3. AIX

diff command

Discussion started by: ariec

4. Shell Programming and Scripting

diff command help

Discussion started by: jaydeep_sadaria