Avatar

rcorbett

  • Total activity 14
  • Last activity
  • Member since
  • Following 0 users
  • Followed by 0 users
  • Votes 0
  • Subscriptions 3

Activity overview

Latest activity by rcorbett
  • Avatar

    rcorbett created a post,

    GATK-DRAGEN standalone pipeline?

    Answered

    Hi folks, I have been asked to try out your GATK-dragen pipeline.   I found this and related pages: https://gatk.broadinstitute.org/hc/en-us/articles/4411716682011-Full-release-of-open-source-DRAGE...

  • Avatar

    rcorbett commented,

    OK.  I think we're making some progress here.   I am starting with a refseq .gtf file where the annotation for my rRNA looks like this: chrUn_GL000220v1 BestRefSeq transcript 105424 118780 . + . ge...

  • Avatar

    rcorbett commented,

    Many thanks for your help Genevieve-Brandt-she-her. It is more that I am confused by the output.  I am not expecting a high UTR amount as my data are almost exclusively ribosomal and 45S, to my kno...

  • Avatar

    rcorbett created a post,

    CollectRnaSeqMetrics GL000220 test not counting ribosomal reads

    Answered

    If you are seeing an error, please provide(REQUIRED) : a) GATK version used: gatk-4.1.9.0 b) Exact command used: gatk-4.1.9.0/gatk CollectRnaSeqMetrics -I A56781.STAR.genome.sorted.bam.dup.bam.chrU...

  • Avatar

    rcorbett commented,

    OK. I uploaded a file called filterSamReads_output_format.tar.gz I was surprised to not be asked for a password at the following address, so it might be worth checking if the file is in fact on you...

  • Avatar

    rcorbett commented,

    Thank you. I tried the same command with 4.1.9.0 and got the same errors.  I will create a bug report.

  • Avatar

    rcorbett commented,

    Good afternoon. gatk ValidateSamFile -I passed.bam   -M SUMMARY Summary output: ## HISTOGRAM java.lang.StringError Type CountERROR:MISMATCH_FLAG_MATE_NEG_STRAND 14964150WARNING:ADJACENT_INDEL_IN_CI...

  • Avatar

    rcorbett commented,

    Hi Genevieve-Brandt-she-her, For 1. above:gatk PrintReads -I stdout.bam -O tmp.tmp --java-options '-DGATK_STACKTRACE_ON_USER_EXCEPTION=true'Using GATK jar /gsc/software/linux-x86_64-centos6/gatk-4....

  • Avatar

    rcorbett commented,

    Hi again, we tried the same command on a BWA 0.7.6a mem aligned bam and ran in to the same issue.

  • Avatar

    rcorbett commented,

    Thanks Genevieve, If I looks through the file I don't see anything human readable.  Similarly, if I try  strings stdout.bam I don't see anything of note. Just a little more information on the above...