Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie and samtools - FASTA formatting and genomic coordinates

    I'm in the process of aligning some RNA-Seq reads against a transcriptome. My transcriptome FASTA has headers in the following format:

    >13405 chr4:34234-34507

    As you might have guessed, the first part is a transcript ID. It seems that bowtie2 and samtools all use the transcript ID as the chromosome name and print out coordinate locations correspondingly. Is there any way to format the fasta files such that bowtie and/or samtools will use genomic coordinates instead? I could write a small program that converts the coordinates, but I'd like to know if there is an easier way.

  • #2
    This is a complicated problem if there is any splicing involved (introns to skip), or worse trans-splicing.

    Since you apparently have chromosome sequences, why not just map to them directly as normal? That should make the downstream analysis much easier.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Addressing Off-Target Effects in CRISPR Technologies
      by seqadmin






      The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...
      08-27-2024, 04:44 AM
    • seqadmin
      Selecting and Optimizing mRNA Library Preparations
      by seqadmin



      Sequencing mRNA provides a snapshot of cellular activity, allowing researchers to study the dynamics of cellular processes, compare gene expression across different tissue types, and gain insights into the mechanisms of complex diseases. “mRNA’s central role in the dogma of molecular biology makes it a logical and relevant focus for transcriptomic studies,” stated Sebastian Aguilar Pierlé, Ph.D., Application Development Lead at Inorevia. “One of the major hurdles for...
      08-07-2024, 12:11 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 08-27-2024, 04:40 AM
    0 responses
    16 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 08-22-2024, 05:00 AM
    0 responses
    293 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 08-21-2024, 10:49 AM
    0 responses
    135 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 08-19-2024, 05:12 AM
    0 responses
    124 views
    0 likes
    Last Post seqadmin  
    Working...
    X