Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • VCF calling for individual chromosomes, domain reduction

    Hello.

    I am just seeing if this method for VCF creating for variant calls is correct.
    Right now I keep the fastq as an entire file, and then align the whole fastq file to individual chromosomes, and this produces a VCF file for each chromosome that I later merge together.

    However this can cause errors if the alignment has a poor map to the individual chromosome because it will generate a forced alignment, even if it maps poorly to the chromosome. (ie the VCF call will show a score to the poor map of the fastq to the chromosome if there is only a single chromosome the fastq is being aligned to:: (ie. reducing the domain of the fastq to a single chromsome at a time)


    Is it better not to reduce the alignment map to a single chromosome at a time and instead map to the whole chromosome at once?

  • #2
    Yes, it is much better to align to the entire genome in a single step. If you want to reduce the memory load then separate the fastq into a few smaller files and align them separately instead of separating the genome like that.

    Comment


    • #3
      I am not sure how to do this.

      how do you split the fastq files? how to align the whole genome?

      when I align the fastq using Tophat, I align both reads at the same time.

      Comment


      • #4
        A fastq is just a text file with one read every 4 lines, if you have paired end data it then two text files. The reads should be in the same order in both files by default so just split the files into smaller files making sure that the number of lines is a factor of 4. Actually, splitting probably isn't necessary anyway.

        Comment


        • #5
          Thank you
          I am familiar with the alignment practice, and am interested in the variant calls.
          How do I create the VCF files, without doing the variant calling for a single chromosome at a time?

          Comment


          • #6
            Not much information to go on there. But the first step is to map against the whole genome not single chromosomes, so if you are using tophat then make a bowtie reference of the whole genome, map to that then use the resulting bam to call variants using whatever program you use which will call variants for the whole genome.

            Comment


            • #7
              Vcf

              yes. this is the protocol I follow.

              I do have an HG19 human reference genome that I align the fastq files to.
              after I use same tools to index the accepted hits.bam file and use mpileup -g command to call the variants.

              what parameters do you use to call the variants? do you use mpileup -g ?

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Genetic Variation in Immunogenetics and Antibody Diversity
                by seqadmin



                The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
                11-06-2024, 07:24 PM
              • seqadmin
                Choosing Between NGS and qPCR
                by seqadmin



                Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                10-18-2024, 07:11 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 11-08-2024, 11:09 AM
              0 responses
              42 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 11-08-2024, 06:13 AM
              0 responses
              31 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 11-01-2024, 06:09 AM
              0 responses
              34 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-30-2024, 05:31 AM
              0 responses
              23 views
              0 likes
              Last Post seqadmin  
              Working...
              X