Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • KaiYe
    replied
    Originally posted by dingxiaofan1 View Post
    How about bwa-maq (http://bio-bwa.sourceforge.net/bwa.shtml), Pindel(seems not available now)
    Pindel is still being actively updated.

    New version of Pindel is being tested with the following additional functions:
    1. Allow sequence errors/SNPs in the same reads containing INDELs/SVs
    2. non-template sequence in deletions
    3. inversions
    4. tandem duplications
    5. breakpoints of large insertions

    Leave a comment:


  • drio
    replied
    Originally posted by joa_ds View Post
    Well, I can tell you the following.

    I have simulated loads of data. I have generated datasets with more than a billion reads where error profiles/PCR errors/sequencing errors were simulated and acted as a background for real SNPs that needed to be discovered. The real SNPs always come out on top. So far so good, but with deletions/insertions, I am completely confused.

    Indeed, the manual/documentation is not clear. Did you have any success contacting the authors?
    In your simulations do you introduce indels or it is only single base changes?
    How does that affects your results?

    Leave a comment:


  • dingxiaofan1
    replied
    Bwa-maq

    How about bwa-maq (http://bio-bwa.sourceforge.net/bwa.shtml), Pindel(seems not available now)

    Leave a comment:


  • joa_ds
    replied
    I just checked out their site and found the following thing.

    SOAPindel: SOAPindel is developed to find the insertion and deletion specially for re-sequence technology : Coming soon.

    I am just wondering if indel detection is already included in SOAP2? Has anyone ever seen an indel being detected by using SOAP2?

    Leave a comment:


  • francesco.vezzi
    replied
    Originally posted by joa_ds View Post
    Well, I can tell you the following.
    Indeed, the manual/documentation is not clear. Did you have any success contacting the authors?
    Yes I contact them more than once.... but they never reply...
    It is really strange that the most used aligner (or at least one of the most used) is so obscure.....

    Leave a comment:


  • joa_ds
    replied
    Well, I can tell you the following.

    I have simulated loads of data. I have generated datasets with more than a billion reads where error profiles/PCR errors/sequencing errors were simulated and acted as a background for real SNPs that needed to be discovered. The real SNPs always come out on top. So far so good, but with deletions/insertions, I am completely confused.

    Indeed, the manual/documentation is not clear. Did you have any success contacting the authors?

    Leave a comment:


  • francesco.vezzi
    replied
    Hi
    I can't help you so much but I can report another strange behaviour of SOAP2. Despite the fact that you vary the -v option that is the total amount of errors (mismatches) allowed in one read the output that you obtain is always the same. I Try to change several parameters like the seed length without understand how to specify an exact number of errors.

    Despite the fact that SOAP2 is one of the most used aligners it has some troubles: there is no a clear documentation, the source code is not available and hence we actually don't know what kind of heuristics are applied in order to speed up the query time...

    Francesco

    Leave a comment:


  • joa_ds
    started a topic SOAP2 and deletions/insertions

    SOAP2 and deletions/insertions

    I was told that SOAP should be able to detect deletions and insertions up to 3bp.

    So i made a reference file and a file with a normal read and some reads with 1,2,3 bp deleted and inserted. Just to test how the output format etc would look like and to get optimal parameter settings.

    i have used different -e settings etc, but it appears that not a single read with deletion/insertions gets aligned...

    All i get is

    Begin Alignment ...
    8 ok 0.00 sec
    Total Reads: 8
    Alignment: 1 (12.50%)

    -> Only the normal read gets aligned...

    Is there a secret parameter i should set somewhere?

Latest Articles

Collapse

  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM
  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 02-28-2024, 06:12 AM
0 responses
22 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-23-2024, 04:11 PM
0 responses
70 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
78 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
67 views
0 likes
Last Post seqadmin  
Working...
X