Hi,
How do i count total and unique reads (may be with 1 mm) on each chromosome from bam file ?
Thanks,
How do i count total and unique reads (may be with 1 mm) on each chromosome from bam file ?
Thanks,
You are currently viewing the SEQanswers forums as a guest, which limits your access. Click here to register now, and join the discussion
awk '{if (substr($1,1,1)!="@") print $3}' file.sam | sort | uniq -c
awk '{if (substr($1,1,1)!="@") print $3,$10}' DNT_Input.sam | sort | uniq | awk '{print $1}' | uniq -c
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:37 PM
|
0 responses
8 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:37 PM
|
||
Started by seqadmin, Yesterday, 06:07 PM
|
0 responses
8 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:07 PM
|
||
Started by seqadmin, 03-22-2024, 10:03 AM
|
0 responses
49 views
0 likes
|
Last Post
by seqadmin
03-22-2024, 10:03 AM
|
||
Started by seqadmin, 03-21-2024, 07:32 AM
|
0 responses
67 views
0 likes
|
Last Post
by seqadmin
03-21-2024, 07:32 AM
|
Comment