File: maqview.pod

package info (click to toggle)
maqview 0.2.5-7
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 1,076 kB
  • ctags: 1,771
  • sloc: ansic: 13,310; cpp: 1,295; java: 178; perl: 116; sh: 41; makefile: 22
file content (159 lines) | stat: -rw-r--r-- 4,544 bytes parent folder | download | duplicates (4)
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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
=head1 NAME

Maqview - Maq alignment indexer and viewer

=head1 SYNOPSIS

B<maqindex> B<-i>|B<-v>|B<-b> [B<-c> F<in.cns>] F<in.map>
[I<chr>[:I<start>[-I<end>]] [...]]

B<maqview> [B<-c> F<in.cns>] I<in.map>

=head1 DESCRIPTION

B<Maqview> is a set of programs that achieves fast random access to the
alignment file generated by `B<maq> map', and displays the alignment
with a nice GUI (Graphical User Interface). It is highly efficient in
speed, memory and disk usage. B<Maqview> is based on OpenGL and is known
to work on both Mac OS X and Linux. Porting to Windows is in principle
easy.

=head2 B<maqindex>

B<maqindex> B<-i>|B<-v>|B<-b> [B<-c> F<in.cns>] F<in.map>
[I<chr>[:I<start>[-I<end>]] [...]]

=over 9

=item B<-i>

Index the alignment file

=item B<-v>

Print the alignment in the specified region in the `maq mapview' format

=item B<-b>

Dump the alignment in the specified region in the binary .map format

=item B<-c> F<FILE>

Index the Maq consensus file for viewing in maqview at the same time

=back

Program B<maqindex> indexes a maq alignment file F<in.map> or quickly
retrieves all the reads in one or multiple regions. One of B<-i>, B<-v>
and B<-b> must be used. For B<-v> or B<-b>, at least one region like
`chrX', `chrX:1000' or `chrX:1,000-2,000' must be specified. Multiple
regions are allowed.

=head2 B<maqview>

B<maqview> [B<-c> F<in.cns>] I<in.map>

=over 9

=item B<-c> F<FILE>

The Maq consensus file

=back

Program B<maqview> displays the read alignment in a graphical
window. When F<in.cns> is specified, the top sequence is the reference,
followed by the Maq consensus sequence; otherwise, the majority-rule
consensus will be calculated.

Maqview has two views: sequence view and box view. In the sequence view,
read sequences will be printed on the screen. Darker bases indicate
lower base qualities and red ones show the differences in comparison to
the reference if F<in.cns> is specified or to the majority-rule
consensus if not. In the box view, different types of nucleotides are
represented as colour boxes with green for A, cyan for C, orange for G,
red for T and dark gray for N. The saturation of colours indicates the
base qualities and the thickness lines of reads shows the mapping
qualities of read alignments. Zooming in/out is supported only in the
box view. In both views, the status bar at the bottom of the window will
show some information about the key touches, and read names and base
qualities pointed by the mouse.

Navigating the alignment is accomplished by vaious key bindings.

=over 2

=item

Key bindings:

 ?                Display help
 <Esc>            Exit
 <F1>             Sequence view
 <F2>             Box view
 p                Switch to the previous reference sequence
 n                Switch to the next reference sequence
 r                Refresh

 <digits> <Enter> Move to the position pointed by <digits>
 h / <Left>       Move left by one base
 l / <Right>      Move right by one base
 k / <Up>         Move up by one line
 j / <Down>       Move down by one line
 >                Move to the next read
 <                Move to the previous read
 g / <Home>       Move to the beginning of the current reference
 G / <End>        Move to the end of the current reference
 u / <PageUp>     Move left by one page
 <Space>/<PageDn> Move right by one page
 <Shift> <Arrow>  Move by 100 bases
 <Ctrl> <Arrow>   Move by 1000 bases

 <Ctrl> +         Add a new view
 <Ctrl> -         Remove the current view

 <Enter>          Toggle on/off right auto-scrolling
 ^b               Toggle on/off left auto-scrolling
 +                Zoom in
 -                Zoom out
 0                Zoom to the default scale
 q                Toggle show/hide mapping qualities
 e                Toggle display SE/PE mapping qualities
 o                Open a new alignment file

=item

Enter an integer and then press <Enter> to jump to the required
position.

=item

Enter an integer and then press a key other than <Enter> to speed up
scrolling by the specified factor.

=item

Left click and drag the yellow cross in the square to move along the
sequence. Right click the reference to highlight a column.

=back

=head1 LICENSE

GNU General Public License (GPL)

=head1 AVAILABILITY

L<http://maq.sourceforge.net>

=head1 AUTHORS

Jue Ruan <ruanjue@genomics.org.cn> for writing and maintaining the
whole software.

Heng Li <lh3@sanger.ac.uk> for testing the program and drafting the
documentations.

Mengyao Zhao <zhaomengyao@gmail.com> for designing the color schemes.

=cut