Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • someperson
    Junior Member
    • Jul 2013
    • 9

    #1

    Velvet + meta-velvet:which fastq format for paired-end infos?

    Hi,
    I've been doing some metagenomic assemblies using meta-velvet. I've been using paired-end data for mor reliable assemblies. However, when I wanted to check those assemblies (e.g. using tablet) I could not find the paired end info any-more.
    I've recently learned that fastq-formats have changed not only with respect to quality values, but also with respect to paired-end info.

    Old format:
    @HWUSI-EAS100R:6:73:941:1973#0/1
    New Format (that I am using):
    @EAS139:136:FC706VJ:2:2104:15343:197393 1:Y:18:ATCACG
    (paired-end info bold and underlined)

    So Mira for example needs the old format and will just ignore paired-end info it doesn't recognize (without warning!) if you use the new Format.

    Does anybody know which Format Velvet and Meta-Velvet use? Are my assemblies now all in reality single-read assemblies?
  • mastal
    Senior Member
    • Mar 2009
    • 666

    #2
    Velvet + meta-velvet:which fastq format for paired-end infos?

    I have been using velvet with the new formats of MiSeq and HiSeq PE data, and it works fine.

    By the way, in the read that you have used to illustrate the new format, the Y (1:Y:18:ATCACG) indicates a read that has not passed Illumina's QC filter step.

    Usually the sequence providers remove those reads.

    Comment

    • muol
      Member
      • Jun 2012
      • 10

      #3
      Velvet ignores fastq headers and replaces them with its own (fasta) headers anyway. You just have to make sure that input read pairs are either shuffled in one file, or synced if separate files are used.

      Comment

      • someperson
        Junior Member
        • Jul 2013
        • 9

        #4
        @mastal
        Thanks, but the example reads I've posted are not mine. they are from the wikipedia page illustrating the differences between the fastq-formats.
        @mastal+@muol:
        But then how can I check the assemblies with regard to paired-end Infomation? With Mira, for example, I know that you still get assemblies (without error-messages) if you use the wrong format. The assemblies will just in reality be based on single-read assemblies. So how can one really make sure if everything is ok, if the paired-end info is left out from the AMOS-file, the sequence-list and everything? (This just intrigues me. For some contigs i would simply like to check how well supported they are by paired reads. Also a lot of reads are not assembled or dropped from the assembly. I would like to see what has happened to their corresponding partners)
        Last edited by someperson; 07-04-2013, 08:29 AM.

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
          by SEQadmin2



          Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

          Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
          ...
          10-14-2026, 09:55 AM
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 10-22-2026, 11:13 AM
        0 responses
        40 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-22-2026, 10:48 AM
        0 responses
        29 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        41 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        45 views
        0 reactions
        Last Post SEQadmin2  
        Working...