I have both short (illumina) and (pacbio) long reads of the same sample and I want to align it to a reference genome. I know there are many aligners, but I was wondering if there is one that uses a hybrid method, i.e. uses both short and long reads to improve the accuracy. I know that unicycler, for example, does a hybrid de novo assembly. However, I am looking for a program that does hybrid alignment to a reference genome. My google search did not find anything. Any comment would be appreciated.
1 answer
No such hybrid aligner exists since aligners perform read alignment independently for each fragment. That is: per read for long read, and per read pair for paired-end short read sequencing. A 'hybrid' approach would give the same results as aligning the short and long reads independently with the same aligner.
Short and long reads can be combined by 'hybrid' methods either prior to alignment (using short reads to correct long reads, or perform de novo assembly), or after alignment (reference-guided assembly/alignment adjustment/variant calling based on the aligned short+long reads).
Log in to answer this question.
-
Fish_genome_assembly
written by manaswiniparija3 •
-
Rotate circular contigs to replication initiator protein (ex. dnaA, repA)
written by SergeW •
-
Quallity trimming by phred score, how far is to far?
written by robert.murphy •
-
PacBio long-reads impact in transcriptome de novo assembly
written by Manuel Mendoza •
-
Variant calling from short/long read hybrid BAM
written by Vitis
-
Illumina polishing after PacBio assembly and polishing?
written by ikangkim •
-
MUMmer alignments: Trying to understand the average identity result
written by caro-ca •
-
Unicycler generating short contigs during hybrid assembly
written by hermoinejean •
-
Improve illumina short read assembly using PacBio long reads
written by shachigahoimbi •
-
Both long and short read of the same sample
written by angela.tithi23 •
SPAdes or MIRA? you can pass short and long reads plus a reference sequence
@JC I guess you are referring to reference based assembly but OP is asking about mapping.
MIRApage says:Technically
minimap2could be used for both reads but not at the same time.Not sure how that will improve mapping accuracy. Reads are aligned one at a time.
I am confused about what are the requirements then. OP, please clarify
Sorry for not being clear. @genomax is right. I am looking for a program that accepts and uses both short and long reads to improve mapping accuracy.
Thanks. I will take a look at minimap2.