Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • houyong
    Junior Member
    • Apr 2010
    • 2

    question about perl

    I am a new comer in this field.recently,i have to finish a bioinformatics task.the contents are as follow:
    basic data(from IIlumina GA II or GA IIx):
    1.long pair reads form sequencing:1.fq and 2.fq;
    2.the scaffold constructed by the reads;
    task:
    after the process of sequencing a long pair fragment,using perl to assess the distribution of the length of mate-pair inserts,and find out how many sequences can be assembled at last.
    Can somebody explain what the task want me to do?As a beginer,I really don't know what the length of mate-pair inserts mean?and How to identify the sequences can be used to assembled.
    Thanks very much!
  • simonandrews
    Simon Andrews
    • May 2009
    • 870

    #2
    With Illumina paired end sequencing you produce two short (normally somewhere between 36-100bp) sequence tags which come from either end of an original sequence fragment in your library, and which sit on opposing strands.

    The length of inserts in your library can vary but would normally be between 100-500bp. I guess therefore that the question you're being asked is to figure out what the actual distribution of fragment sizes in the original library was.

    Normally you'd do this by mapping the reads back against a known reference sequence and seeing how far apart they mapped. In your case you say you have a scaffold constructed by the reads which suggest you've assembled them. You could use this assembly to see how far apart the sequence pairs assembled. For this to work though you'd have needed to have assembled your sequences already. I can't think of any way of figuring out the insert size distribution from just a set of fastq files without using a reference or building an assembly.

    Comment

    • houyong
      Junior Member
      • Apr 2010
      • 2

      #3
      thanks very much

      Originally posted by simonandrews View Post
      With Illumina paired end sequencing you produce two short (normally somewhere between 36-100bp) sequence tags which come from either end of an original sequence fragment in your library, and which sit on opposing strands.

      The length of inserts in your library can vary but would normally be between 100-500bp. I guess therefore that the question you're being asked is to figure out what the actual distribution of fragment sizes in the original library was.

      Normally you'd do this by mapping the reads back against a known reference sequence and seeing how far apart they mapped. In your case you say you have a scaffold constructed by the reads which suggest you've assembled them. You could use this assembly to see how far apart the sequence pairs assembled. For this to work though you'd have needed to have assembled your sequences already. I can't think of any way of figuring out the insert size distribution from just a set of fastq files without using a reference or building an assembly.
      thank you very much! I've got it.Just aligne the 1.fq and 2.fq to the reference sequence(the assembled scaffold).then calculate the distance between the two pair-reads.and about the second question,may i have to align the reads each other to remove the dupicated one?

      Comment

      Latest Articles

      Collapse

      • seqadmin
        New Genomics Tools and Methods Shared at AGBT 2025
        by seqadmin


        This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

        The Headliner
        The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
        03-03-2025, 01:39 PM
      • seqadmin
        Investigating the Gut Microbiome Through Diet and Spatial Biology
        by seqadmin




        The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
        02-24-2025, 06:31 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 03-20-2025, 05:03 AM
      0 responses
      17 views
      0 reactions
      Last Post seqadmin  
      Started by seqadmin, 03-19-2025, 07:27 AM
      0 responses
      18 views
      0 reactions
      Last Post seqadmin  
      Started by seqadmin, 03-18-2025, 12:50 PM
      0 responses
      19 views
      0 reactions
      Last Post seqadmin  
      Started by seqadmin, 03-03-2025, 01:15 PM
      0 responses
      186 views
      0 reactions
      Last Post seqadmin  
      Working...