I'm calling these attributes of the feature class: feature.location, feature.type, feature.strand

from the record.features in SeqIO.parse(gbk_file,"genbank")

Output:

[0:6183918](+) source 1
join{[1871:2015](-), [1776:1872](-), [1751:1775](-), [1419:1752](-), [1304:1418](-), [1051:1303](-), [810:1050](-), [722:809](-), [634:721](-), [336:633](-)} CDS -1
[2126:2723](-) CDS -1
[2674:2959](-) CDS -1
[2906:3869](-) CDS -1

I know the (-1) goes with the (-) strand and (1) with the (+) strand.

Is the slice object referencing the (+) strand even when the feature is on the (-) strand? In other words, the 3rd entry, is [2126:2723] on the + strand but the reverse-complement is the feature OR is it in relation to the opposite end of the strand?

Also, what is going on with this join object?