RN Sarma
- Total activity 18
- Last activity
- Member since
- Following 0 users
- Followed by 0 users
- Votes 0
- Subscriptions 7
Activity overview
Latest activity by RN Sarma-
RN Sarma created a post,
No output from GATK
gatk HaplotypeCaller \ -R "$REFERENCE_GENOME" \ -I "$BAM_FILE" \ -O "$GVCF_FILE" \ -ERC GVCF \ -G StandardAnnotation \ --tmp-dir "$TMP_DIR"
-
RN Sarma created a post,
GATK Command for variant discovery
Is the script correct for GTK version 4.5 for SNP discovery? #!/bin/bash#@(#)GATK_4.50-slurm_bioconda3_rns.sh 2025-01-06 R.N.Sarma # Define pathsREFERENCE_GENOME="/home/rnsarma/Desktop/Tea_Reseq...
-
RN Sarma commented,
The usage of gatk GatherVcfs is:gatk GatherVcfsCloud \ -I cohortA_chr1.vcf.gz \ -I cohortA_chr2.vcf.gz \ -O cohortA_chr1chr2.vcf.gz 1. Why is there is GAtherVcfCloud, as my output files...
-
RN Sarma commented,
My script, as shown below, has produced folders like "genomicsdb_1" to "genomicsdb_11" in the directory "/home/rnsarma/Desktop/Mungbean/Combined_GVCF". How to get one VCF file for all chromosomal...
-
RN Sarma commented,
Java version: openjdk version "11.0.23" 2024-04-16OpenJDK Runtime Environment (build 11.0.23+9-post-Ubuntu-1ubuntu122.04.1)OpenJDK 64-Bit Server VM (build 11.0.23+9-post-Ubuntu-1ubuntu122.04.1, mix...
-
RN Sarma commented,
The script ended abruptly with this line at the end " A USER ERROR has occurred: genomicsdb-workspace is not a recognized option ********************************************************************...
-
RN Sarma created a post,
GenomiCDBImport issue
REQUIRED for all errors and issues: a) GATK version used:GATK4This script is asking me some more parameter. Kindly check the script and suggest. #!/bin/bash# "combined_GVCF_genomicDB_normal_rns_...
-
RN Sarma created a post,
Consolidating GVCF
GATK vesion 4.8 Which one should I use to to combine 107 GVCF of a diploid genome of 2n=450 Mb size, GenomicsDBImport or CombineGVCFs?
-
RN Sarma created a post,
Idetify biallelic loci in GATK 4,8
Is my codes correct to identify biallelic loci of diploid genome? #!/bin/bash gatk HaplotypeCaller \-R /home/rnsarma/Desktop/Mungbean/mungbean_genome/Vigna_radiata.Vradiata_ver6.dna.toplevel.fa \-...
-
RN Sarma created a post,
Checking of GATK script
Is this script will retain only biallelic loci and remove indels in gataGATKversion 4.5.0? Organism is diploid, purpose is to identify bialelic SNP, not indels. gatk --version being used is:Using ...