Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • very small N50 values from preliminary assembly of bacterial genome

    Greetings, I am trying to assemble a small microbial genome that is approximately 1.6 Mbp in size. We used Illumina HiSeq technology with 72bp paired end reads. As a first pass at assembly, i used the velvet package. As well, i read through a few tutorials on pre-processing of Illumina data.
    Downloading and installation was quite simple.

    I initially followed nick loman's suggestions for pre-processing sequences.


    Overall, our read qualities displayed median Q-values of 38, with slight decreases towards the 5' prime ends of the sequences.

    However, I found it quite strange that I was not seeing any contigs greater than 60 bp, and the N50 values was around 24. Below are the initial velveth and velvetg commands I used.

    tonybert$
    velveth run1velveth_01022012/ 31 -fasta -shortPaired COLLAPSED.fasta

    tonybert$
    velvetg run2velveth_01022012/ -ins_length 300 -exp_cov 227

    I assessed the contigs.fa file, and found only sequences of ~60bp. Quite disappointing.

    Since then, I have tried using unfiltered reads raw reads, both paired and and single end, and I continually am getting the same results.

    As well, i have tried adjusting the kmer length to 21, as well as running the data through with with no insert length estimation or expected coverage.

    If anyone has any ideas of what the issues might be, i would sincerely appreciate any comments.

  • #2
    If there were problems with your run, and you have the budget to resequence with pacbio (probably < $1k), it may make more sense to simply send some DNA to Expression Analysis. With the new consensus calling module (quiver), the consensus accuracy can be higher than Sanger.

    Comment


    • #3
      being new to illumina, i don't know how exactly to guide whether the run was quality or not. I assume there must be some files in the output from the sequencing center that would have this. Are there any specific characteristics you look for to assess run quality? To be honest, i thought everything went fine; most of the read positions had really high median Q values, i was under the impression this was OK for assessing run quality.

      Thanks again for the comments. -tony

      Comment


      • #4
        Did you try running Velvetg with "-exp_cov auto" instead of "-exp_cov 227"?

        Comment


        • #5
          Where did you get COLLAPSED.fasta from and are you sure that it's properly formatted (i.e. are the read pairs properly interleaved and read 2 is in the correct orientation)? My initial guess would be that there is a problem with your input file such that reads aren't correctly paired or in the right orientation and that is what's causing your bad assemblies.

          My suggestion would be to take your raw read files and run them through without any pre-processing. The pre-processing scripts, if run on read 1 and read 2 separately, most likely resulted in one read of a pair being thrown away at any given step, so when you created your merged file for input into Velvet, the reads were no longer properly paired together, which resulted in a very poor assembly. Not doing any pre-processing will mean that your read pairs stay together correctly and will give you a starting point for what you should expect in terms of # of contigs and N50.

          Once you have the raw assembly done, you can mess around with various pre-processing options, taking care to keep the read files properly paired, and then redo your assembly.


          And if you think you really need to do more sequencing, don't even think about a PacBio run. If you already have an Illumina library made, then find someone with a MiSeq to do a 2x250 PE run for you.

          Comment


          • #6
            That said, mcnelson is right in that you should fork your workflow and attempt to use your paired end files without pre-processing in parallel with testing pre-processing. If you're trying to reduce runtime, I've seen decreases in runtime using velveth with -create_binary, and if you're going to be running velvet and tuning your assembly, then any decrease in time and resources consumed will be helpful.
            Last edited by winsettz; 01-09-2013, 11:22 AM.

            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, 07-10-2024, 07:30 AM
            0 responses
            29 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 07-03-2024, 09:45 AM
            0 responses
            201 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 07-03-2024, 08:54 AM
            0 responses
            212 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 07-02-2024, 03:00 PM
            0 responses
            193 views
            0 likes
            Last Post seqadmin  
            Working...
            X