Skip to main content

Table 4 Assembly statistics for the CamDro2; CamDro3 (Pilon) using one round of Pilon [51] for polishing; and CamDro3 (BBMap) using one round of variant calling with BBMap (https://sourceforge.net/projects/bbmap/) for polishing. Note that CamDro3 (BBMap) was chosen over CamDro3 (Pilon) as the final version of CamDro3 because of better BUSCO and RNA-Seq mapping percentages

From: Nucleotide diversity of functionally different groups of immune response genes in Old World camels based on newly annotated and reference-guided assemblies

 

Assembly

CamDro2

CamDro3

(Pilon)

CamDro3

(BBMap)

Total size

2,154,386,959

2,194,229,671

2,169,346,739

Gap length

20,603,579

17,930,821

17,043,352

Scaffolds

 Number

23,439

21,070

21,070

 Longest

124,992,380

125,472,505

124,715,342

 N90a

4,922,612

25,062,887

24,767,672

 L90b

31

32

32

 N50a

75,021,453

70,557,636

70,369,702

 L50b

11

12

11

Contigsc

 Number

45,969

41,934

53,085

 Longest

9,490,880

14,412,615

2,012,572

 N90

177,587

202,272

49,444

 L90

1944

1436

10,023

 N50

1,333,162

1,961,815

236,380

 L50

423

303

2637

 Single-copy BUSCOsd

3851

3853

3852

 Duplicated BUSCOs

24

23

25

 Fragmented BUSCOs

133

132

134

 Missing BUSCOs

96

96

93

 RNA-Seq Mapping Percentagee

88.30

90.36

92.04

  1. aN90/N50 are the scaffold or contig lengths such that the sum of the lengths of all scaffolds or contigs of this size or larger is equal to 90/50% of the total assembly length
  2. bL90/L50 are the smallest number of scaffolds or contigs that make up at least 90/50% of the total assembly length
  3. cUsing minimum gap length of 25 bp
  4. dBUSCOs: Benchmarking Universal Single-Copy Orthologs [19] are mammalian BUSCOs from OrthoDB v. 9.1 genes [20]
  5. eOverall mapping rates using HiSat v. 2.1.0 [53] of dromedary RNA-Seq reads from Sequence Read Archive accession: SRP017619 and Alim et al. [54]