Error in PGAP pangenome pipeline?
Hi
Have anybody used PGAP pangenome pipeline, if it so please help me to fix an issue. I have used both 2012 PGAP version and recent PGAP-X version aswell and both shows error for user input. But works well for the given test data.
sequence
genome
sequencing
gene
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
The GATK "The given bam input has no sample names." error
written by Begonia_pavonina 22for f in MINIMAP BWA ; do ~/gatk-4.2.0.0/gatk HaplotypeCaller --reference /home/tmichel/projects/rbge/HybSeq_thibauld/reference_genomes/Begonia_loranthoides_scaffold.fasta --input Hillebrandia_sorted.$f.bam --output Hillebrandia.$f.g.vcf.gz --emit-ref-confidence GVCF ; done I have used GATK to call …
-
how to define open-pangenome?
written by Kumar 12I have done `pan-genome` analysis using `roary` pipeline and I determined that my bacterial genome datasets posses `open-pangenome` based on roary pipeline generated `pan-genome progress …
-
error in Islandviewer4 API.
written by Kumar 12I have been predicting genomic island using [IslandViewer4][1] batch upload option (HTTP API), I have followed the following command as they mentioned but doing the …
-
error while executing panX script (KeyError: 'db_xref').
written by Kumar 12I have been working with [panX][1]. I ran the test dataset without any error. But, while running my genome of interest ([AF391541.1][2] [AY585228.1][3] [AY700211.1][4] ), …
-
Variant calling using GATK::ERROR
written by Kumar 17Dear All, I am working to setup pipelines for variants calling. I used SOAPnuke, HISAT, for pre-process. Finally, I am using GATK to call snps/variants. …
-
samtools mpileup: bam file to consensus sequence generation error?
written by Kumar 12I have been doing reference based assembly of illumina based paired end reads. For which, I indexed my reference genome and generated sam, bam and …
-
devtools Installation error in R (version 3.5.3)?
written by Kumar 12I am trying to install devtools in R (R version=3.5.3, OS= Ubuntu 18.04.2 LTS) by using following commands, ```r install.packages("devtools") install.packages("devtools") devtools::install_github("r-lib/devtools") ``` Its getting …
-
SSR`s identification and primer designing
written by moohit21 0i am using GMATA software for SSR`s identification and primer designing, i install the programme as per direction given in the manual but while running …
-
Badly formed genome location? GATK GenotypeGVCFs
written by jtwalker 2I'm trying to genotype >1000 samples using the GATK pipeline. I've already created gVCFs for my samples, but when attempting to use GATK's GenotypeGVCFs tool, …
-
PseudoPipe is very hassle to configure
written by Kumar 17Hi, I tried [PseudoPipe][1]. It is very hassle to configure, even with its example data itself. there is no any README and Manual file to …
So it is likely that there is an error in your input. Have you compared your inputs with the test data to see if you are able to spot what may be wrong?
Dear @genomax, Initially I used the test data and it works fine as I mentioned earlier. Then I have used my data including the assembled draft genome (.fna) of bacterial strains and their corresponding .ptt files as specified in the manual genome files (.fna) kept in GenomeSequence folder and annotation files (.ptt) kept in GenomeAnnotation folder. Both folder have been kept in input folder. I tried to run the executive file
pgaxand started my analysis by clickingGet Startbutton and then I have clickedGenome Alignmentbutton, the process have been completed successfully. Further, I followedOrthologs analysisby clicking the same, and chosen Nucleotide Sequence Based with 70 % Identity and Coverage as well. But this time it ran without giving result. The resulting folder namelyGenomeOrthologsClusterfolder having 1kb sizedOrthologsClusterfile, which has the following strings onlyCluster ID citri oryzae phas. Even I tried with the complete genome strains, but the result is same. I ran this inwindows 10 OS.Please give us a reproducible example that shows evidence of the pipeline working with test data but failing on different data, including the test data, your input data, the exact error message as well as your thoughts on why the error could be happening.
Dear @RamRS, I checked my input data, which is same like test data. It differs in terms of size only, my each genome file has 5MB, whereas the test data is around 2.5 to 5 MB. I tried both complete and draft genome sequences as well, but both stops working in
Orthologs analysisprocess. If it possible I can upload my dataset for your reference. Please let me know, how to upload my dataset here.