SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Adapter trimming figo1019 RNA Sequencing 2 07-17-2018 04:00 AM
About adapter trimming starchildcn Illumina/Solexa 1 04-28-2015 07:28 AM
Best adapter trimming? ronton Bioinformatics 12 08-21-2014 03:04 PM
Adapter trimming and trimming by quality question alisrpp Bioinformatics 5 04-08-2013 04:55 PM
3' Adapter Trimming caddymob Bioinformatics 0 05-27-2009 12:53 PM

Reply
 
Thread Tools
Old 10-04-2019, 06:44 AM   #1
Anastaziel
Junior Member
 
Location: Moon

Join Date: Jun 2015
Posts: 7
Default miRDeep2 and adapter trimming

I am trying to align a miRNA-seq FASTQ file using the mapper.pl from miRDeep2 but I am not having any luck.

The reads in the FASTQ file are ~50bp.

I am not sure where my mistake is as I get 46380 reads mapped (less than 1.6% of total reads). This is my code

Code:
mapper.pl /mnt/d/work/NNNN_1_SAMPLE1_ATCACGA_S2_L001_R1_001.fastq -e -h -j -k ATCACGA -l 16 -m -p GRCh38_primary_assembly_genome -s /mnt/d/work/tools/mirdeep2/test/2.collapsed_reads.fa -
t /mnt/d/work/tools/mirdeep2/test/2.grch38.arf -n -v
I think my mistake is with the adapter sequence. The adapter for this sample is ATCACGA (it is in the sample name too). I looked at this website illumina-adapter-and-primer-sequences/ and my samples (I do have more) follow the TruSeq™ DNA Sample Prep Kit v2 but if I try to put a longer adapter (TruSeq Adapter, Index 1) I just get +- 10000 mapped reads.

I am not sure how to proceed. Any help would be appreciated.

Last edited by Anastaziel; 10-04-2019 at 11:11 AM.
Anastaziel is offline   Reply With Quote
Old 10-04-2019, 11:10 AM   #2
Anastaziel
Junior Member
 
Location: Moon

Join Date: Jun 2015
Posts: 7
Default

I think I fixed it. The adapter for my files should be TGGAATTCTCGGGTGCCAAGG. From what-sequences-do-i-use-for-adapter-trimming under the category of TruSeq Small RNA.

Now I get ~60% mapping. I am not sure if this is sufficient, though.
Anastaziel 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 02:08 AM.


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