Holiday Notice:
The Frontline Support team will be slow to respond December 17-18 due to an institute-wide retreat and offline December 22- January 1, while the institute is closed. Thank you for your patience during these next few weeks. Happy Holidays!

SplitNCigarReads unrecognized options

Hi,

I've been following the "Calling Variants in RNASeq" workflow, using the most recent version of gatk (v4.0.2.1).

When I run the SplitNCigarReads command I get the following error:

$gatk SplitNCigarReads -R /data/fasta/ensembl-Mus_musculus.fa -I /clscratch/kerrgr1/rnaseqVar/AB1-A1-medium/2Pass/Aligned.dedupped.bam -O /clscratch/kerrgr1/AAB312_rnaseqVar/AB1-A1-medium/2Pass/Aligned.split.bam -RF ReassignOneMappingQuality -RMQF 255 -RMQT 60 -U ALLOW_N_CIGAR_READS
Using GATK jar /home/kerrgr1/gatk-4.0.2.1/gatk-package-4.0.2.1-local.jar
Running:
java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=1 -jar /home/kerrgr1/gatk-4.0.2.1/gatk-package-4.0.2.1-local.jar SplitNCigarReads -R /db/nibrgenome/NG00012.0/fasta/ensembl-Mus_musculus.fa -I /clscratch/kerrgr1/AAB312_rnaseqVar/AB1-A1-medium/2Pass/Aligned.dedupped.bam -O /clscratch/kerrgr1/AAB312_rnaseqVar/AB1-A1-medium/2Pass/Aligned.split.bam -RF ReassignOneMappingQuality -RMQF 255 -RMQT 60 -U ALLOW_N_CIGAR_READS
USAGE: SplitNCigarReads [arguments]
....


A USER ERROR has occurred: U is not a recognized option


I also noticed in the usage options printed to screen that "ReassignOneMappingQuality" isn't an option for parameter "RF"

--read-filter,-RF:String Read filters to be applied before analysis This argument may be specified 0 or more
times. Default value: null. Possible Values: {AlignmentAgreesWithHeaderReadFilter,
AllowAllReadsReadFilter, AmbiguousBaseReadFilter, CigarContainsNoNOperator,
FirstOfPairReadFilter, FragmentLengthReadFilter, GoodCigarReadFilter,
HasReadGroupReadFilter, LibraryReadFilter, MappedReadFilter,
MappingQualityAvailableReadFilter, MappingQualityNotZeroReadFilter,
MappingQualityReadFilter, MatchingBasesAndQualsReadFilter, MateDifferentStrandReadFilter,
MateOnSameContigOrNoMappedMateReadFilter, MetricsReadFilter,
NonZeroFragmentLengthReadFilter, NonZeroReferenceLengthAlignmentReadFilter,
NotDuplicateReadFilter, NotOpticalDuplicateReadFilter, NotSecondaryAlignmentReadFilter,
NotSupplementaryAlignmentReadFilter, OverclippedReadFilter, PairedReadFilter,
PassesVendorQualityCheckReadFilter, PlatformReadFilter, PlatformUnitReadFilter,
PrimaryLineReadFilter, ProperlyPairedReadFilter, ReadGroupBlackListReadFilter,
ReadGroupReadFilter, ReadLengthEqualsCigarLengthReadFilter, ReadLengthReadFilter,
ReadNameReadFilter, ReadStrandFilter, SampleReadFilter, SecondOfPairReadFilter,
SeqIsStoredReadFilter, ValidAlignmentEndReadFilter, ValidAlignmentStartReadFilter,
WellformedReadFilter}

I also noticed that the online documentation specifies "-o" for the output file, whereas in practice it requires "-O".

Am I doing something blatantly wrong that I can't see?

Thanks for any help,
gkerrs

Answers

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @gkerrs
    Hi gkerrs,

    I think those first two options are not necessary in GATK4. Let me know if you get an error when running without them. You may also find the WDLs here helpful.

    -Sheila

    P.S. Thanks for pointing out the doc error. I will fix it now.

  • FPBarthelFPBarthel HoustonMember

    Hi @Sheila ! I am also trying to follow this workflow. Do I understand correctly that (to prepare input RNAseq BAM file for variant calling) in GATK 4 one no longer needs to supply the parameters "-RF ReassignOneMappingQuality -RMQF 255 -RMQT 60 -U ALLOW_N_CIGAR_READS" to gatk SplitNCigarReads? Thus, I only need to supply reference (-R), input BAM (-I) and output BAM (-O)?

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @FPBarthel
    Hi,

    It does seem that way from the documentation :smile: Please let us know if you have issues running with just the basic inputs.

    Thanks,
    Sheila

  • FPBarthelFPBarthel HoustonMember

    Thank you @Sheila ! It does seem to work fine and thanks for confirming. I wanted to be sure we are not missing something here.

Sign In or Register to comment.