Best Alignment tool of Denovo assembly Results.
1
0
Entering edit mode
7.0 years ago
clear.choi ▴ 30

I have Denovo Assembly sequence from PacBio sequence.

It is around 2.9 GB size of sequences.

I'd like to know what is best alignment tool to analyze this data.

I am usually using this "blasr" to align HG19 to Sub reads. However, When I align HG19 to Denovo assembly results, It throwing the error. It looks like too large sequence to align it.

Could anyone know what is best practice to do it? I want to see and have bam file to visualize in IGV program.

Thank you in advance.

Assembly sequence alignment SNP • 3.0k views
ADD COMMENT
1
Entering edit mode

If you are trying to align 2.9 GB (genome sized) sequence to another reference genome (hg19) then you should look at a tool like LASTZ which is meant for aligning chromosome sized fragments.

ADD REPLY
0
Entering edit mode

Bowtie and BWA are popular and widely used tools. Samtools will also become your friend. Example workflow.

ADD REPLY
0
Entering edit mode

I don't think this really addresses the question, which is about whole-genome alignment. Neither bowtie nor bwa are suitable for that purpose. Mummer IS suitable but I don't think it produces bam files. If Blasr causes problems, you might try daligner.

Alternatively, you could just align the error-corrected reads to the human genome. That would be much simpler and likely more productive.

ADD REPLY
0
Entering edit mode
7.0 years ago

You can probably achieve what you are trying to do easiest by aligning the pacbio contigs to the genome, not the genome to the contigs. Brian's comment is also very applicable, as is genomax's.

If you have really big contigs Bwa mem may still work, but use samtools and not Picard for converting and sorting contig size alignments from SAM to bam.

ADD COMMENT

Login before adding your answer.

Traffic: 848 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6