Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • laura
    replied
    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/

    Leave a comment:


  • nexgengirl
    replied
    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!

    Leave a comment:


  • Alex Renwick
    replied
    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

    Leave a comment:


  • swbarnes2
    replied
    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.

    Leave a comment:


  • mez
    started a topic 1000 genomes site BAM Files

    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.

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
26 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