Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bam2mpg output

    Dear all,
    I'm trying BAM2MPG software to analyze the goodness of my variant call.
    But I can't understand the output when i run
    > bam2mpg hg19.fasta sample.bam

    This is the output:
    MPG_SNV chr1 824017 g AA 0 1 1
    MPG_SNV chr1 835499 A GG 0 1 1
    MPG_SNV chr1 846489 T CC 0 1 1
    MPG_SNV chr1 852758 G CC 0 1 1
    MPG_SNV chr1 862093 T CC 0 1 1
    MPG_SNV chr1 862124 A GG 0 1 1
    MPG_SNV chr1 862383 C TT 0 1 1
    MPG_SNV chr1 862389 A GG 0 1 1
    MPG_SNV chr1 865219 G AA 0 1 1
    MPG_SNV chr1 869303 C TT 0 1 1
    MPG_SNV chr1 869323 T CC 0 1 1
    MPG_SNV chr1 869369 G TT 0 1 1
    MPG_SNV chr1 869373 G CC 0 1 1
    MPG_SNV chr1 869401 G CC 0 1 1
    MPG_SNV chr1 879563 A GG 0 1 1
    MPG_SNV chr1 895706 G AA 0 1 1

    Anyone could explain me what are the values in last three columns?
    Thankx a lot
    ME

  • #2
    Hi ME--

    The final three columns are:

    genotype score: a reasonable cutoff for a high quality genotype is 10, corresponding to expected accuracy of roughly 99.9%.

    ref/nonref column: this column is one when the genotype contains an allele that doesn't match the reference (useful for pulling variants with "awk").

    coverage: the number of reads used to call the genotype (if you are using the --qual_filter option, low quality reads won't be counted).

    It looks from your output like you have only one read in this region, and that the bases don't match the reference. My best guess is that you're specifying the wrong reference fasta file. What does this region look like if you view it in samtools tview or IGV?

    --Nancy

    Comment


    • #3
      Thank you Nancy, you are right! I was using the wrong reference fasta file. Now bam2mpg run correctly and your explanation about the column data are very useful for me!

      Comment


      • #4
        Great! I'm sorry I didn't notice your post earlier--glad you're getting good use of the program. There's actually fairly extensive documentation within the perl script, so type "perldoc bam2mpg" to get more information on options and output.

        --Nancy

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 06:37 PM
        0 responses
        11 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        51 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        67 views
        0 likes
        Last Post seqadmin  
        Working...
        X