miRNA seq trimming
2
0
Entering edit mode
7.1 years ago
bioLife ▴ 50

I am analyzing miRNA sequencing data derived from a typical Illumina protocol for small RNAs. The problem is that although I have trimmed the reads for the adaptors, there are a lot of reads with length longer than 26bp.

What are these long reads? Should I trim the reads for hairpin sequence too? What's been actually sequenced? Else, what are the steps needed before aligning the reads to the reference sequence?

many thanks

miRNA NGS trimming • 3.3k views
ADD COMMENT
0
Entering edit mode

which species? I have found a lot of tRNA, that could or couldn't be functional. If you are interested to know what they are, just map to ncrnadb, and see what are the longer reads. But if not, just follow Devon advice. Just make sure it happens the same everywhere.

ADD REPLY
3
Entering edit mode
7.1 years ago

There's more than just miRNAs in a typical smallRNAseq experiment. You'll also have piRNAs, snRNA, snoRNAs, and so on. If all that you care about are the miRNAs, then just focus on them and ignore the rest.

ADD COMMENT
2
Entering edit mode
7.1 years ago
Manvendra Singh ★ 2.2k
######## First step is to trim your file quite carefully, I do it with Cutadapt
/usr/local/bin/cutadapt --discard-untrimmed --minimum-length="Number" --maximum-length="Number" -a <adapter_sequence> In_seq.fastq > your_trimmed_file.fastq

###### then remove also the reads mapping to ribosomal and tRNA sequences
bowtie --seedlen=23 --un output_file.fastq /path_to/bowtieindex/r_tRNA your_trimmed_file.fastq > /dev/null

your output_file.fastq should look better to align.

HTH

ADD COMMENT

Login before adding your answer.

Traffic: 2334 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6