This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MAFTOOLS HELP - Mutational Signature

Hi All,

I am trying to run mutational signatures on my MAF file, using maftools (https://github.com/PoisonAlien/maftools). while trying to run this i am constantly My the following error. It would be great if you can help me. My knowledge in this area is elementary level.

I am using Homo_sapiens.GRCh38.dna_sm.primary_assembly.fa as my reference file. It would be of a GREAT HELP if you can help me sorting out this issue.

Please let me know the solution for this roadblock !!! Error as follows!!!

contigs in maf: [1] "chr10" "chr12" "chr22" "chr3" "chr1" "chr2" "chr4" "chr9" "chr7" [10] "chr16" "chr17" "chr19" "chr6" "chr13" "chr11" "chr20" "chr18" "chr5" [19] "chr14" "chr15" "chrX" "chr8" "chr21" "chrY" missing reference contigs from fasta file. [1] "chr10" "chr12" "chr22" "chr3" "chr1" "chr2" "chr4" "chr9" "chr7" [10] "chr16" "chr17" "chr19" "chr6" "chr13" "chr11" "chr20" "chr18" "chr5" [19] "chr14" "chr15" "chrX" "chr8" "chr21" "chrY"

Contig names in MAF must match to contig names in reference fasta. Ignorinig 29573 single nucleotide variants from chr10, chr12, chr22, chr3, chr1, chr2, chr4, chr9, chr7, chr16, chr17, chr19, chr6, chr13, chr11, chr20, chr18, chr5, chr14, chr15, chrX, chr8, chr21, chrY

Error in trinucleotideMatrix(maf = laml, ref_genome = "/home//mutation/centos/gr37.fasta", : Zero mutations to analyze! Maybe add or remove prefix?

Thanks a lot, Dav

mutational signature maftools
• 3,331 views • link
updated by shubhra.bhattacharya shubhra.bhattacharya • written by David_emir David_emir

1 answer

Hey!

Please check whether the chromosomal naming nomenclature matches of the reference matches with the input file. Such as that your reference has chromosome named as chr10,chr1 etc. I guess chromosomes are named as 1 ,2 and so on in the reference.

• 1 views • link
written by shubhra.bhattacharya shubhra.bhattacharya

Log in to answer this question.

More posts like this
  • How can I get this sorted by ascending order? (when the number and character is together)
    written by woojoy14 • woojoy14
    samtools view NA12878_dedup.bam | cut -f 3 | sort -u > chr.list8 > chr1 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr2 …
  • Removing a chromosome from a fasta file
    written by t86dan • t86dan
    Hello, I have a fasta file for human Genome GRCh37 Reference Assembly, but for some reason chromosome 20 is repeated and I would like to …
  • Error with correctGCBias in deepTools
    written by clmin • clmin
    Hello, I'm trying to correct the GC bias with deepTools. I successfully used computeGCBias > computeGCBias -b file.bam --effectiveGenomeSize 2150570000 -g mm9.2bit -l 200 --GCbiasFrequenciesFile …
  • diferent contigs chrM
    written by cristina_sabiers cristina_sabiers
    Hi! I try to call snps to learn how to do it. I got myt bam file, I hadn't made it. I have the snp, …
  • coverageBed with -g option ?
    written by bioguy24 bioguy24
    I am running a coverageBed and getting the below error. ``` cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -d -sorted -a /home/cmccabe/Desktop/NGS/bed/bedtools/xgen_targets_sorted.bed -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ERROR: Database file /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam …
  • problem with methylKit of the annotation
    written by hua.peng1314 hua.peng1314
    I am dealing with RRBS data follow [the manual][1]. I met some problem when doing annotation.The problem showed below. ``` > read.transcript.features("mm9.refseq.genes.bed",remove.unsual=TRUE,up.flank=1500,down.flank=500,unique.prom=TRUE) > diffAnn <- …
  • Splitting chromosomes in bins of 100kb
    written by ChIP ChIP
    Hi! I have the genomic table file of hg19, describing the lengths of each chromosomes. How can I split it up into bins of 100kb …
  • CONTRA problem. It does not stop.
    written by mafonso • mafonso
    Hi, I am trying to run an algorithm for CNV detection and was testing CONTRA. The test file works but using my own files the …
  • How To Retrieve Chromosome Lengths Using Bsgenome Package?
    written by James Ashmore • James Ashmore
    <p>Hello, I am following along with the <a href='http://www.ebi.ac.uk/training/online/sites/ebi.ac.uk.training.online/files/user/42/documents/embo_chipseq_course_1101026.pdf'>EBI's online course for ChIP-Seq analysis</a> and have encountered an error when trying to retrieve chromosome lengths:</p> …
  • Error Import File With "Import" Of Rtracklayer
    written by ska007 • ska007
    <p>Hi, This is related to a the question found here (http://biostar.stackexchange.com/questions/3414?sort=newest#sort-top). I am trying to import a gff file using import("file.gff") which contains a set …