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!
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!
Comment