Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SNP calling problems

    Hi

    I am carrying out alignments using BWA and Bowtie and then using samtools to produce some fairly large SAM files (>100G) followed by calling SNPs and producing VCF files.

    I have found that BWA seems to be taking many days longer to run than Bowtie. However, for test sequences BWA seemed to be sufficient for our purposes, the read depth of the aligned sequences to the reference sequences was sufficient to call SNPs which pass our statistical tests. Bowtie, however, was much faster but results often only had a read depth of 5 or so and much fewer SNPs were called and those that were did not pass stats tests.

    I thought reducing the stringeny of Bowtie somehow might help; I used -v 3 and -y with resulted in a read depth which seemed to be slightly higher. A few more SNPs were called but this was not sufficient.

    Are there any other tools I can use which will be faster than BWA or other Bowtie parameters I could explore? (I am expecting a read depth more towards 50-60) as was originally possible for the test sequences when using BWA.

    For certain datasets am I also receiving errors in the sai -> sam conversion stage. With BWA taking so long to run and then producing these errors it is taking much longer to carry out these tests and then re-run them than originally expected giving our success with the test sequences. There are certain datasets where I just cannot get past this stage. Now the input files contain multiple sequences instead of just one, could this be an issue? Would perhaps piping the output of BWA to samtoools for the production of the SAM file help to stop this error?

    I am really at loss as to why scaling up the process from the test sequences has resulted in these problems. I can give more biological context if required. Hope you can help!

    Thanks

    Jom

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-10-2024, 07:30 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 09:45 AM
0 responses
201 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 08:54 AM
0 responses
211 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-02-2024, 03:00 PM
0 responses
193 views
0 likes
Last Post seqadmin  
Working...
X