Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie preference for earlier alignments?

    So, I aligned my small RNA reads using the following command in bowtie.

    bowtie -S -a -v 1 -m 1 --best --strata Reference/Release4.3/dmel_rel4_3 2012/fastq2012/Adap_len_rRNA/A_final.fastq --max A_BDGP4_3_multireads.fastq A_BDGP4_3.sam

    The reference is for Drosophila chromosomes arm sequences. If I grep for the title I get this.

    >2L type=chromosome_arm; loc=2L:1..22407834; ID=2L; release=r4.3; species=Dmel;
    >4 type=chromosome_arm; loc=4:1..1281640; ID=4; release=r4.3; species=Dmel;
    >2R type=chromosome_arm; loc=2R:1..20766785; ID=2R; release=r4.3; species=Dmel;
    >3R type=chromosome_arm; loc=3R:1..27905053; ID=3R; release=r4.3; species=Dmel;
    >X type=chromosome_arm; loc=X:1..22224390; ID=X; release=r4.3; species=Dmel;
    >3L type=chromosome_arm; loc=3L:1..23771897; ID=3L; release=r4.3; species=Dmel;
    >dmel_mitochondrion_genome type=chromosome; loc=dmel_mitochondrion_genome:1..195


    I then used multiBamCov to find read counts for my regions of interest (~500) which are evenly distributed throughout the chromosomes. But, I find that ALL the alignments which match my region of interest are in 2L and 4. I know it could be that the coordinates are off for all the chromosomes except 2L and 4, but then I wonder whether there could be anything forcing a preference for alignments that happen earlier in the reference file (since 2L and 4 are the first chromosomes)? Any other problems I should check for?

Latest Articles

Collapse

  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM
  • seqadmin
    Current Approaches to Protein Sequencing
    by seqadmin


    Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
    04-04-2024, 04:25 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 11:49 AM
0 responses
13 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 08:47 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-11-2024, 12:08 PM
0 responses
61 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
60 views
0 likes
Last Post seqadmin  
Working...
X