Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Mithril
    Member
    • Apr 2012
    • 12

    SNP detection by samtools returns dots under the reference columns

    Hello,
    I have a very peculiar problem. I used samtools to get SNPs using the mpileup command. The commands that i used are as follows:

    samtools mpileup -uf cpp34_contigs_nr_prinseq_good_NAF4.fasta aln-sorted.bam | bcftools view -bvcg - > var.raw.bc

    bcftools view var.raw.bcf > var.flt.vcf

    The .vcf file that i got had the follwing results:
    #CHROM POS ID REF ALT QUAL
    PT_34_209 159 . T C 7.8
    PT_34_209 164 . C T 7.8
    PT_34_209 655 . A G 7.8
    PT_34_209 658 . G C 7.8
    PT_34_209 659 . C G 7.8
    PT_34_238 533 . T C 7.8
    PT_34_238 600 . T C 7.8
    PT_34_238 786 . A T 7.8
    PT_34_238 789 . A C 7.8
    PT_34_238 790 . T C 7.8
    PT_34_238 791 . T G 7.8
    PT_34_238 792 . T C 7.8
    PT_34_238 793 . G T 7.8
    PT_34_238 795 . T G 7.8
    PT_34_258 561 . A G 7.8
    PT_34_549 520 . T C 7.8

    Could anybody please help me with this? and let me know how i can correct the problem between the reference column and the quality column?
  • Alex Renwick
    Member
    • Jul 2011
    • 44

    #2
    I don't see any problem here. For example, the first row of you results says chromosome PT_34_209 has reference allele T at position 159 and a C snp was observed with quality 7.8, and there is no snp id for that position. What are you expecting to see that is not here?

    Comment

    • Mithril
      Member
      • Apr 2012
      • 12

      #3
      Ow... My bad. i just realised my mistake. thank you
      Last edited by Mithril; 04-13-2012, 06:49 AM.

      Comment

      • Alex Renwick
        Member
        • Jul 2011
        • 44

        #4
        Again I'll say that there is no problem. In the vcf format the entries are separated by tabs. A long entry (the chromosome, in this case) can make the columns look misaligned, but the downstream analysis tools will have no trouble parsing the file.

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
          by SEQadmin2


          Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


          The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
          ...
          06-02-2026, 10:05 AM
        • SEQadmin2
          Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends
          by SEQadmin2


          With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.


          Introduction

          Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
          05-22-2026, 06:42 AM
        • SEQadmin2
          Environmental Genomics in the Age of NGS: From Microbes to Conservation Strategies
          by SEQadmin2

          Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.


          Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...
          05-06-2026, 09:04 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Yesterday, 08:59 AM
        0 responses
        13 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 06-02-2026, 12:03 PM
        0 responses
        22 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 06-02-2026, 11:40 AM
        0 responses
        19 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 05-28-2026, 11:40 AM
        0 responses
        31 views
        0 reactions
        Last Post SEQadmin2  
        Working...