Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Nomijill
    replied
    This is great thank you.

    Leave a comment:


  • joa_ds
    replied
    Hi there, i am under some pressure to deliver results from a solexa run. They did a run, and all i had to do was give them a list of variants.

    Problem is, i ran bowtie, give them the list. They verify with Sanger sequencing, and the results are different. So i used MAQ then. Results different than BowTie results (but also different than the 'reality' obtained with Sanger).

    So I have a lot of questionmarks... If results vary depending on the mapping software, well, my neighboor's genome has not changed since i pressed the MAQ button, so that is a problem.

    You guys still think using Eland will give me some 'strength' in explaining i did the best possible mapping i could do?

    Problem is that they have multiple genotypes mixed and are looking for QTL's, so a 40% or 60% variant, in the 'normal world' is just a heterozygous variant, in QTL world, it is a big difference...

    Leave a comment:


  • der_eiskern
    replied
    Originally posted by regyre View Post
    Hi lh3,

    Thanks for the original benchmarking! I'm actually looking for such a benchmarking including the latest tools, like ZOOM!, bowtie, R Biostings pairwiseAlignment(), etc. Has anybody heard of that?

    Cheers,

    N.
    I second this. I'd be interested in a comparison of erroneous mappings with the dominant WGS alignment-to-ref tools.

    Leave a comment:


  • dan
    replied
    Originally posted by ECO View Post
    I'm definitely out of my league in this discussion, but if anyone needs hosting for some of these sample datasets, let me know!
    Can we make a wiki to store all the datasets and results? A nice table would be soo cool! Also, I think people contributing data and analysis to a wiki page is more 'cite-able' than forum posts, although this is clearly an effective mechanism for community discussion. A wiki page to summarize this thread could be a good 'citation'. I tried to find a suitable spot on Wikipedia, but they don't even have an article for next-generation sequencing yet, never mind short read aligners to cope with the data.

    I'd like to start benchmarking the de-novo aligners (where accuracy is much harder to pin down). This will become more and more of an issue as all the genome sequencing projects move over to next-gen technology.


    BTW, anyone seen this:

    I’ll be attending the coveted Marco Island meeting early next month (February 4-8), where I’ll present a poster on my evaluations of short read aligners for next-gen sequencing data



    Looks relevant to the topic of this thread.
    Last edited by dan; 01-22-2009, 06:09 AM. Reason: added a note about looking for a suitable page on Wikipedia

    Leave a comment:


  • regyre
    replied
    Where to find a recent benchmarking?

    Hi lh3,

    Thanks for the original benchmarking! I'm actually looking for such a benchmarking including the latest tools, like ZOOM!, bowtie, R Biostings pairwiseAlignment(), etc. Has anybody heard of that?

    Cheers,

    N.

    Leave a comment:


  • valeu
    replied
    Hi Heng!

    Have you heard about SOCS (http://bioinformatics.oxfordjournals...tract/btn512v1, http://socs.biology.gatech.edu/)? In their article they say that "The overall algorithm is similar to that used by software tools developed for analysis of Illumina-Solexa data (Li et al, 2008; Smith et al, 2008)

    I'm interested in alignment of SOLiD data and I'd like to know your opinion what to use.. Maq, Mosaik, SHRiMP, ZOOM or this new tool SOCS ..

    Best regards,
    Valentina

    Leave a comment:


  • jhui
    replied
    SeqMap (http://biogibbs.stanford.edu/~jiangh/SeqMap/) - work like ELand, can do 3 or more bp mismatches and also insdel

    Leave a comment:


  • BaCh
    replied
    Originally posted by Amit View Post
    ...
    MIRA
    ...
    Although the version 2.9.95 doesnt support SNP analysis yet but its compact.
    ...
    It does since 2.6 or something. Alas, the docs need to catch up.

    Leave a comment:


  • zee
    replied
    I am keen to find out what the optimum set of parameters are for MAQ in situations where we expect to find more indels as welll as the 0-2 mismatch hits.
    The latest version sounds like it would be better for minimizing false positives.

    Leave a comment:


  • ShaunMahony
    replied
    Heng,
    If you wanted to test RMAPQ, you could always convert the FASTQ files into an approximation of the PRB files:
    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

    Shaun

    Leave a comment:


  • cgb
    replied
    Eland is going to be hard to beat. It has had a few years of hard work, optimisation and thought put into it by Anthony (SSAHA) Cox at Solexa. It was designed from day 1 for aligning far more reads than you currently get from a GAII, to a full human reference on a desktop computer.

    Leave a comment:


  • lh3
    replied
    Originally posted by myrna View Post
    -H FILE dump multiple/all 01-mismatch hits to FILE [null]
    If you specify this option, maq will dump all hits to a gzip file. -C specifies how many hits to output.

    Leave a comment:


  • myrna
    replied
    Originally posted by lh3 View Post
    The latest version, 0.6.6, can output ALL hits with 0- or 1-mismatch in the seed.
    Hi Heng.
    I was excited to see this feature added to MAQ in the latest version as much of my work is applied to RNA (hence, it is quantitative). This should (hopefully) allow me to reduce some biases introduced by losing reads which map to multiple locations. Now, I am wondering, how do I go about using this feature? Is there a new option when running maq map? Or mapview? I have been unable to find it in the manpage.

    Thanks,

    Ryan

    FOLLOWUP:

    I found out the answer to this in the latest doc provided with version 0.6.6 (not the version on the sourceforge page).

    Usage: maq map [options] <out.map> <chr.bfa> <reads_1.bfq> [reads_2.bfq]

    Options: -1 INT length of the first read (<64) [0]
    -2 INT length of the second read (<64) [0]
    -m FLOAT rate of difference between reads and references [0.001]
    -e INT maximum allowed sum of qualities of mismatches [70]
    -d FILE adapter sequence file [null]
    -a INT max distance between two paired reads [250]
    -n INT number of mismatches in the first 24bp [2]
    -M c|g methylation alignment mode [null]
    -u FILE dump unmapped and poorly aligned reads to FILE [null]
    -H FILE dump multiple/all 01-mismatch hits to FILE [null]
    -C INT max number of hits to output. >512 for all 01 hits. [250]
    -s INT seed for random number generator [random]
    -N record mismatch positions (max read length<=55)
    -t trim all reads (usually not recommended)
    -c match in the colorspace
    Last edited by myrna; 05-05-2008, 02:00 PM.

    Leave a comment:


  • Amit
    replied
    Mira

    Hello everybody,

    I am new to this group and couldnt resist myself to join this exiciting discussion

    Has anybody heard of MIRA ?



    There is this guy quitely working on another software tool for Next Gen assembly .

    The USP of this tool is it can perform a true hybrid assembly SAnger+454 or 454+Solexa which I believe will solve the Next gen assembly issues.

    Although the version 2.9.95 doesnt support SNP analysis yet but its compact.

    This tool might be on slower side becuase it performs assembly iterative correcting errors on the way.

    I hope somebody evaluates this new version becuase I dont have the much needed hardware to run this program.

    regds,
    Amit

    Leave a comment:


  • lh3
    replied
    The latest version, 0.6.6, can output ALL hits with 0- or 1-mismatch in the seed.

    Leave a comment:

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM
  • seqadmin
    Exploring Human Diversity Through Large-Scale Omics
    by seqadmin


    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
    06-25-2024, 06:43 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 07-19-2024, 07:20 AM
0 responses
119 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-16-2024, 05:49 AM
0 responses
105 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-15-2024, 06:53 AM
0 responses
103 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
43 views
0 likes
Last Post seqadmin  
Working...
X