Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • how to convert BLAT output (pslx or maf) to gff, sam, wig to load in browser

    I was able to view BLAT output (psl format) in the IGB (igb.bioviz.org) but would like to be able to use the BROAD IGV instead which doesn't recognize any of the BLAT outputs.

    Is there a way to convert a psl, pslx, or maf format to sam, wig, gff or any other format that will allow me to view my sequences using IGV?

    Thank you!

  • #2
    I'm pretty new to the field of bioinformatics, but I have been dealing with a very similar situation, taking the PSL output of BLAT and putting it into a GFF file.

    I've been processing the PSL files on Galaxy (http://main.g2.bx.psu.edu/). It pretty easy and only takes a minute to convert it over into a GFF.

    Steps:
    1.) Upload the PSL file and select the file type as tabular. (get data>upload file)
    2.) Remove the header, if there is one (text manipulation>remove beginning). There are 5 lines in the header.
    3.) Add columns for the source and feature fields. (text manipulation>add column)
    4.) Add columns for the score, frame, and group (text manipulation>add column)
    For these I just put "."
    5.) Now that you have all the fields necessary for the GFF it is just a matter of cutting the columns in order. (text manipulation>cut).
    6.) Click on the pencil to edit the name of the file and to change it into GFF format and the columns will be named with the proper GFF fields.

    The first time it took me a while to figure out, but now I can go from PSL to GFF in about 1-2 minutes using Galaxy. There is probably an easy way to do it with programming and scripts, but this way has been working well for me.

    Hope this helps.


    -Brandon

    Comment


    • #3
      Thanks for your reply. I haven't tried it but I will. I also found convenient to get to a sam output which could be sorted and indexed with the IGV tools and then displayed in the browser.

      Comment


      • #4
        Could you please describe the method you used to convert Blat to SAM?

        Comment


        • #5
          blat (psl) to sam

          Samtools (v0.1.7) comes with a tool to convert psa to sam format called psl2sam.pl.

          Comment


          • #6
            We are trying to use Cufflinks to analyze our data. This requires CIGAR matching and it appears that psl2sam.py is unable to do this. Are there any other scripts or flags we can use for this particular application?

            Comment


            • #7
              Hi David - I would also like to use cufflinks to analyze BLAT aligned sequences. Did you find a solution to converting psl to sam with CIGAR matching?
              Thanks!
              Lindy McBride - Rockefeller University

              Comment


              • #8
                how to convert BLAT output (pslx or maf) to gff, sam, wig to load in browser

                Originally posted by rdeborja View Post
                Samtools (v0.1.7) comes with a tool to convert psa to sam format called psl2sam.pl.
                Dear all,
                I've been trying to use this perl script (psl2sam.pl) to convert a .psl file (from BLAT).
                However, this is not generating an output file (.sam)... instead, it just prints (what I think is) the output on shell

                $path/psl2sam.pl myfile.psl

                I'm not using any of the options in usage, but they are all related to the conversion itself and not about the output file...
                I must add that I never used SAMtools before...
                what am I missing here?

                thank you in advance

                Comment


                • #9
                  Originally posted by pbarros View Post
                  Dear all,
                  I've been trying to use this perl script (psl2sam.pl) to convert a .psl file (from BLAT).
                  However, this is not generating an output file (.sam)... instead, it just prints (what I think is) the output on shell

                  $path/psl2sam.pl myfile.psl

                  I'm not using any of the options in usage, but they are all related to the conversion itself and not about the output file...
                  I must add that I never used SAMtools before...
                  what am I missing here?

                  thank you in advance
                  ok!
                  just found the answer for that!
                  never mind

                  Comment


                  • #10
                    I'm pretty late to this thread, missed the original post. However, for the record IGV displays blat output in the form of .psl files, the extension must be .psl. There's no need to convert it to another format.

                    Jim

                    Comment


                    • #11
                      Hi, got to the same problem.
                      If you don't need mapping quality, I would get maf
                      (`blat -out=maf`) and then convert maf to sam (`maf-convert.py sam`) from last aligner last/src.

                      Comment


                      • #12
                        Why don't you just load the blat output (.psl file)?

                        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, 03-27-2024, 06:37 PM
                        0 responses
                        12 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-27-2024, 06:07 PM
                        0 responses
                        11 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-22-2024, 10:03 AM
                        0 responses
                        53 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-21-2024, 07:32 AM
                        0 responses
                        69 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X