SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   puzzle with gmapper -M mirna (shrimp) (http://seqanswers.com/forums/showthread.php?t=36042)

Sergio.pv 11-21-2013 05:55 AM

puzzle with gmapper -M mirna (shrimp)
 
Hi everyone,
i am mapping short reads from ion torrent PGM. Reads are QV~20 and 16-155 length.
In FastQC appears this as the highest represented read:
TGAGGTAGTAGGTTGTGTGGTT 19528counts

Runing a quick blast it matches 100% bta-let-7b (bta my sp of interest):
UserSe 1 ugagguaguagguugugugguu 22
bta-let-7b 1 ugagguaguagguugugugguu 22

Of course I can't do this with all my reads, so i mapped them to all ncRNAs and mirna-stem-loops for bos taurus, using shrimp in mirna mode:

gmapper-ls myfile.fa P1_adapter.fa bta-ncRNAs.fa -M mirna >myfile.out 2>myfile.log

After counting the repeated elements in the RNAME coulumn of the SAM file, bta-let-7b appears nowhere and neither the corresponding miRNAs for the other highly counted reads.

I have tried this again using same parameters with only stem loops and then only mature mirna sequences.

In all trials the resuls were the same.

Hope you can come up with soemthing, i will keep trying other settings for gmapper and I will let you know how it went.

Thaks!

Sergio.pv 11-21-2013 07:16 AM

I did the following, as a test:
first, i use the above metioned sequences (identical as you can see)
>bta-let-7b MIMAT0004331
UGAGGUAGUAGGUUGUGUGGUU
>dummy
UGAGGUAGUAGGUUGUGUGGUU

gmapper-ls dummy-seq.fa dummy-let-7b.fa -N 12 -o 10 -h 80% >map.out 2>map.log

Output:
HD VN:1.0 SO:unsorted
@SQ SN:bta-let-7b LN:22
@PG ID:gmapper VN:2.2.3 CL:gmapper-ls dummy-seq.fa dummy-let-7b.fa -N 12 -o 10 -h 80%
dummy 0 bta-let-7b 1 250 22M * 0 0 NGAGGNAGNAGGNNGNGNGGNN *

Which is ok.

But if I change the dummy sequence (U -> T), into the original output of the sequencer (Us as Ts) i don't get results:
@HD VN:1.0 SO:unsorted
@SQ SN:bta-let-7b LN:22
@PG ID:gmapper VN:2.2.3 CL:gmapper-ls dummy-seq.fa dummy-let-7b.fa -N 12 -o 10 -h 80%

It seems that the use of Ts is not computed. Is that correct?, is there an option to modifify that?

Best!

Sergio.pv 11-22-2013 02:45 AM

Ok, to make it work i converted U into T and it worked.
if it is any use, with fastx tool kit:

fasta_nucleotide_changer -v -i bta-stem-loops.mirbase20.fa -o bta-stemloops.fa -d


All times are GMT -8. The time now is 11:54 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.