Extract the sizes of all sequences in a paf alignment

chrom_sizes(ali)

Arguments

ali

pafr or tibble containing the genome alignment (as returned by read_paf)

Value

list with two elements (tlens and qlens) Each element is a dataframe with one column of sequence names and another column containing the length of each sequence

Examples

ali <- read_paf( system.file("extdata", "fungi.paf", package="pafr") ) chrom_sizes(ali)
#> $qlens #> qname qlen #> 1 Q_chr1 6273420 #> 388 Q_chr2 5852505 #> 792 Q_chr3 4992232 #> 1226 Q_chr4 4808386 #> 1654 Q_chr5 4708905 #> 2005 Q_chr6 4343484 #> 2315 Q_chr7 3141109 #> 2486 Q_chrm 85999 #> #> $tlens #> tname tlen #> 1 T_chr2 5893594 #> 2 T_chr5 4620715 #> 59 T_chr6 3360039 #> 64 T_chr4 4649703 #> 75 T_chr1 6377790 #> 76 T_chr3 5539885 #> 79 T_chr7 3267564 #> 2486 T_chrm 111040 #>