I have performed SNPeff with RNA-seq results. Finally, I want to design primers but I can't know the sequences near SNP positions. I want to compare the sequences of the selected transcripts within samples or SNP between samples and the reference genome. I don't know what files and programs should be used. fastq.gz files? SAM files? BAM files? or VCF files? Is there an easy way to perform this?
1 answer
You could use a combo of bedtools getfasta to get the sequence you're interested in then use samtools view to get the mapped reads that fall into that region. I would then convert the output of samtools view into fastq files then do a multiple sequence alignment over your reference sequence. This method assumes you know where your SNPs are (approximately).
Log in to answer this question.
-
Geneious can't find all the sequences it needs in the selected file.
written by Timotheus •
-
Primer Design - how to choose primer for SNPs in Primer 3
written by Human •
-
Program/script to count how many sequences have mismatches on each position of my primers, and also…
written by gastonviarengo •
-
perl script for BWA-mem on multiple different files
written by brunobsouzaa
-
Aligning transcript with known position for SNP to draft genome to find the position of SNPs
written by dasakan •
-
RNA-seq samples along different time points, how to run Stringtie?
written by iraia.munoa
-
How to design primers when I only know the position of the mutation given in the VCF file?
written by jpxaviersantos90 •
-
Ideal practice to manage huge NGS dataset/Intermediate files of Whole Genome Human sequencing
written by ravi.uhdnis
-
Vcf-Compare Vcftools Error
written by Leandro Batista •
-
Qpcr: How And What Softwares You Used To Get All Transcripts Of One Gene?
written by Cheng Zhongshan
