# HG changeset patch # User jackcurragh # Date 1654618836 0 # Node ID 1be4bbe849e3819ee47cbfd2e31726a6a7b3e796 # Parent 8b059b0a19dc16f939111ef3b51b2c0e76fb3a0e Uploaded diff -r 8b059b0a19dc -r 1be4bbe849e3 trips_bam_to_sqlite/bam_to_sqlite.py --- a/trips_bam_to_sqlite/bam_to_sqlite.py Fri May 27 10:14:00 2022 +0000 +++ b/trips_bam_to_sqlite/bam_to_sqlite.py Tue Jun 07 16:20:36 2022 +0000 @@ -288,9 +288,8 @@ master_offset_dict = {"fiveprime": {}, "threeprime": {}} master_metagene_stop_dict = {"fiveprime": {}, "threeprime": {}} - os.system(f'samtools sort -n {bam_filepath} -o {bam_filepath}_n_sorted.bam') pysam.set_verbosity(0) - infile = pysam.Samfile(f"{bam_filepath}_n_sorted.bam", "rb") + infile = pysam.Samfile(f"{bam_filepath}", "rb") header = infile.header["HD"] unsorted = False if "SO" in header: diff -r 8b059b0a19dc -r 1be4bbe849e3 trips_bam_to_sqlite/trips_bam_to_sqlite.xml --- a/trips_bam_to_sqlite/trips_bam_to_sqlite.xml Fri May 27 10:14:00 2022 +0000 +++ b/trips_bam_to_sqlite/trips_bam_to_sqlite.xml Tue Jun 07 16:20:36 2022 +0000 @@ -1,10 +1,9 @@ - + Convert BAM file to SQLITE for Trips-Viz pysam sqlitedict sqlite - samtools