Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • whsqwghlm
    Member
    • Jun 2009
    • 14

    #16
    Smallest NM value? Sorry - you lost me...

    The idea is to record the hit(s) with the longest identical 5' match(es) to the genome, the theory being that primer artefacts, sequencing errors and RNA edits are all concentrated at the 3' end. We also assume that natural variation is absent for miRNAs. If we get multiple matches with the same score, then all of the matches are recorded.

    Comment

    • staylor
      Member
      • Feb 2009
      • 17

      #17
      Clearly not using NM then!:-) In SAM format, NM = the number of nucleotide differences to the reference sequence. I thought this may be a useful tag for filtering. Or do you just count the length of the match?

      Do you reject sequences longer than 22bp?

      Comment

      • didymos
        Junior Member
        • Jun 2010
        • 8

        #18
        I am new in the miRNA field and I am wondering why you are using -k 200 or 101 option? In other words why you want to have 200 alignment with 0 mismatches, rather than one unique with 0 mm?
        Thanks!

        tomek

        Comment

        • whsqwghlm
          Member
          • Jun 2009
          • 14

          #19
          Each read may map exactly to many places in the genome. We want to capture all of these locations to a threshold promiscuity, typically 100, over which we discard all of the mappings (i.e. if 101 alignments are returned from the search).

          Comment

          • yjhua2110
            Member
            • Nov 2009
            • 68

            #20
            Originally posted by whsqwghlm View Post
            Each read may map exactly to many places in the genome. We want to capture all of these locations to a threshold promiscuity, typically 100, over which we discard all of the mappings (i.e. if 101 alignments are returned from the search).
            you can use the options: -a -m 100. Specifying -m 100 instructs bowtie to refrain from reporting any alignments for reads having more than 100 reportable alignments.

            Comment

            • dukevn
              Member
              • Apr 2009
              • 50

              #21
              Anybody knows how I can trim out the adapter sequence with a certain number of mismatch (1 or 2)?

              Thanks,

              D.

              Comment

              • didymos
                Junior Member
                • Jun 2010
                • 8

                #22
                Hi,

                In Bowtie you can trim n nucleotides from 3' or 5' end with -3 n or -5 n option - so you can play with this and -v 2 for two mm. However for trimming I am using blast with smaller word size -W option, and then choosing only those mappings which are on the 3' end.

                tomek

                Comment

                • dukevn
                  Member
                  • Apr 2009
                  • 50

                  #23
                  Originally posted by didymos View Post
                  Hi,

                  In Bowtie you can trim n nucleotides from 3' or 5' end with -3 n or -5 n option - so you can play with this and -v 2 for two mm. However for trimming I am using blast with smaller word size -W option, and then choosing only those mappings which are on the 3' end.

                  tomek
                  Thanks, but I want to trim the adapter sequence, not n nucleotides.

                  Comment

                  • mgogol
                    Senior Member
                    • Mar 2008
                    • 197

                    #24
                    fastx_clipper can trim adapter sequence.

                    Comment

                    • jdanderson
                      Member
                      • Sep 2010
                      • 45

                      #25
                      Hello Everyone,

                      I would like to second mgogol's suggestion. It is part of a larger set of tools called FASTX Toolkit: http://hannonlab.cshl.edu/fastx_toolkit/index.html

                      I've found the tools to be very useful and user friendly.

                      Cheers,
                      Johnathon

                      Comment

                      • zeam
                        Member
                        • Oct 2010
                        • 43

                        #26
                        Hi everyone,
                        When you were processing smRNA-seq data,which software did you use?I have attempted SOAP2 and bowtie,both of them didn't work well enough,many of the reads could not map to the reference.But when I used blat,it works well except consuming time.Can you guys give me some suggestions?My data is from plants,and I want to sepate miRNA and siRNA from smRNA-seq data.
                        Thanks!

                        Comment

                        • NicoBxl
                          not just another member
                          • Aug 2010
                          • 264

                          #27
                          did you allow mismatches in your alignment ?

                          Comment

                          • zeam
                            Member
                            • Oct 2010
                            • 43

                            #28
                            No mismatch allowed

                            Originally posted by NicoBxl View Post
                            did you allow mismatches in your alignment ?
                            What I want is to dissect the correlation between smRNA and DNA methylation.Previous studies reveals that both of miRNA and siRNA involed in DNA methylation.The problem I met is 1)mapping;2)separating miRNA and siRNA from smRNA-seq data.Can you give me some suggestions if you have some experience on these subjects?

                            Comment

                            • unique379
                              Member
                              • Aug 2012
                              • 27

                              #29
                              Dear all,

                              as above discussions and suggestion, could i use bowtie something like this for miRNA alignment ?? just please make me sure.

                              bowtie --sam -q -n 0 -l 15 -e 99999 -k 200 -a --chunkmbs 1024 --best ....

                              Comment

                              • JonB
                                Member
                                • Jan 2010
                                • 84

                                #30
                                Originally posted by yjhua2110 View Post
                                in our deepBase database, we use options: –k 200 –v 0. the Specifying the parameters (–k 200 –v 0) instructs Bowtie to report up to 200 perfect hits for each read.

                                deepBase is a platform for annotating and discovering small and long ncRNAs from next generation sequencing data. It is available at http://deepbase.sysu.edu.cn
                                I am mapping reads from a small RNA sequencing, and I am curious how you treat the results after mapping with the -k 200 option? Do you assemble the transcripts? Or do you report all the locations where a read has been mapped?
                                Last edited by JonB; 06-10-2013, 05:26 AM. Reason: lacked quote

                                Comment

                                Latest Articles

                                Collapse

                                • seqadmin
                                  Pathogen Surveillance with Advanced Genomic Tools
                                  by seqadmin




                                  The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
                                  Yesterday, 11:48 AM
                                • seqadmin
                                  New Genomics Tools and Methods Shared at AGBT 2025
                                  by seqadmin


                                  This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                                  The Headliner
                                  The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                                  03-03-2025, 01:39 PM
                                • seqadmin
                                  Investigating the Gut Microbiome Through Diet and Spatial Biology
                                  by seqadmin




                                  The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
                                  02-24-2025, 06:31 AM

                                ad_right_rmr

                                Collapse

                                News

                                Collapse

                                Topics Statistics Last Post
                                Started by seqadmin, 03-20-2025, 05:03 AM
                                0 responses
                                26 views
                                0 reactions
                                Last Post seqadmin  
                                Started by seqadmin, 03-19-2025, 07:27 AM
                                0 responses
                                33 views
                                0 reactions
                                Last Post seqadmin  
                                Started by seqadmin, 03-18-2025, 12:50 PM
                                0 responses
                                25 views
                                0 reactions
                                Last Post seqadmin  
                                Started by seqadmin, 03-03-2025, 01:15 PM
                                0 responses
                                190 views
                                0 reactions
                                Last Post seqadmin  
                                Working...