This is a long shot, but could the Target class include - strand features? My pacbio_amplicon.gb includes a single barcode and two genes in reverse orientations. My work-around is to run the pipeline twice with each gene on the + strand, then merge outputs on the barcodes.
file: alignparse/targets.py
lines: 162-166
if bio_feature.location.strand != 1:
raise ValueError(
f"feature {feature_name} of {self.name} is - "
"strand, but only + strand features handled"
)
This is a long shot, but could the Target class include - strand features? My
pacbio_amplicon.gbincludes a single barcode and two genes in reverse orientations. My work-around is to run the pipeline twice with each gene on the + strand, then merge outputs on the barcodes.file: alignparse/targets.py
lines: 162-166