Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie Mapping Wrong Strand

    Hello,

    I ran bowtie 0.12.7 on a FASTA file of microarray probe sequences using the hg19 index as supplied on the website, and I noticed from manually inputting into BLAT a couple of sequences that sequences mapping to the negative strand map to + strand with BLAT. For the sequences mapping to + strand with bowtie, they map to the + strand with BLAT.

    Can someone check this ? How could it be

    Code:
    bowtie -f -v 0 -m 1 --un probesNotMapping.txt hg19 probeSeqs.txt probesMapped.txt
    e.g.

    Bowtie output :
    - chr16 54901845 AGTATTAAACAAGAATGCACGTAAAGTGTGCAATGCATATAGCAGACAGTCAAGAAATGGTGT

    BLAT output :
    100.0% 16 + 54901846 54901908 63

    For reproducibility, the entire sequences file is here.

  • #2
    bowtie: reads aligned to reverse strand are returned in reverse complement

    just now I thought I had run into the same problem. I was fooled by the fact that the output of bowtie contains the read sequences corresponding to the forward strand. So, when the flag in the SAM output of bowtie places the read on the reverse strand, the reported read sequence will now align on the forward strand.
    Check your fasta file you used for input and my guess is that the reads mapped to the reverse strand by bowtie are reverse complement of the reads in the bowtie output.
    The bowtie manual (http://bowtie-bio.sourceforge.net/bo...tml#sam-output) indeed specifies that reads in SAM output are reverse-complemented in aligned to the reverse strand.

    Ludo

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Genetic Variation in Immunogenetics and Antibody Diversity
      by seqadmin



      The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
      11-06-2024, 07:24 PM
    • seqadmin
      Choosing Between NGS and qPCR
      by seqadmin



      Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
      10-18-2024, 07:11 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Today, 11:09 AM
    0 responses
    24 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Today, 06:13 AM
    0 responses
    23 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 11-01-2024, 06:09 AM
    0 responses
    31 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-30-2024, 05:31 AM
    0 responses
    21 views
    0 likes
    Last Post seqadmin  
    Working...
    X