Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • IGV-Compatible FASTA Nomenclature to Indicate a Split Sequence?

    Hi! New poster here, so please tell me if I should make any changes in post style or should have used a different subforum!

    This is a x-post from StackBioinformatics, but I'm very interested in an answer and so I figured this might be the best place.

    I have a nucleotide sequence where I am interested in aligning reads distinctly to one or another portion of the sequence. However, I would like to use a visualization tool such as IGV to view these alignments in the context of the parent sequence. For example, I would like to split something along the lines of:

    >seq1 parent_sequence
    ACTGACTGACTGACTGACTGGTCAGTCAGTCAGTCAGTCAACACACACACACACACACAC

    Into:

    >seq1_1:20 parent_sequence subsequence_1_to_20
    ACTGACTGACTGACTGACTG
    >seq1_21:40 parent_sequence subsequence_21_to_40
    GTCAGTCAGTCAGTCAGTCA
    >seq1_41:60 parent_sequence subsequence_41_to_60
    ACACACACACACACACACAC

    Then perform alignments using the split subsequences as a reference for alignment of reads, but visualize these read alignments on the parent sequence.

    While I'm sure it would be possible to write a script to manipulate all the alignments to all of the "split" subsequences back into the context of their original parent sequence, I am looking for a more elegant solution and trying to avoid this if possible.

    Is there a standard FASTA sequence description nomenclature for subsequences to accomplish what I'm describing, ideally that is recognized by programs like IGV?

    If there is not, I will likely still be splitting these sequences as described. Is there a suggested standard/common method for annotating/identifying these split sequence identifiers that maximizes readability?

    If relevant, the biological context for this question is analysis of RNA-Seq data collected after pulldown of a specific RNA-binding protein and digestion of overhanging read ends. The data should be enriched for a specific subsequence of interest, but will also contain the "parent" transcript. Alignments to both portions are of interest to me, but in different ways, and the alignments to the subsequence and the remainder of the parent sequence should be considered distinct. In terms of visualization it is helpful to view the alignments to all pieces in context of the "parent" transcript.

    I realize that splitting the parent sequence into two subsequences for alignment may result in loss of reads straddling the two portions, but that is acceptable given the necessity of having the two sections be distinct for alignment purposes.

    Thanks very much! -Dan

Latest Articles

Collapse

  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM
  • 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-28-2024, 06:12 AM
0 responses
27 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-23-2024, 04:11 PM
0 responses
74 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
81 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
69 views
0 likes
Last Post seqadmin  
Working...
X