Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 1000 genomes site BAM Files

    Please,

    I would like anyone who has had experience working with the .bam alignment files from the 1000 genomes site to help me out with this situation. I have been trying for long to get a complete FASTA sequence of the sample DNA from the .bam files for a given genomic region. I had used the samtools view, sort, index as well as pile up commands to get a pile up file and eventually extracted the consensus DNA sequence from the fourth column of the pile up file, creating the FASTA file.

    The problem I am having with this approach, is that the raw DNA sequence seem to be obviously shorter in each case than the DNA sequence from the reference genome, and for different samples I get different DNA sequence length. I do not know if the issue is because of the low coverage sequence approach of the 1000 genomes site, or if I am doing something wrong.

    Thanks and will appreciate any guidance.

  • #2
    You'd think that there woud be an out-of-the-box answer to making fastas from .bams, but there isn't.

    The first thing I see wrong with your approach is that pileup will simply not show any line if it is totally uncovered. Are you padding your fasta with N's when there is no corresponding line in the pileup?

    Older versions of samtools came with a samtools.pl that had a utility to make consensus fastq of pileups. You could try using that. perl scripts are flat text files, you can read them like an ordinary text file.

    Search around seqanswers a bit, I know I've written previously on how to convert a perl script that writes a fastq into one that writes a fasta. It's just a matter of changing a few characters in the script.

    Comment


    • #3
      Originally posted by mez View Post
      ... extracted the consensus DNA sequence from the fourth column of the pile up file, creating the FASTA file.
      ...
      By "fourth column" I suppose you mean "third column", since the fourth column of the pileup format is the count of reads covering the site. See pileup format.

      The sequence you get from the third column is the reference sequence, not a consensus sequence. It will be shorter because sites with no read coverage won't appear.

      See http://samtools.sourceforge.net/mpileup.shtml for instructions on using samtools to generate a consensus sequence:

      Code:
      samtools mpileup -uf ref.fa aln.bam | bcftools view -cg - | vcfutils.pl vcf2fq > cns.fq

      Comment


      • #4
        Check out this tool:



        It will convert a sam or a bam to a fastq. That would get you at least half way to a fasta file.

        I'm pretty sure there have been other posts on how to convert from fastq to fasta+qual so I would check on seqanswers for more about that.

        Good luck!

        Comment


        • #5
          If you are after all the sequence for a particular individual from the 1000 genomes project you are much better using our fastq files rather than our bam files

          ftp://ftp.1000genomes.ebi.ac.uk/vol1....sequence_data
          ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/sequence.index

          If you want just the sequence which has been aligned to a particular region then a solution like picards SamtoFastq which has already been suggested would be best

          If you want the reference genome annotated for all the variants in a particular individual then you should look at solutions like



          or

          vcfutils.pl vcf2fq from samtools

          amd our vcf files which you can get from

          ftp://ftp.1000genomes.ebi.ac.uk/vol1...ease/20110521/

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Recent Developments in Metagenomics
            by seqadmin





            Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
            09-23-2024, 06:35 AM
          • seqadmin
            Understanding Genetic Influence on Infectious Disease
            by seqadmin




            During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

            Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
            09-09-2024, 10:59 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 04:51 AM
          0 responses
          8 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-01-2024, 07:10 AM
          0 responses
          12 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 09-30-2024, 08:33 AM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 09-26-2024, 12:57 PM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Working...
          X