forked from tseemann/PEAR
-
Notifications
You must be signed in to change notification settings - Fork 0
/
README
82 lines (62 loc) · 2.6 KB
/
README
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
____ _____ _ ____
| _ \| ____| / \ | _ \
| |_) | _| / _ \ | |_) |
| __/| |___ / ___ \| _ <
|_| |_____/_/ \_\_| \_\
....Paired-End reAd mergeR....
Authors: Jiajie Zhang, Kassian Kobert, Tomas Flouri, Alexandros Stamatakis
License: Creative Commons license
with
Attribution-NonCommercial-ShareAlike 3.0 Unported
Introduction:
-------------
PEAR assembles Illumina paired-end reads if the DNA fragment sizes are smaller
than twice the length of reads. PEAR can assemble 95% of reads with 35-bp mean
overlap with a false-positive rate of 0.004. PEAR also works with multiplexed
data sets where the true underlying DNA fragment size varies. PEAR has an
extremely low false-positive rate of 0.0003 on data sets where no overlap exists
between the two reads (i.e. when DNA fragment sizes are larger than twice the
read length).
For more information, requests and bug-reports visit our website at
http://www.exelixis-lab.org/web/software/pear
How to compile:
---------------
1. git clone https://github.com/xflouris/PEAR.git
2. cd PEAR
3. make
4. make install
How to run self-tests:
----------------------
1. Make sure you have python 2.4 (or newer) installed
2. Go to the "test" folder
3. type: ./test.py
This will run PEAR on several simulated data sets with various options to check
if the newly compiled program works properly.
How to use:
-----------
PEAR can robustly assemble most of the data sets with default parameters. The
basic command to run PEAR is
./pear -f forward_read.fastq -r reverse_read.fastq -o output_prefix
The forward_read file usually has "R1" in the name, and the reverse_read
file usually has "R2" in the name.
Output files:
-------------
PEAR produces 4 output files:
1. output_prefix.assembled.fastq - the assembled pairs
2. output_prefix.unassembled.forward.fastq - unassembled forward reads
3. output_prefix.unassembled.reverse.fastq - unassembled reverse reads
4. output_prefix.dicarded.fastq - reads which did not meet criteria specified in options
Advanced usage:
---------------
For further options and fine-tuning type
./pear -h
Important information:
----------------------
1. The input files must be in FASTQ format
2. PEAR does not check the paired-end reads names. PEAR assumes that the reads
in both files are in the same flowcell position if they appear on the same line
number. Therefore, the validity of the input files is left as a user
responsibility.
How to cite:
------------
J. Zhang, K. Kobert, T. Flouri, A. Stamatakis. PEAR: A fast and accurate Illimuna Paired-End reAd mergeR