This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to get percent identity between two RNA-seq fastq files?

Hello, Does anybody have an idea about the way to get percent identity between two or more RNA-seq fastq files ? or to show how many of the sequences in one of them are found in another one?

thanks in advance

percent identity rna-seq
• 1,633 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by lkianmehr lkianmehr

1 answer 

 comm \
     <(gunzip  -c 11.fq.gz | paste - - - - | cut -f 2 | sort)  \
     <(gunzip  -c f2.fq.gz | paste - - - - | cut -f 2 | sort) | \
     awk -F '\t' '{if($1=="" && $2=="") {C+=1.0;};} END {print C/NR;}'
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I have run that command on two fastq files that are sequenced on the same library, I got one value = 0.119817, what does it mean? it means they have just 0.1 percent difference?

• 0 views • link
written by lkianmehr lkianmehr

Log in to answer this question.

More posts like this