Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • ulz_peter
    replied
    1) In the directory of the bam files you'll find a .bas file including a MD5 field

    2) that should actually work, providing the header of the first bam file is correct: to be sure try:

    Code:
    samtools view -H in.bam > header.sam
    samtools merge -h header.sam out.bam in1.bam in2.bam ....
    But be sure that the first in.bam is the mapped bam file

    3) Your header does NOT look ok. It should mention every reference sequence so at least 24 entries (including X, Y and MT)...

    Leave a comment:


  • zhanglu295
    replied
    Thank you very much for your help.

    1. Can you tell me where I can find the MD5 code for bam files from 1000 genome website? I can only find the fastq MD5.

    2. I used the samtools merge function to combine this three files(./samtools merge out.bam in1.bam in2.bam in3.bam), is it appropriate here?

    3.The head appeared as
    "@SQ SN:GL000192.1 LN:547496 AS:NCBI37 UR:ftp://ftp.1000genomes.ebi.ac.uk/vol1...k_v37.fasta.gz M5:325ba9e808f669dfeee210fdd7b470ac SP:Human"
    Is it correct?

    Sorry for trouble you so much.
    I do not have much experience in this field.

    Thank you.

    Leave a comment:


  • ulz_peter
    replied
    Did you check the MD5 checksum of your downloaded files. Maybe there was something wrong with the download?

    How did you combine the three files (just using cat or similar tools wouldn't work as all of them contain headers)?

    Did you check the headers of the file. Is every reference sequence in the header section (./samtools view -H)?

    Leave a comment:


  • zhanglu295
    replied
    Yes, I have used the reference sequence from 1000 genome project form the website:

    ftp://ftp.sanger.ac.uk/pub/1000genom...k_v37.fasta.gz

    This time I can successfully call SNPs from some chromosomes but not all the chromosomes.

    And when I used the samtool index function

    such messeage appeared:

    [bam_header_read] EOF marker is absent.
    [bam_index_core] truncated file? Continue anyway. (-2)

    Leave a comment:


  • ulz_peter
    replied
    I assume you already changed the reference sequence file names, did you?

    Leave a comment:


  • zhanglu295
    started a topic 1000 genome data download

    1000 genome data download

    I have downloaded the three bam files from

    ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data/

    for each individual such as

    HG00119.chrom20.ILLUMINA.bwa.GBR.low_coverage.20101123.bam
    HG00119.mapped.ILLUMINA.bwa.GBR.low_coverage.20101123.bam
    HG00119.unmapped.ILLUMINA.bwa.GBR.low_coverage.20101123.bam

    Then I combined them together to use samtools calling SNP.

    But the result (pileup file) did not contain all the chromosomes.

    Is it only a truncated file now?

Latest Articles

Collapse

  • seqadmin
    Quality Control Essentials for Next-Generation Sequencing Workflows
    by seqadmin




    Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

    Nucleic Acid Quality Control
    Preparing for NGS starts with isolating the...
    02-10-2025, 01:58 PM
  • seqadmin
    An Introduction to the Technologies Transforming Precision Medicine
    by seqadmin


    In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
    01-27-2025, 07:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 02-07-2025, 09:30 AM
0 responses
66 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-05-2025, 10:34 AM
0 responses
101 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-03-2025, 09:07 AM
0 responses
81 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-31-2025, 08:31 AM
0 responses
45 views
0 likes
Last Post seqadmin  
Working...
X