#
Align Sequences to a Genome
After your sequences (in FASTQ format) have been checked for quality, you will need to align them to a reference genome before you can call variants. Harpy offers several aligners for this purpose:
Only aligned sequences
The aligners in this module only output aligned sequences, meaning unmapped reads will not appear in the results. Keep that in mind when you read the alignment summary report and see "100%" aligned reads. That just means there are no unmapped reads in the BAM file, which is the result of aligner output configuration and not biological processes.
Despite the fact that EMA is the only barcode-aware aligner offered, when using BWA or strobealign, Harpy retains the barcode information from the sequence headers and will
assign molecule identifiers (MI:i
SAM tags) based on these barcodes and the molecule distance threshold.