Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • jamalshah888
    replied
    Here's an example of how you can sort a SAM file in Python using the pysam library:

    code:
    import pysam # Open the input SAM file samfile = pysam.AlignmentFile("input.sam", "r") # Create an output file for the sorted SAM data sorted_samfile = pysam.AlignmentFile("sorted.sam", "wb", template=samfile) # Sort the input SAM file by read name for read in samfile.fetch(until_eof=True): sorted_samfile.write(read) # Close the input and output files samfile.close() sorted_samfile.close()
    This code will open an input SAM file called "input.sam", sort the reads by read name, and write the sorted data to an output file called "sorted.sam".


    You can also sort the SAM file by other criteria such as reference name, position, etc. by using the appropriate sort function from pysam. For example, to sort the SAM file by reference name and position, you can use the following code:


    code:
    sorted_samfile = pysam.AlignmentFile("sorted.sam", "wb", template=samfile) sorted_samfile.sort("-o", "sorted.sam", "input.sam") sorted_samfile.close()
    This code will use the pysam sort function to sort the input SAM file by reference name and position and write the sorted data to the output file "sorted.sam".

    Regards
    Jamal Shah

    Leave a comment:


  • elviszboncak
    started a topic Help with sorting sam files in python!

    Help with sorting sam files in python!

    Hello All!

    I am new to coding. I have written a python script to analyze a bunch of fastq files from RNASeq analysis. I am using the subprocess module to call bowtie2 in order to index the genome and align the unpaired short reads to the genome.My script generates a sam file which needs to be sorted before running the cuffdiff module. The linux sort command works perfectly on my mac, however, in spite of installing the Unxupdates on my windows at work, I am unable to make my script work. I am hoping for some help/advice here to sort my sam file in python. Here is the script I tried.

    #!/usr/bin/env python2.7

    import os, sys

    inputFile = open(‘Outputfile1.txt’, 'r')

    lineList = inputFile.readlines()

    for line in sorted(lineList):

    print(line.rstrip());

    With this script, the stdout is written to the screen. But this recognizes only the text file. Sam files are not recognized.

    Any help is much appreciated.

    Thanks in advance

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