Hi all, I had done a illumina read mapping based on a particular gene sequence (reference sequence), which is highly repetitive in nature by using the following command
bowtie2-build -f ref_gene.fasta ref_index && bowtie2 -p 4 -N 1 -t -x ref_index -1 genome_R1.fastq -2 genome_R2.fastq -S tal.sam && samtools view -bS tal.sam | samtools sort - -o tal.bam && samtools mpileup -uf ref_gene.fasta tal.bam | bcftools call -c | vcfutils.pl vcf2fq > cns.fastq && seqtk seq -aQ64 -q20 -n N cns.fastq > contig.fasta
It is works perfectly fine and yielded exactly same length of mapped sequence as ref_gene. However, I need to get both end (reverse and forward end) extended overlapping reads in the assembled contig (Please see the image for better understanding) Is there any way to do it the same?
Thanks in advance.
1 answer
If I understood your question correctly you are looking at unaligned bases at the end of reads. That means you have to extract the soft clipped bases. Googling for that brings up this thread: extracting the soft clipped seq only from a sam file
I suspect you're looking at the wrong file.
Log in to answer this question.
-
Bcftools consensus generates mismatched consensus sequence
written by Duy •
-
Name of the consensus fasta sequence from a Minimap + Samtools alignment
written by lordbleys •
-
Topcoat error during resuming topcoat pipeline with unmapped reads
written by lingling •
-
TopHat error during aligning reads
written by lingling •
