fixed bug where the wrong base was assumed for bam files

RickGelhausen · Jan 20, 2020 · 0eb7614 · 0eb7614
1 parent 36e89e9
commit 0eb7614
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/scripts/Ribo_seq_occupancy.py b/scripts/Ribo_seq_occupancy.py
@@ -73,7 +73,7 @@ def process_ribo(inputFile,occupancy,min_read_len,max_read_len,outputFileS,outpu
     for read in inFile.fetch():
         bit_flag = int(read.flag)   # FLAG;
         chromosome = read.reference_name  # chromosome
-        genomic_position = read.pos   # left-most position
+        genomic_position = read.pos+1   # left-most position
         cigar = read.cigarstring   # CIGAR string [only M and S observed in bacteria]
 
         length = len(read.query_sequence)     # footprint length