Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Does BAM file from SRA contain aligned reads?

    Hello, everyone. I need to do SNP calling for the first time, so I have a couple questions.
    For example, I want to use this sample https://trace.ncbi.nlm.nih.gov/Trace...run=SRR1295568

    I can download from there reads in FASTQ format, align them on the reference genome using BWA, for example. After that I will get SAM file, then BAM file and so on.

    But also, I can download BAM file from SRA(you choose this option on the bottom line).
    As I understand, this file already contains aligned on the reference genome reads. So you don't need to align them, this work is already done, just launch SAMTools, bcftools and do SNPcalling.

    So, my questions:
    1) Am I right?
    2) If Yes, is it better to align reads by yourself, using BWA, or just download BAM file? Would be there some differences in result?

  • #2
    2) If Yes, is it better to align reads by yourself, using BWA, or just download BAM file? Would be there some differences in result?
    That depends on your experience, available time, and domain knowledge.

    I generally prefer doing my own read cleaning and alignment, because then I know precisely the process that was used for read filtering and mapping. From my own memory of how submissions work, I think the BAM files in SRA are aligned using whatever method the experimenter decided to use.

    Comment


    • #3
      As mentioned, the BAM files you download were prepared and uploaded by the submitter (FASTQ, while supported, is not a preferred format for SRA and is accepted as reference sequences for the aligned read files). So, if you read their methods and are happy with their QC cutoffs and mapping approach, you could use their BAM files. If not, you need to get the actual read files and do your own reference mapping.

      Also take note of how old any submission is and what reference release they used. You may want to re-align the reads to a more recent reference build anyway.

      A lot of times, yes, it would make little if any difference how you proceed. Sometimes though, it might. Look through the information provided about the submission and decide for yourself how you wish to proceed.
      Michael Black, Ph.D.
      ScitoVation LLC. RTP, N.C.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Essential Discoveries and Tools in Epitranscriptomics
        by seqadmin


        The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
        Yesterday, 07:01 AM
      • seqadmin
        Current Approaches to Protein Sequencing
        by seqadmin


        Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
        04-04-2024, 04:25 PM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 04-11-2024, 12:08 PM
      0 responses
      37 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-10-2024, 10:19 PM
      0 responses
      41 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-10-2024, 09:21 AM
      0 responses
      35 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-04-2024, 09:00 AM
      0 responses
      55 views
      0 likes
      Last Post seqadmin  
      Working...
      X