Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Increasing contig lengths

    I'm working on a project to identify by sequence viruses infecting grapevines.

    I have single end Illumina reads (50 bp) and have been trying to assemble them using a combination of Velvet and PRICE. I've been able to get to a max contig length of around 1500 with Velvet and an n50 of 46. After putting this output through PRICE, I can increase the n50 to 195. However, I am having trouble increasing my contig length after this. Do you have any advice regarding contig extension with single end reads?

  • #2
    Have you tried varying the kmer length when assembling? Also, it would be helpful to know more about your data, like the read length and total amount, and quality metrics.

    I encourage you to read this thread:
    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

    Comment


    • #3
      ... and the amount of contaminating sequences?

      Originally posted by Brian Bushnell View Post
      Have you tried varying the kmer length when assembling? Also, it would be helpful to know more about your data, like the read length and total amount, and quality metrics.

      I encourage you to read this thread:
      http://seqanswers.com/forums/showthread.php?t=42555

      Comment


      • #4
        Thanks for the replies.

        I used VelvetOptimiser to determine optimal k-mer length. Our data contains a mixture of grape and virus reads, but we removed the reads that aligned to the grape reference genome. Our read length is 50 bp and we have 7,764,190 reads after filtering out the grape reads.

        Here is the quast output from the optimal velvet run:

        All statistics are based on contigs of size >= 100 bp, unless otherwise noted (e.g., "# contigs (>= 0 bp)" and "Total length (>= 0 bp)" include all contigs).

        Assembly contigs
        # contigs (>= 0 bp) 3547
        # contigs (>= 1000 bp) 1
        Total length (>= 0 bp) 326445
        Total length (>= 1000 bp) 1073
        # contigs 941
        Largest contig 1073
        Total length 156559
        GC (%) 46.84
        N50 162
        N75 122
        L50 305
        L75 584
        # N's per 100 kbp 0.00

        Comment


        • #5
          Are you trying to assemble genomic data or transcriptomic data?

          What is the expected genome size of the virus genome you are trying to assemble?

          What kmer length have you used?

          As Brian already mentioned above, I would play around with the kmer length
          when using velvet, to see what kmer length gives you the best n50.

          Have you done any QC, adapter trimming or quality trimming on your reads?

          Comment


          • #6
            Do you think the viral genome will be divergent within a sample from replication errors? That could cause issues for assembly if there are lots of related kmers at a location instead of just one or two alleles and a low level of sequencing error.
            Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

            Comment


            • #7
              I was able to assemble my data using IDBA_UD. I set it to cycle through k mers that were less than my read size and it produced a 15000bp sequence: idba_ud -r ../data/trimmed-reads/LV89-02.fa -o ../results/contigs/008 --mink 19 --maxk 49 --step 2

              Comment


              • #8
                quast quality stats of the assembly?

                Comment


                • #9
                  All statistics are based on contigs of size >= 100 bp, unless otherwise noted (e.g., "# contigs (>= 0 bp)" and "Total length (>= 0 bp)" include all contigs).

                  Assembly contig
                  # contigs (>= 0 bp) 295
                  # contigs (>= 1000 bp) 37
                  Total length (>= 0 bp) 199556
                  Total length (>= 1000 bp) 86293
                  # contigs 295
                  Largest contig 15124
                  Total length 199556
                  GC (%) 46.27
                  N50 759
                  N75 426
                  L50 53
                  L75 141
                  # N's per 100 kbp 0.00

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Exploring the Dynamics of the Tumor Microenvironment
                    by seqadmin




                    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                    07-08-2024, 03:19 PM
                  • seqadmin
                    Exploring Human Diversity Through Large-Scale Omics
                    by seqadmin


                    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                    06-25-2024, 06:43 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Today, 06:53 AM
                  0 responses
                  4 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-10-2024, 07:30 AM
                  0 responses
                  30 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-03-2024, 09:45 AM
                  0 responses
                  202 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 07-03-2024, 08:54 AM
                  0 responses
                  212 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X