Announcement

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

  • figo1019
    replied
    Originally posted by Torst View Post
    The issue seems to be in the "bcftools" step in your pipe ie. [bcf_sync].

    Googling suggests two things:

    1. there was an bug in previous bcftools when using "-" stdin, so maybe you have an old bcftools.

    2. this occurs when you have indels, so try the mpileup "-I" option to not call indels and see if the error still occurs.
    Dear Torst

    Thank you for your reply.

    1. I am using the samtools 0.1.18 version.
    2. I tried what you suggested and I didn't get the error this time but it again stops at the same point and I miss the same no. of chromosomes in the consensus sequence file e.g. when use grep to count the no. of chr in the .fq file I find the few chromosomes missing which was the same with the error also.

    Regards

    Leave a comment:


  • Torst
    replied
    The issue seems to be in the "bcftools" step in your pipe ie. [bcf_sync].

    Googling suggests two things:

    1. there was an bug in previous bcftools when using "-" stdin, so maybe you have an old bcftools.

    2. this occurs when you have indels, so try the mpileup "-I" option to not call indels and see if the error still occurs.

    Leave a comment:


  • figo1019
    started a topic Samtools consensus sequence error

    Samtools consensus sequence error

    Dear All

    I am again posting this question to the Bioinformatics forum as I didn't get the answer from the Illumina forum.

    I am trying to make consensus sequences using command
    samtools mpileup -uf ref.fa aln.bam | bcftools view -cg - | vcfutils.pl vcf2fq > cns.fq.
    But i am getting this error:

    [bcf_sync] incorrect number of fields (4 != 5) at 15:16294867
    [afs] 0:49470.732 1:50.052 2:33.216

    Can some one help me to sort out this problem. What exactly is going wrong.

    Regards

Latest Articles

Collapse

  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
30 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
18 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-14-2024, 09:19 AM
0 responses
52 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-12-2024, 03:37 PM
0 responses
451 views
0 likes
Last Post seqadmin  
Working...
X