Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Going from FASTA scaffolds and FASTQ reads to AMOS BANK

    Nowadays, most NGS assembly pipelines end in a scaffolds FASTA file. Typically, we have no assembly files (ace, contig, afg, among others). Thus, all we have to work is a scaffolds file (FASTA) and the original reads (typically FASTQ). SAM files could be easily generated by aligning the reads to reference. The question is: is it posible to use these information to load an AMOS BANK and run some very useful pipelines, such as amosvalidate? Has anyone scripted something (sam2ace, sam2afg, sam2amos, bam2ace, bam2amos or anything) and is willing to share it?

    I think it would be very interesting to develop such piece of software and it would be very useful, so we could take advantage of the powerful AMOS tools. I understand that there are technical problems in going from SAM to an assembly format. For example, what to do with reads that more than one mapping? Proper read pairing information is also crucial to get an assembly format that is useful to amosvalidate, for example. There might be many other problems in this conversion that I don't know yet.

    Any other ideas do you want to share? If this tool does not exists, I would be very willing to develop and share it.

  • #2
    Hi Glacerda. I've been trying to figure this out recently as well. You posted some time ago but this may still be relevant to yourself/others:

    #ABYSS
    KAligner --seq -mk *K* reads.fa contigs.fa | abyss2afg contigs.fa >out.afg

    #AMOS
    bank-transact -cb out.bnk -m out.afg
    hawkeye out.bnk



    supposedly toAmos.pl (AMOS) or addSolexaReads (Consed) also work - but I haven't successfully used them

    Comment

    Latest Articles

    Collapse

    • 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
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 11:49 AM
    0 responses
    15 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-24-2024, 08:47 AM
    0 responses
    16 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    61 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    60 views
    0 likes
    Last Post seqadmin  
    Working...
    X