Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • parameters for Denovo assembly

    for bacterial genome

    how the parameters bubble size , K-mer , Length size , contig length , mimatch cost , insertion cost , deletion cost , length fraction , similarity fraction effects for Illumina data in Denovo Assembly

  • #2
    The documentation of the program should explain each of these parameters in great detail.

    Comment


    • #3
      Guys, I am new to NGS and have a problem. I have been provided with fastq sequences of an yeast which are 66 in numbered and of the order abcd_2_NoIndex_L002_R1_001.fastq upto abcd_2_NoIndex_L002_R33_001.fastq and similarly abcd_2_NoIndex_L002_R2_001.fastq upto abcd_2_NoIndex_L002_R2_033.fastq and have been told it's a paired end data. How do i carry out de novo assembly of the same as i am interested in finding out the presence of certain genes of interest. I had been told that velvet is one such assembler for these, but it accepts only one pair of r interleaved files for processing. How do i modify my files according to velvet standards. Kindly guide me on the commands and steps required for the same. Thanks a ton in advance.

      Comment


      • #4
        diptarka,
        Within your velvet directory there should be a perl script called shuffleSequences_fasta.pl

        This will order your 2 fastq files into a single, correctly ordered file.

        Usage From Velvet Manual:
        ./shuffleSequences_fasta.pl forward_reads.fa reverse_reads.fa output.fa

        Comment


        • #5
          @jgibbons1:
          actually i have in total 66 fastq files. i know the shufflesequences.pl script. what i am asking is the 33 fastq files for the forward end that i have, can i combine them all into 1? and so for the rest 33 for paired end?otherwise the shufflesequences.pl script won't run. I was asking whether the procedure is correct or not? cos, i just don't have two fastq files, i have 33+33 fastq files that are paired end.

          Comment


          • #6
            Ahhh...I see, sorry I miss understood.

            It sounds like you have 2 options:

            (1) Use the shufflesequences.pl script on each of the 33 pairs
            Velvet can handle multiple input files, so you can list all files in velveth, or simply use a wildcard -shortPaired *.fastq

            (2) Merge the 33 forward and 33 reverse fastq files and then run the shufflesequences.pl script

            You can merge files in a unix environment using the 'cat' command. If you go with this route, be sure to keep the files in the same order.

            Comment


            • #7
              @jgibbons1:
              i can merge the 33 files into one single abcd_1.fastq and similarly the rest into abcd_2.fastq. But, the read info won't be hampered is it?I have actually done that, but was unsure whether further processing via this method will hamper or not.

              Comment


              • #8
                Merging the files shouldn't change any information in the files themselves. The key part is making sure the files in both reads are merged in the same order to maintain their pairing.

                Comment


                • #9
                  by same order you mean sequentially right from 1 to 33 for one end and 1 to 33 for the other. I guess the cat command in linux will merge sequentially.

                  Comment


                  • #10
                    Do, any filtering of the reads necessary before the files are merged and interleaved for velvet to be run?any qc checks required?

                    Comment


                    • #11
                      Exactly what kcchan said. Just make sure you merge the files in the same order (1, 2, 3...33) in both read sets.

                      Filtering is a different story. If you are trimming the file as a whole, (i.e. trim 10 bp from beginning and end) there should be no problem. If you are filtering reads individually you need to be careful that bad reads are removed from both files. A quick way to check if this worked correctly is by checking the number of reads in each file at the end of filtering:

                      wc -l yourfile.fastq | awk '{seq=$1/4} END {print seq}'

                      Comment

                      Latest Articles

                      Collapse

                      • seqadmin
                        Strategies for Sequencing Challenging Samples
                        by seqadmin


                        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                        03-22-2024, 06:39 AM
                      • seqadmin
                        Techniques and Challenges in Conservation Genomics
                        by seqadmin



                        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                        Avian Conservation
                        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                        03-08-2024, 10:41 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, 03-27-2024, 06:37 PM
                      0 responses
                      12 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-27-2024, 06:07 PM
                      0 responses
                      11 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-22-2024, 10:03 AM
                      0 responses
                      53 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-21-2024, 07:32 AM
                      0 responses
                      68 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X