Skip to content

Instantly share code, notes, and snippets.

@inodb
Created May 16, 2014 15:14
Show Gist options
  • Save inodb/80557dadf991fd4a038b to your computer and use it in GitHub Desktop.
Save inodb/80557dadf991fd4a038b to your computer and use it in GitHub Desktop.
Map contigs against itself with MUMmer. Then find duplicate contigs
for row in $(cat nucmer.coords | awk -v OFS=, '{if ($10 == 100.00 && $11 == 100.00 && $7 == 100.00 && $12 != $13) {print $12,$13}}')
do
echo $row | tr ',' '\n' | sort | tr '\n' ',' | sed 's/,$//'
echo
done | cut -d, -f2 | sort -u > duplicatecontigs.txt
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment