Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Counting reads in miRNA experiment sara_ Epigenetics 3 01-26-2016 04:48 PM
miRNA seq analysis - large numbers of non-aligning reads lbeltrame Bioinformatics 10 05-03-2013 09:30 AM
miRNA aligning/counting kmcarr Bioinformatics 28 02-08-2013 02:35 PM
aligning miRNA-seq reads to other short sequence data impeachme2 Bioinformatics 2 07-03-2012 08:34 AM
GATK to discover Single Nucleotide Variation in mature miRNA from miRNA-Seq Bioinfo83 Bioinformatics 0 01-31-2012 05:11 AM

Thread Tools
Old 03-26-2018, 09:09 AM   #1
Junior Member
Location: Blacksburg VA

Join Date: Mar 2018
Posts: 1
Default Question: counting after aligning miRNA reads to miRBase mature miRNA instead of refe

I have human miRNA-seq data. In the past I have aligned these data using Bowtie2 to the reference genome (Homo_sapiens.GRCh38.dna.toplevel.fa), and I have subsequently performed counting using featureCounts with the annotation file hsa.gff3 from miRBase. . Now I have aligned the reads to the mature miRNA from miRBase (mature.fa), but when I look at the resulting bam files, the reads have a flag of 4 (segment unmapped) and there is no position information in mature.fa.

So can I use mature.fa as the reference for alignment and if so how is counting performed?

Thanks, Ina is offline   Reply With Quote
Old 03-26-2018, 01:39 PM   #2
Location: Santa Cruz, CA

Join Date: Feb 2017
Posts: 17

Hi Ina,

I use samtools and picard-tools to get the counts for miRNAs.

Here is what I would do for bowtie2 alignment to miRBase
bowtie2 -L8 --local -x miRBase-mature-hsa-index -U FILE.fastq -S FILE.sam
samtools view -b -S FILE.sam > FILE.bam
samtools sort FILE.bam > FILE.sorted
picard-tools BuildBamIndex I=FILE.sorted
picard-tools BamIndexStats I=FILE.sorted > FILE.txt

Makes sense?

sbarberan is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:41 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO