I have illumina short reads which I want to call SNP without using a reference genome. I want to make a consenses sequence of some selected reads and then align the other reads to it. so the basic idea is that the consensus sequence need to be act like a "dummy genome" Is there any way to do this?
•
written
by
shynilasanthi •
2 answers
I think SNP is at the base position that the base in one genome (ex: genome A) is different from the base in the other genome (ex: genome B).
From your question, it seems that the reads derived from genome A and genome B are mixed in the data, and
genome A = sequences produced by some selected reads;
genome B = sequences produced by other reads.
I wonder how you extract reads of genome A from the original data, but if you can do that, I agree with toralmanvar.
You should assemble the reads derived from genome A, and you will get contigs that you called “dummy genome”.
Then you align the reads derived from genome B to the contigs of genome A, and you will detect SNPs.
I recommend that you also align the reads derived from genome A to the contigs of genome A at that time, to compare the allele frequency of two genomes.
Log in to answer this question.
-
Can BWA align genome fasta to fasta?
written by kgwkk2 •
-
When to use Genome assembly and when to use short read alignment for consensus genration??
written by Anisur Rahman •
-
how to align many single end fastq files from the same study on each other without using any refere…
written by FadyNabil •
-
Best way to assemble a genome with a reference
written by genseq •
-
No-clip alignment of scaffold against reference
written by adhb •
-
SNP calling from a scaffold.fa file
written by shynilasanthi •
-
constructing a supercontig from selected reads
written by shynilasanthi •
-
Assembly using two reference genomes
written by deepti1rao •
-
Making a custom reference genome
written by biogirl
-
How To Call Snp And Indel Based On A Sanger Sequences Based Bac Contig Vs Reference?
written by William
