Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Effect
    Junior Member
    • Feb 2013
    • 2

    #1

    Build VCF for exons

    I have long reads for human exome and also I have reference exome. This reference isn't a whole human genome but only a collection of exons, totally about 400000 sequences.
    So reference file looks like:
    Code:
    >ENSG00000187_DAZ3_chrY_26920497-26920568_exon
    AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
    >ENSG00000187_DAZ3_chrY_26930284-26930355_exon
    CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
    >ENSG00000187_DAZ3_chrY_26932661-26932732_exon
    GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG
    I want to build VCF file for my reads and this reference.
    I wonder is it possible to work with such multi-sequence reference? Or should I megre sequences from one chromosome together and then work with this file?

    I've tried use bwa for build .bam file at first stage and faced to the trouble I can't browse the result .bam file. Samtools -tview shows me only first sequence of reference and then put NNNNN letters. Integrative Genomics Viewer also shows me nothing for my .bam file.
  • swbarnes2
    Senior Member
    • May 2008
    • 910

    #2
    I think you want to be aligning to genome, not a bunch of exonic sequence. You did some kind of exon capture? I which case, you'll have contaminating genomic DNA, and you want that to align to where it belongs. Also, you likely have some reads that hang off of the end of the exon, and those won't align to the right exon without enough flanking genomic sequence.

    Comment

    • Effect
      Junior Member
      • Feb 2013
      • 2

      #3
      Thank you, for whole genome sequence it works good.

      Comment

      Latest Articles

      Collapse

      • kinpict.com
        Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by kinpict.com
        The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
        Today, 06:37 PM
      • SEQadmin2
        16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by SEQadmin2





        NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
        10-30-2026, 06:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 11-12-2026, 07:59 AM
      0 responses
      27 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 11-05-2026, 09:25 AM
      0 responses
      35 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 11-03-2026, 11:54 AM
      2 responses
      56 views
      0 reactions
      Last Post kinpict.com  
      Started by SEQadmin2, 10-28-2026, 09:57 AM
      0 responses
      175 views
      0 reactions
      Last Post SEQadmin2  
      Working...