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
                          Genetic Variation in Immunogenetics and Antibody Diversity
                          by seqadmin



                          The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
                          Today, 07:24 PM
                        • seqadmin
                          Choosing Between NGS and qPCR
                          by seqadmin



                          Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                          10-18-2024, 07:11 AM

                        ad_right_rmr

                        Collapse

                        News

                        Collapse

                        Topics Statistics Last Post
                        Started by seqadmin, 11-01-2024, 06:09 AM
                        0 responses
                        24 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 10-30-2024, 05:31 AM
                        0 responses
                        21 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 10-24-2024, 06:58 AM
                        0 responses
                        25 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 10-23-2024, 08:43 AM
                        0 responses
                        56 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X