Thanks @snyquist2!! So removing the deltas and gammas from the reference fasta got me past that previous error, but now I'm stuck in call-recover_cdr3s:
================================================================================
18:58:58, 2021-06-21 - TCRGO: Initialized logger, ready for TCRGO.
18:58:58 - INFO: Parsing input data...
18:58:58 - INFO: No index was detected, building index...
18:59:01 - INFO: Reading /cromwell_root/fc-d8b64044-a999-46bf-846b-b52b2a9a1b7c/33d380fe-6f7e-4362-8472-660f91a7f896/TCRGO/ab90f7a0-15d5-4bf3-9329-dce024f1babe/call-alignment/out/pregnancy_test2_repairedsorted.bam...
18:59:01 - INFO: Reading CDR3 positions file and FASTA file...
18:59:01 - ERROR! AE000658|TRAV1-1*01| does not begin with 'TRA' or 'TRB'!
How is the reference fasta you're using syntactically constructed? Or better yet, could you share the human reference fasta you're using?
Originally posted by @ajwilk in #1 (comment)
Thanks @snyquist2!! So removing the deltas and gammas from the reference fasta got me past that previous error, but now I'm stuck in
call-recover_cdr3s:How is the reference fasta you're using syntactically constructed? Or better yet, could you share the human reference fasta you're using?
Originally posted by @ajwilk in #1 (comment)