Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • localized consensus build instead of scaffolding

    I am trying to make something out of a denovo, heterozygous plant WGS.
    it won't go past a contig stage (scaffolding improved average contig length by 1 nt!) for several reasons.
    However, the ref seq (species of one of the parents) is covered to 85% with contigs up a few hundred contigs deep in places. I would like to create localized consensuses by collapsing all those reads into one stretch of consensus, without "inserting" the refseq (like GATK alternate ref maker would do). this should be possible by checking bedtools genome coverage for stretches with coverage, isolating all reads for that region from SAM file, then creating a consensus for that stretch, essentially making a scaffold. This would make it easier to make some use of the assembly attempt, for instance for blasting against to find certain genes or promotor regions etc. however for an entire plant genome even I will run out of patience trying to do this by hand!!

    Anybody good with python or other languages that could write a script? like I said, it should be possible by finding uninterrupted coverage stretches from bedtools genome coverage, grouping reads with their alignment info from SAM files and then using something like GATK alternate reference maker on each stretch.
    Maybe all that's needed would be an executable Linux wrapper script?
    Any suggestions? maybe there is a tool already out there that I missed stumbling upon?

  • #2
    Was it a hiseq 2x100 bp lane run?

    So what sequencing technology had been used and total volume of data after adapter/quality trimming?

    Assuming it is some sort of the illumina run and no much improvement after scaffolding, I assume it is truseq or nextera shotgun library. Ideally it should have been 2x250 on hiseq or a few miseq runs (2x300), but in quite a few cases people just do 2x100 or 2x125 cycle hiseq, and then wonder what to do with millions of contigs with N50<1kb...

    If you are really keen of a few particular genes -

    I would start with making legacy blast database from the preprocessed raw reads, aiming for the coverage of 25-50x, (divide and manage it with *.nal files, if formatdb volume count is over 99), stick it on NVMe SSD :-) and try some blasts of 1-2 kb probe dna or aa sequences (use megablast -n T option and up the wordsize to 32 - 64 bp...)

    If longer probe sequence - then divide onto 1-2kb chunks and conquer.

    If the range is still to short - try nextera matepair, Pacbio, Molecullo, 10K, BAC libraries.

    If you wish than you can pull out the matching reads and try assembling them with any assembler (even PHRAP works ok on up to 100K illumina reads).

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      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...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X