This is a test version of Biostars. For the public version, visit https://www.biostars.org.
identify the coordinate for coding and non_coding region.

Hi,

I would like to calculate the beginning and end positions for the coding and non coding regions in my genome sequence. is there any tool or script to do this ? my consensus sequence differ than the NCBI sequence. It has N stretch at the beginning.

Any help would be much appreciated. Thanks in advance

enter image description here

coding non_coding
• 927 views • link
written by G.S • G.S

Why do you have those N's at the beginning of the sequence? If the remainder of the sequence matches 100% then the initial N's may be wrong in your assembly.

• 0 views • link
written by GenoMax GenoMax

mmmm I am not sure. This is how I generate my consnsus sequence

 # Get consensus fastq file
samtools mpileup -uf  KT992094.1.fasta  seq-89_markup.bam | bcftools call -c | vcfutils.pl vcf2fq > seq-89_markup_sorted.fastq

# Convert .fastq to .fasta 
seqtk seq  seq-89_markup_sorted.fastq > seq-89.fasta
• 0 views • link
written by G.S • G.S

0 answers

No answers yet.

Log in to answer this question.

More posts like this