Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Makoe
    Junior Member
    • Jul 2014
    • 2

    Masking A/C/G/T with N?

    I'm working with a tetraploid soybean species, and I want to determine which DNA-seq reads exist on only one of the two diploid progenitor genomes. I have two consensus fasta files (one for each diploid genome) and want to map reads to them to determine which reads map to only one genome based on mapping scores. I already have a pipeline for accomplishing this, but I'm running into a problem. One consensus file has 30 million more Ns than the other. Reads that should map equally to both species favor one because it has fewer Ns, and therefore a higher mapping score. Is there a way to replace A/C/G/T's in one consensus file with an N where the other consensus file has an N? The consensus sequences were made by mapping RNA-seq and DNA-seq reads to the Glycine max genome, so they should be comparable.
  • N311V
    Member
    • Jul 2013
    • 34

    #2
    I'm not aware of a ready-made tool to do this job but it sounds like something you could do with Perl. I'm sure there are online tutorials you could use otherwise I found this book useful when first starting (http://shop.oreilly.com/product/9780596000806.do)

    Are both consensus sequences the same length for each chromosome/scaffold? Doesn't really matter if they aren't, would just make the perl script to do this job easier.

    Comment

    • Makoe
      Junior Member
      • Jul 2014
      • 2

      #3
      Ah, I guess I'll just have to build up my Perl skills. Probably a good idea, anyway. And yes, the consensus sequences are the same size, so it shouldn't be too difficult.

      Thanks for your help!

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Pathogen Surveillance with Advanced Genomic Tools
        by seqadmin




        The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
        Yesterday, 11:48 AM
      • 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
      26 views
      0 reactions
      Last Post seqadmin  
      Started by seqadmin, 03-19-2025, 07:27 AM
      0 responses
      33 views
      0 reactions
      Last Post seqadmin  
      Started by seqadmin, 03-18-2025, 12:50 PM
      0 responses
      25 views
      0 reactions
      Last Post seqadmin  
      Started by seqadmin, 03-03-2025, 01:15 PM
      0 responses
      190 views
      0 reactions
      Last Post seqadmin  
      Working...