This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Spades Assembler output

hello,
I used Spades for assembling short reads.
in the Spades_output, i have two files "final_contigs.fasta" and "final_contigs.fastg".
what's the fastg format?
and what file should I take for my research?

thanks

assembly spades contig
• 5,821 views • link
updated by MathGon • MathGon • written by midox midox

you should really have a contigs.fasta(g) and scaffolds.fasta(g) in your output. The final_ are usually in the K(kmer) folders.

• 0 views • link
written by Adrian Pelin Adrian Pelin

2 answers

If you read carefully the SPAdes manual, you will find a link to the FASTG format specification. Its first paragraph:

This document introduces a prototype format for faithfully representing genome assemblies in the face of allelic polymorphism and assembly uncertainty. It is called FASTG, like FASTA, but the G stands for 'graph'.
• 0 views • link
written by h.mon h.mon

Yes that's true.
I have not seen.
thank you

• 0 views • link
written by midox midox

As mentioned, .fastg is a graph file. You can open it with Bandage to show lins between your contigs.

• 0 views • link
written by MathGon • MathGon

Log in to answer this question.

More posts like this
  • Error parsing attribute transcript_id (" " required)
    written by Yuyayuya • Yuyayuya
    <p>Hello everyone,</p> <p>I was using StringTie to get gene count and I got the error message: <br/> <a href="https://imgur.com/wdyZrPL"><img src="https://i.imgur.com/wdyZrPL.png" title="source: imgur.com" /></a></p> <p>After checking, …
  • Denovo assembly of Putative Plasmid
    written by jigarnt • jigarnt
    Hi All, I have illumina sequenced a putative plasmid with paired end library. Through sequencing I fetched two files, namely Plasmid_R1.fastq and Plasmid_R2.fastq. I am …
  • BLAT error sequnces
    written by midox midox
    <p>Hello,<br /> I use a BLAT alignment with two files.<br /> But here is the error &quot;<strong>Error: sequence name channel_231_read_27_template is repeated in the database, …
  • paired end illumina reads
    written by midox midox
    <p>Hello,<br /> I want to know if paired end sequencing we have two files and R2.fq and R1.fq.<br /> do R1.fq is the first strand …
  • short reads Assembly in contigs
    written by midox midox
    Hello, I made assembling illumina short reads of size 2 * 300bp and I used Velvet. After the assembly, I got contigs with size less …
  • unpacking zipped file in R
    written by zizigolu zizigolu
    <p>hello,</p> <p>i have 96 arrays in a zipped file, i did so to unpacking but recieved error, what happened please?</p> <p>i changed directory and typed:</p> …
  • How to remove the duplicate lines retaining first occurences
    written by oec.rakesh • oec.rakesh
    <p>Let&#39;s say a input text file &quot;input_msg.txt&quot; file ( file size is 70,000 kb ) contains following records..</p> <p>Jan 1 02:32:40 hello welcome to python …
  • assembly of large Genome(tool)
    written by midox midox
    <p>hello,<br /> I wanted to know what is the most appropriate tool assembler for data from a large &quot;large genome&quot; as the plant genome? according …
  • aligner for CORRECTED pacbio long reads
    written by cmo cmo
    <p><em><strong>After </strong></em>error-correcting PacBio long reads using Illumina short-reads, what aligners are adept for aligning the <strong>corrected </strong>PacBio long reads against the genome?</p> <p>Should I think …
  • Genomic Coverage - Samtool's undocumented "depth" verses the poorly documented pileup.
    written by John John
    <p>Hey all :]</p> <p>I use samtools&#39;s depth, and occasionally samtool&#39;s pileup commands to calculate coverage of my reads to the genome before binning for coverage. …