forked from biogo/hts
-
Notifications
You must be signed in to change notification settings - Fork 0
/
flag.go
59 lines (52 loc) · 1.96 KB
/
flag.go
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
// Copyright ©2012 The bíogo Authors. All rights reserved.
// Use of this source code is governed by a BSD-style
// license that can be found in the LICENSE file.
package sam
// A Flags represents a BAM record's alignment FLAG field.
type Flags uint16
const (
Paired Flags = 1 << iota // The read is paired in sequencing, no matter whether it is mapped in a pair.
ProperPair // The read is mapped in a proper pair.
Unmapped // The read itself is unmapped; conflictive with ProperPair.
MateUnmapped // The mate is unmapped.
Reverse // The read is mapped to the reverse strand.
MateReverse // The mate is mapped to the reverse strand.
Read1 // This is read1.
Read2 // This is read2.
Secondary // Not primary alignment.
QCFail // QC failure.
Duplicate // Optical or PCR duplicate.
Supplementary // Supplementary alignment, indicates alignment is part of a chimeric alignment.
)
// String representation of BAM alignment flags:
// 0x001 - p - Paired
// 0x002 - P - ProperPair
// 0x004 - u - Unmapped
// 0x008 - U - MateUnmapped
// 0x010 - r - Reverse
// 0x020 - R - MateReverse
// 0x040 - 1 - Read1
// 0x080 - 2 - Read2
// 0x100 - s - Secondary
// 0x200 - f - QCFail
// 0x400 - d - Duplicate
// 0x800 - S - Supplementary
//
// Note that flag bits are represented high order to the right.
func (f Flags) String() string {
// If 0x01 is unset, no assumptions can be made about 0x02, 0x08, 0x20, 0x40 and 0x80
const pairedMask = ProperPair | MateUnmapped | MateReverse | MateReverse | Read1 | Read2
if f&1 == 0 {
f &^= pairedMask
}
const flags = "pPuUrR12sfdS"
b := make([]byte, len(flags))
for i, c := range flags {
if f&(1<<uint(i)) != 0 {
b[i] = byte(c)
} else {
b[i] = '-'
}
}
return string(b)
}