Lep-Anchor Code

Linkage map guided genome anchoring

Brought to you by: lep-map

[78dfa4]: / longsubread.awk Maximize Restore History

Download this file

23 lines (22 with data), 318 Bytes

#pick the longest subread from pacbio data
{
	if (/^>/) {
		split($1, a, "/")
		nn = a[1]"/"a[2]
		if (nn != name && longest != "") {
			print longest
			maxLen = 0
		}
		name = nn
		fname = $1
	} else {
		len = length($0)
		if (len >= maxLen) {
			longest = fname"\n"$0
			maxLen = len
		}
	}
}
END{
	print longest
}