Hi Jason,
since you are adding a very large number of junction (~7M), you would need to do it at the genome generation, specifying
--limitSjdbInsertNsj 10000000 . This should solve the GstrandBit problem.
Such a large number of junctions may lead to a significant slowdown and increase in multi-mapping rate, so I generally recommend to filter the junctions to bring their number to <1M.
Cheers
Alex