I blasted a gene sequence against a whole genome and found out that the subject start site:31636628 is higher than the subject stop site:31626793. I need the coordinates to excise a portion of a bamfile but since the subject start site is greater than the stop site, samtools reports invalid region. Can someone help me as to how best to specify these coordinates in my bamfile?
↧