Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • maq match command

    Dear all,

    I have a text file generated from the maq match command aligning paired end short reads to a reference genome. Any ideas how to filter out poor quality reads from this file, i.e. those reads that have been mapped to more than 1 location in the genome?

    How are people generally dealing with multiple hits from single read in the genome?

    Thanks for any help

    L

  • #2
    It all depends!

    Which text file are you referring to? Pileup? Dumped hits? If you specify how exactly you generated it, that will help others help you ...

    Also, for the map and downstream steps (consensus, SNP calling), maq puts a read that maps equally well (and satisfies the specified cutoffs) to multiple locations in the reference, in one of those locations ... randomly. This may or may not suit your needs, but as far as I know, there's no way to change it. You should be able to determine which reads map multiple times (and thus exclude them in a second round of mapping) by parsing the dumped hits file .. specified in the '-H' option to the map/match command ...

    Hope that helps
    ~Joe

    Comment


    • #3
      The file I was referring to was generated using the maq match command. I have re-run the maq match command using the -H and -u options. Do you think the multiple matches should be removed and the maq match command run again or would it be sufficient to remove the multiple matches and move on to do a pileup.

      At this stage I shall focus on correctly paired reads (flagged 18), remove multiple hits (flag of 0) and also low mapping quality scores (<30).

      Any other suggestions or comments how people would go about cleaning their chip-seq data?

      Cheers

      L

      Comment


      • #4
        Sounds like you're talking about the mapview output ... generated from the "mapview" command, using the binary map file generated by the map/match command. I'm not as familiar with that file - for instance, I didn't know that there was a flag for multiply mapped reads in mapview's output - but it sounds like you've got a good strategy for parsing that file and filtering your pairs.

        Comment


        • #5
          Thanx Joe..

          How did you convert the file created using the -H option in ./maq match command. The -H option was to generate the multiple hits and created a binary file. The ./maq mapview conversion does not work as it does for out.map. Is there a way to convert this binary file to text?

          Cheers
          L

          Comment


          • #6
            Layla, the file created with the -H option is not actually a binary file but a gzipped text file with information about the multiply mapped reads. I had the same confusion and finally figured this out!

            Comment


            • #7
              I figured the same fact as Owen explains, its gzipped!
              For multiply mapped reads from mapview result, the reads with 0 quality are mapped to multiple locations, using -q 1 should do the trick in excluding multiply-mapped reads

              Thoughts?
              --
              bioinfosm

              Comment


              • #8
                I think q -1 sounds like a valid option. Or you can simple grep for reads with the 0 flag and remove them before down-stream processing.

                Whilst on this note, anyone used SISSR instead of Maq? And if so, any thoughts on what to do with the data after SISSR gives 80,000 potential binding sites with p values < 0.001, high tag counts and fold changes?

                The data never simplifies!!!!

                L

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Strategies for Sequencing Challenging Samples
                  by seqadmin


                  Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                  03-22-2024, 06:39 AM
                • seqadmin
                  Techniques and Challenges in Conservation Genomics
                  by seqadmin



                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                  Avian Conservation
                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                  03-08-2024, 10:41 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Yesterday, 06:37 PM
                0 responses
                8 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                8 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                49 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                66 views
                0 likes
                Last Post seqadmin  
                Working...
                X