This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Non IUPAC error on hg38 GATK FASTA

hi all, i downloaded the Homo_sapiens_assembly38.fasta from GATK cloud including all the indices.

when i try to run GATK TargetRealignerCreator I am getting error, non IUPAC code found in reference: 70

what does this mean? I have not manipulated the reference FASTA in any way.

also i have tried dos2unix Homo_sapiens_assembly38.fasta and it not make any difference.

any ideas anyone?

ngs
• 3,457 views • link
written by win win

output of

$ file Homo_sapiens_assembly38.fasta
$ grep -v '>' Homo_sapiens_assembly38.fasta | tr -d "\n\r ATGCNatgcn"

?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

The first commands output is ASCII Text The second commands output is FMRY@RFYWYYKWMRYWBYYWB@FYRYSYRYYURRRYWMWWYMKYYWYKKMYRRRWRRYYBRRRRKWRS@YYMOYKYRRSKYWYYRYRYKMMRRWYRYYYRSW

what does this mean and how may it be removed?

• 0 views • link
written by win win

it's not a fasta file. I would say it's just a FASTQ file.

what is the output of 

$ head Homo_sapiens_assembly38.fasta
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum 
    >chr1  AC:CM000663.2  gi:568336023  LN:248956422  rl:Chromosome  M5:6aef897c3d6ff0c78aff06ac189178dd  AS:GRCh38
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
• 0 views • link
written by win win

@win -- did you find a solution to your issue?

• 0 views • link
written by Cricket • Cricket

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Help with error in GATK variant calling
    written by Chris • Chris
    Hi all, I try some reference genome such as Homo_sapiens_assembly38.fasta and Homo_sapiens.GRCh38.dna_sm.primary_assembly.fa but I still got the error below. Would you please have a suggestion? …
  • GATK3 error message
    written by shuksi1984 • shuksi1984
    I ran following command: java -Xmx4g -jar /path/to/gatk3.8.jar -T RealignerTargetCreator -R /path/to/Homo_sapiens_assembly38.fasta -I SP1_markduplicate.bam -o SP1_realignedtargets.list I got following error message: ERROR MESSAGE: Bad input: …
  • Reference allele is too long message GATK
    written by win win
    hi all, i generated a gVCF file using HaplotypeCaller. When i run ValidateVariants from GATK i get the message "Reference allele is too long (108) …
  • QD Annotation error in GATK in VariantRecalibration step
    written by win win
    hi all, i downloaded a BAM file from 1000 genomes to work on and converted it to FASTQ and aligned with bwa mem to hg38, …
  • CosmicCodingMuts.vcf.gz for hg38?
    written by Tania Tania
    Hi all I downloaded `CosmicCodingMuts.vcf.gz` from Cosmic site. Then I use for my RNAseq somatic calling `Homo_sapiens_assembly38.dbsnp138.vcf` from GATK resource bundle. And also use reference: …
  • Human reference files in HG38 GATK resource bundle
    written by win win
    hi all, maybe i missed something obvious in the GATK resource bundle (cloud bucket) at the following site. `https://console.cloud.google.com/storage/browser/genomics-public-data/resources/broad/hg38/v0?pli=1` they have this > Homo_sapiens_assembly38.fasta file …
  • BWA-MEM alignment error with hg38 from GATK
    written by win win
    hi all, i downloaded a bam file from 1000 genomes HG00100.chrom20.ILLUMINA.bwa.GBR.exome.20121211.bam then i converted this to fastq using sudo java -jar algorithms/picard/picard.jar SamToFastq VALIDATION_STRINGENCY=SILENT I=data/HG100.bam …
  • RealignerTargetCreator (GATK) contig issue
    written by win win
    So i downloaded a .bam from phase 3 1000 genomes and i want to run the RealignerTargetCreator on it using GATK 3.8. The following command …
  • GATK "We encountered a non-standard non-IUPAC base in the provided reference: '13"
    written by Phismil • Phismil
    Hello Team, I am trying to run GATK java -jar GenomeAnalysisTK.jar -T RealignerTargetCreator -R reference.fasta -I header.bam -o realigned.intervals My reference is a trinity assembled …
  • Reference File From Bam
    written by win win
    Hi all, Hope someone can help. I will be getting several BAM files from users and they may not inform us which reference file to …