Filtering sam or bam file with maximum matching region
1
0
Entering edit mode
28 days ago
analyst ▴ 50

Hi all,

I am analysing NGS amplicon reads of crisper gene samples. I do not know exactly which genes do they belong to because if I align reads against gene sequence (that was supposed to be) using BWA it gives 0% identity however by aligning reference genome I get 99% alignment rate. Now how can I find about the gene through reads data. Is there any way to find through sam or bam file that which chromosomal position are maximum reads aligning to?

Thanks

crisper-edited alignment match • 603 views
ADD COMMENT
3
Entering edit mode
28 days ago

. Is there any way to find through sam or bam file that which chromosomal position are maximum reads aligning to?

samtools depth in.bam | awk -F '\t' 'BEGIN{BEST_DEPTH=0;BEST="";} {DEPTH=int($3);if(BEST_DEPTH<DEPTH) {BEST_DEPTH=DEPTH;BEST=$0;} } END{print BEST;}'
ADD COMMENT
0
Entering edit mode

Thanks Pierre!

I got this output

Can you please interpret the output

enter image description here

ADD REPLY
2
Entering edit mode

Please do not paste screenshots of plain text content, it is counterproductive. You can copy paste the content directly here (using the code formatting option shown below),

ADD REPLY
0
Entering edit mode

Okay sure Pierre I will take care of it next time

ADD REPLY
2
Entering edit mode

the max depth is 3520 at chr1:80660513

ADD REPLY
0
Entering edit mode

Thanks for the interpretation Pierre!

ADD REPLY

Login before adding your answer.

Traffic: 1418 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6