Question: SNVer error: chrome not found in Reference
0
gravatar for anjali.gopal91
4.7 years ago by
United States
anjali.gopal9150 wrote:

Hi all,

I'm pretty new to bioinformatics, so please bear with me if I'm using incorrect terminology.

I'm trying to perform variant calling on a series of genes (NOT a whole genome). I have a fasta file that has reference sequences for all of these genes, and a series of fastq files for reads. I aligned these reads with Bowtie, got an output BAM file, and (since these are haploid genes) am now trying to use SNVer to perform variant calls.

The fasta file that I'm using for SNVer is the same fasta file that I used to make a bowtie index. However, whenever I run this, I'm getting the following error:

The chrome M=AD_2|ORF_ID=6222|GENE=10743|RhIY8ADG02_15:F02|ORF_SIZE=2901 is not found in Reference

Except I've checked the fasta file in question and the reference sequence (with the same title) is actually there.

Not sure what's going on. Is there some key information I'm missing? Is fasta not the appropriate format? Do I have to do pre-processing of the fasta file before I can pass it through SNVer? 

Any help would be much appreciated.

Thanks!
Anjali

ADD COMMENTlink modified 2.7 years ago by rse70 • written 4.7 years ago by anjali.gopal9150
0
gravatar for anjali.gopal91
4.7 years ago by
United States
anjali.gopal9150 wrote:

Never mind, I fixed this. Turns out my bam file was not sorted.

If anyone else has this error, you can get around it by:

samtools sort [filename].bam [output]

ADD COMMENTlink written 4.7 years ago by anjali.gopal9150
0
gravatar for rse
2.7 years ago by
rse70
Singapore
rse70 wrote:

SNVer generates a pretty big output file as compared to other variant callers. What can be the reason? Too many false positives?

ADD COMMENTlink written 2.7 years ago by rse70
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1686 users visited in the last hour