forked from minillinim/ionPairer
/
README
75 lines (56 loc) · 4.07 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
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************
d8b 8888888b. d8b
Y8P 888 Y88b Y8P
888 888
888 .d88b. 88888b. 888 d88P 8888b. 888 888d888 .d88b. 888d888
888 d88""88b 888 "88b 8888888P" "88b 888 888P" d8P Y8b 888P"
888 888 888 888 888 888 .d888888 888 888 88888888 888
888 Y88..88P 888 888 888 888 888 888 888 Y8b. 888
888 "Y88P" 888 888 888 "Y888888 888 888 "Y8888 888
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************
Scrimpy way to scaffold an assembly using ion torrent reads.
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************
Usage:
Take your existing assembly, bin out the contigs of interest and then
run bwa-long twice (once for each of the forward and reverse reads) to map your
ion MP reads against your contigs.
Then simply run ionPairer.pl script. There are a few hard-coded options lurking
within but the up and up is that it will detect pairs, wrongly mapped regions
unpaired reads etc. It can determine orientation and insert size automagically.
The most useful output is a .gv file which you can use to visualise the
pairings using graphviz. Use this visualisation to remove spurious pairings and
redraw etc. This script is quite useful for making rough scaffold maps.
Once edited the .gv file, use the scaffolder.pl script to create the scaffolds.
The script will orient the contigs (forward or reverse complement) and scaffold
the contigs with 25 Ns between each contig.
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************
Project home page
Info on the source tree, documentation, issues and how to contribute, see
http://github.com/minillinim/ionPairer
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************
Cite
This software is currently unpublished.
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************
Copyright © 2012 Michael Imelfort. See LICENSE.txt for further details.
*******************************************************************************
*******************************************************************************
*******************************************************************************
*******************************************************************************