Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • metrics on BAM file

    Hi everyone,
    I am doing exome sequencing, paired end, but not used to analyze reads.
    After alignement, how do i proceed the BAM file?
    Thank you

  • #2
    How many samples do you have?
    How did you align the reads?
    Are they tumor and normal samples?
    Are looking for germline characteristics or somatic changes?

    Comment


    • #3
      I have three normal samples and looking for germline characteristics.
      Mapping with BWA-MEM.
      Thank you for you help Richard

      Comment


      • #4
        You will probably want to do some variant calling on the bam files. This may include SNPS, fusions, amplification or deletions.

        There are many tools to do this. Searching for these tools is straightforward:



        You may want to use favorite, known tools already used at your organization, so ask around. A new tool may require some non-trivial installation.

        Comment


        • #5
          Thank you for you reply.
          i need to understand some details, once i have a BAM file, before variant calling i have to remove duplicate, recalibrate BAM? run other metrics ?
          sorry i am completely new to this.

          Comment


          • #6
            You could do some quality checks and remove duplicates and recalibrate; but if you are just starting off, just go ahead and assume the data is good and try and see if you can get a report on variations using a commonly used tool. After you get that working, you can check out the other stuff and re-run a new workflow. Order of importance would be : QC, duplicate removal or marking , then re calibration. Many here at work say recalibration is not worth the effort. Theoretically marking or removing duplicates may be important upstream to a variant caller, and again, theoretically a variant caller could do that work for you so read the docs for your tools and follow their advice.

            Do try and run a QC tool on your fastqs early.
            Absolute easiest QC on bam files is "samtools flagstat" (google it).
            Study the results and see if you can make sense of them.

            Comment


            • #7
              Thank you, i ll give a try. Do i ignore UMi for this analysis?
              Last edited by Moip; 09-28-2018, 10:34 PM.

              Comment

              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-19-2024, 07:20 AM
              0 responses
              40 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-16-2024, 05:49 AM
              0 responses
              52 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-15-2024, 06:53 AM
              0 responses
              63 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-10-2024, 07:30 AM
              0 responses
              43 views
              0 likes
              Last Post seqadmin  
              Working...
              X