SAUTE: sequence assembly using target enrichment.
BACKGROUND
Illumina is the dominant sequencing technology at this time. Short length, short insert size, some systematic biases, and low-level carryover contamination in Illumina reads continue to make assembly of repeated regions a challenging problem. Some applications also require finding multiple well supported variants for assembled regions.


RESULTS
To facilitate assembly of repeat regions and to report multiple well supported variants when a user can provide target sequences to assist the assembly, we propose SAUTE and SAUTE_PROT assemblers. Both assemblers use de Bruijn graph on reads. Targets can be transcripts or proteins for RNA-seq reads and transcripts, proteins, or genomic regions for genomic reads. Target sequences are nucleotide and protein sequences for SAUTE and SAUTE_PROT, respectively.


CONCLUSIONS
For RNA-seq, comparisons with TRINITY, RNASPADES, SPALIGNER, and SPADES assembly of reads aligned to target proteins by DIAMOND show that SAUTE_PROT finds more coding sequences that translate to benchmark proteins. Using AMRFINDERPLUS calls, we find SAUTE has higher sensitivity and precision than SPADES, PLASMIDSPADES, SPALIGNER, and SPADES assembly of reads aligned to target regions by HISAT2. It also has better sensitivity than SKESA but worse precision.
['Souvorov|Alexandre|A|', 'Agarwala|Richa|R|http://orcid.org/0000-0002-5518-9723']
D000465:Algorithms D016678:Genome D023281:Genomics D059014:High-Throughput Nucleotide Sequencing D000081246:RNA-Seq D017422:Sequence Analysis, DNA
Illumina reads de Bruijn graphs RNA-seq Antimicrobial resistance De-novo assembly