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
      Recent Advances in Sequencing Analysis Tools
      by seqadmin


      The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
      05-06-2024, 07:48 AM
    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin




      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
      04-22-2024, 07:01 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 05-14-2024, 07:03 AM
    0 responses
    20 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-10-2024, 06:35 AM
    0 responses
    44 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-09-2024, 02:46 PM
    0 responses
    54 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-07-2024, 06:57 AM
    0 responses
    43 views
    0 likes
    Last Post seqadmin  
    Working...
    X