SEQanswers

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 03:48 PM
miRNA seq analysis - large numbers of non-aligning reads lbeltrame Bioinformatics 10 05-03-2013 08:30 AM
miRNA aligning/counting kmcarr Bioinformatics 28 02-08-2013 01:35 PM
aligning miRNA-seq reads to other short sequence data impeachme2 Bioinformatics 2 07-03-2012 07:34 AM
GATK to discover Single Nucleotide Variation in mature miRNA from miRNA-Seq Bioinfo83 Bioinformatics 0 01-31-2012 04:11 AM

Reply
 
Thread Tools
Old 03-26-2018, 08:09 AM   #1
[email protected]
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
inah@bi.vt.edu is offline   Reply With Quote
Old 03-26-2018, 12:39 PM   #2
sbarberan
Member
 
Location: Santa Cruz, CA

Join Date: Feb 2017
Posts: 16
Default

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?

Cheers,
sergio
sbarberan is offline   Reply With Quote
Reply

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:09 PM.


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