Dear Nash,
You may try out the Subread aligner (http://subread.sourceforge.net), which can align reads as long as 1200bp.
Cheers,
Wei
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
blasr output
Mark,
That's great to know....I have printed out the help pages, but there are still unanswered questions for me...would you like to take this discussion offline or do you want me to post the questions right here? If there's a special PacBio support site for blasr, I can reach you through that...Please let me know your convenience. Thanks much,
Nash
Leave a comment:
-
Originally posted by naragam View PostPerhaps, I should really start a new thread...but, does anybody on this forum have good experience with blasr alignments to discuss the various options for the run and further the several output formats that are available. I have just started playing with some of the balsr runs and I have some pointed questions that I'd like to ask and/or seek detailed docs to refer to in terms of understanding all the options and outputs.
Any help available in this forum?
Thanks much in advance for any pointers,
Nash
Most of the help may be found by running blasr -h, or blasr -help for detailed help. There are many output formats including tabular ones for which you can get column labels with the -header option, human readable output (-m 0), and sam (specified by -sam).
-mark
Leave a comment:
-
PacBio "blasr" questions....
Perhaps, I should really start a new thread...but, does anybody on this forum have good experience with blasr alignments to discuss the various options for the run and further the several output formats that are available. I have just started playing with some of the balsr runs and I have some pointed questions that I'd like to ask and/or seek detailed docs to refer to in terms of understanding all the options and outputs.
Any help available in this forum?
Thanks much in advance for any pointers,
Nash
Leave a comment:
-
A quicker refined flavour of BLAT is BFAST :http://www.plosone.org/article/info:...l.pone.0007767
Leave a comment:
-
blasr compilation
Mark,
Am trying to compile blasr on my machine and am missing some header files in the tar file distribution. Can you please point me to sources who can help me or provide the *.h files needed? Thanks much,
Nash
Leave a comment:
-
Originally posted by naragam View PostYeah, I hope to run blasr soon but, in the meantime, I am trying to learn some of these long read tools that I haven't worked with before. Do you know if you have to have the fastq files for bwa sw?
Nash
You will want the bas.h5 files since they have additional information about subread coordinates.
Leave a comment:
-
Thank you Mark...I don't have access to hd5 files yet....they are with the core sequencing facility and I am not sure they will give me those right now.... But I am working with them to gradually get some of the pipeline tools locally on my new Ubuntu machine that still needs memory upgrades before I can run your pipeline tools...
Yeah, I hope to run blasr soon but, in the meantime, I am trying to learn some of these long read tools that I haven't worked with before. Do you know if you have to have the fastq files for bwa sw?
Nash
Leave a comment:
-
Hi Nash,
You can install blasr on your own using github (https://github.com/PacificBiosciences/blasr).
If you have the hdf files, there are options (-useccsdenovo) to align the ccs sequences instead of the raw subreads.
HTH,
-mark
Originally posted by naragam View PostThank you all for suggesting blat and bwa sw for aligning long reads. I am looking at the docs for bwa sw and it looks like in the command:
bwa bwasw database.fasta long_read.fastq >aln.sam
the parameter "database.fasta" is the reference and the parameter "long_read.fastq" is the sequence being aligned. Right?
So does it absolutely need the fastq file or can it just work w/o the quality data, i.e., just a *.fasta file? Also how about the ccs based output from PacBio? Anybody has tried the PacBio ccs outputs? I'm trying to get "blasr" tool from PacBio pipeline installed here, but I am not there yet...
Thanks in advance,
Nash
Leave a comment:
-
Re: Alignment tool for long reads
Thank you all for suggesting blat and bwa sw for aligning long reads. I am looking at the docs for bwa sw and it looks like in the command:
bwa bwasw database.fasta long_read.fastq >aln.sam
the parameter "database.fasta" is the reference and the parameter "long_read.fastq" is the sequence being aligned. Right?
So does it absolutely need the fastq file or can it just work w/o the quality data, i.e., just a *.fasta file? Also how about the ccs based output from PacBio? Anybody has tried the PacBio ccs outputs? I'm trying to get "blasr" tool from PacBio pipeline installed here, but I am not there yet...
Thanks in advance,
Nash
Leave a comment:
-
Originally posted by pacbio View PostHi All
We have actually developed a fast and accurate aligner named BLASR (Basic Local Alignment with Successive Refinement) - http://www.smrtcommunity.com/SMRT-An...gorithms/BLASR to align our long reads. The source code for this as well as the full analysis software suite is freely available at the same PacBio DevNet site. A publication on this algorithm is also currently in review, so stay tuned.
Leave a comment:
-
Hi All
We have actually developed a fast and accurate aligner named BLASR (Basic Local Alignment with Successive Refinement) - http://www.smrtcommunity.com/SMRT-An...gorithms/BLASR to align our long reads. The source code for this as well as the full analysis software suite is freely available at the same PacBio DevNet site. A publication on this algorithm is also currently in review, so stay tuned.
Leave a comment:
-
I used Blat to do reference based scaffolding by aligning contigs to scaffolds. So that should work with long reads I assume.
Leave a comment:
Latest Articles
Collapse
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Leave a comment: