Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • consensus sequence

    have mapped illumina PE reads to the reference genome using BWA and now, would like to get a consensus sequence of the mapped reads- can samtools be used on the aligned bam file??

  • #2
    See this thread: http://seqanswers.com/forums/showthread.php?t=14565 (and the other thread mentioned in #4 there).

    Comment


    • #3
      samtools mpileup -uf /<reference.fasta> <aligned.bam> | bcftools view -bcg - > <intermediate.bcf>

      gave me an error: -bash: syntax error near unexpected token "<

      I was wondering BCFtools are part of samtools, right?- is this the correct comand?

      Comment


      • #4
        You don't literally use < > around the file names.

        Comment


        • #5
          has anyone tried this command to generate a consensus file

          samtools mpileup -uf reference.fasta sorted_alignment.bam | bcftools view -cg - | vcfutils.pl vcf2fq

          but where should I put the output file name???

          Comment


          • #6
            Originally posted by mmmm View Post

            but where should I put the output file name???
            Code:
            $ samtools mpileup -uf reference.fasta sorted_alignment.bam | bcftools view -cg - | vcfutils.pl vcf2fq > output_file_name

            Comment


            • #7
              yes- thanks.
              I will ignore the abnormal shapes that appeared after the sequence
              ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

              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, Yesterday, 06:37 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              51 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