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, 11:49 AM
|
0 responses
15 views
0 likes
|
Last Post
by seqadmin
Yesterday, 11:49 AM
|
||
Started by seqadmin, 04-24-2024, 08:47 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
04-24-2024, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
61 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
Comment