How To Estimate K And Lambda For Local Alignment'S Significance?
Hi,
I'm wondering how to estimate K and lambda for local alignment's significance for a specific scoring matrix and gap penalty. I know there's a paper (The statistical parameters for local alignment score distributions) about it, but I hope that there's already a tool I can use directly. Does anyone know such a tool? Thanks.
• 4,874 views
•
link
updated
by
Casey Bergman
•
written
by
Tianyang Li •
•
written
by
Tianyang Li •
1 answer
I had this same problem a number of years back. I emailed Stephen Altschul who sent me some C code to calculate K and lambda, which worked fine as I recall. Perhaps you can try this same thing and also see if Altschul would be willing to post this code more publicly.
• 1 views
•
link
written
by
Casey Bergman
Log in to answer this question.
More posts like this
-
Aligning very short sequences using bowtie2, I need fine control over how many mismatches are allow…
written by Mauro •
-
What are recommended parameters for the local protein alignment ?
written by Alexander •
-
Herald: The Biostar Herald for Monday, January 16, 2023
written by Biostar
-
How to find if sequences perfectly overlap
written by gero.knittel •
-
Forum: BWA command guide
written by Chen Sun
-
EBSeq explained to a biologist
written by Rik Verdonck •
-
Tutorial: InsideDNA: Analysing raw sequencing reads with FASTQC for quality control and filtering
written by shilparaopradeep •
-
comparing count data with a vector of quantitative Scores
written by Assa Yeroslaviz
-
Best Existing Tool For Computing The Overlaps Between Each Read In A Set Of Short Reads?
written by Tianyang Li •
-
Local Blast Returns No Result, While Using The Website There Are Some Matches
written by Zhizhong
I asked Dr. Altschul and he suggested that I use by Dr. Spouge's group. The tool they provided (alp1.3) couldn't be compiled with g++ 4.4, so I fixed a few things (added a few headers) and it can now be compiled using g++. The code is here.
I asked Dr. Altschul and he suggested that I use by Dr. Spouge's group for gapped alignment (ALP on http://www.ncbi.nlm.nih.gov/CBBresearch/Spouge/html_ncbi/html/index/software.html). The tool they provided (alp1.3) couldn't be compiled with g++ 4.4, so I fixed a few things (added a few headers) and it can now be compiled using g++. The that works with g++ is here.
I asked Dr. Altschul and he suggested that I use the tool by Dr. Spouge's group for gapped alignment (ALP on http://www.ncbi.nlm.nih.gov/CBBresearch/Spouge/html_ncbi/html/index/software.html). The tool they provided (alp1.3) couldn't be compiled with g++ 4.4, so I fixed a few things (added a few headers) and it can now be compiled using g++. The that works with g++ is here.