After assembly, 63,485 gene models were obtained (average gene size 882 bp; N50 = 1485 bp), among which 38,540 unigenes (60.70 % of total gene models) were annotated by comparing them with four public databases (Nr, Swiss-Prot, KEGG and COG).
After assembly, 63,485 gene models were obtained (average gene size 882 bp; N50 = 1485 bp), among which 38,540 unigenes (60.70 % of total gene models) were annotated by comparisons with four public databases (Nr, Swiss-Prot, KEGG and COG): 38,338 unigenes (60.38 % of total gene models) showed high homology with sequences in the Nr database.
In total, 78,806 contigs greater than 1 kb in length are recovered, and this size, while at its lower limit too small for complete gene models (centipede S. maritima median total gene size: 3,170 bp), is sufficient to span most protein-coding domains.
594 (61%) of gene models showed a more than 10% increase in gene size while only 28 33%) of validated gene models showed a more than −10% change in gene size.
To determine the probability of observing multiple de novo truncating mutations in ASXL3, we developed a model [ 21, 22] accounting for gene size, GC content and de novo rates of SNVs and small insertion/deletions, and of the probability of those mutations causing a truncation of the protein.
We briefly discuss the effects of the other significant gene size and methylation terms in this model.
These criteria result in an observed absolute change in expressed gene size for ∼62% of the Myc-BestModel gene model set.
When the pan-genome size and core genome size of n strains (1 ≤ n ≤ N) were available, a series of mathematics models (Supplementary Material) were used to depict pan-genome size, core genome size and new gene size.
Of the 974 modified gene models, 453 (30%) had no error flags within their FDM sequences, but showed a more that 10% increase in gene size.
First we consider a simple model in which genes are assumed to have uniform size, and then extend it to variable gene sizes.
