Skip to main content

Table 8 Pipeline assembly comparison; Canu assembler is used with MHAP and HISEA as read aligners

From: HISEA: HIerarchical SEed Aligner for PacBio data

Genome

Parameter

Canu + MHAP

Canu + HISEA

  

30x

50x

30x

50x

E.coli

Contig #

7

3

8

1

 

NG50

2,771,323

3,969,196

1,223,211

4,642,165

 

Max contig

2,771,323

3,969,196

1,525,215

4,642,165

 

% Ref

99.85

99.97

99.82

100.00

 

Avg idy

99.97

99.99

99.97

99.99

 

Breakpoints

3

3

3

3

S.cerevisiae

Contig #

43

31

35

29

 

NG50

540,299

687,498

682,168

774,485

 

Max contig

964,505

1,534,125

1,537,586

1,534,133

 

% Ref

98.90

99.35

99.12

99.58

 

Avg idy

99.81

99.88

99.82

99.88

 

Breakpoints

17

14

17

14

C.elegans

Contig #

393

170

127

133

 

NG50

636,401

1,987,017

2,140,282

2,032,954

 

Max contig

2,648,207

4,224,025

4,227,561

5,669,072

 

% Ref

96.00

99.84

99.81

99.80

 

Avg idy

99.76

99.91

99.85

99.91

 

Breakpoints

431

423

390

435

A.thaliana

Contig #

159

99

140

122

 

NG50

3,331,858

6,715,370

5,069,662

8,124,422

 

Max contig

12,892,206

14,177,369

12,890,806

15,940,320

 

% Ref

92.22

92.55

92.37

92.51

 

Avg idy

99.17

99.22

99.17

99.22

 

Breakpoints

2,550

2,693

2,680

2,704

D.melanogaster

Contig #

597

390

553

372

 

NG50

1,933,939

4,983,913

6,417,268

13,672,005

 

Max contig

8,238,062

17,900,724

17,366,974

25,767,672

 

% Ref

95.08

98.55

96.47

98.65

 

Avg idy

99.80

99.89

99.80

99.87

 

Breakpoints

1,039

1,383

1,254

1,461

  1. Two coverage levels, 30x and 50x, are used for each genome. The best results are shown in bold