Dinh/Dinh 2012/NOTES/2012-2-13: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Dinh
(Created page with '==24 ChIP-seq - Stimulated and Resting Neurons== *Analysis at repeats. Downloaded RepeatMasker track fro mm9. Then filtered and get Non-overlapping windows: 2072 awk '{if($12…')
 
>Dinh
mNo edit summary
Line 12: Line 12:
  ./removeOverlapping < sorted.Repeats > tmp
  ./removeOverlapping < sorted.Repeats > tmp
  mv tmp sorted.Repeats (3,187,848 entries)
  mv tmp sorted.Repeats (3,187,848 entries)
[[File:removeOverlapping.txt]]
*Remove chr*_random entries
*Remove chr*_random entries
  awk '{if($1 !~ /_random/) print $0}' sorted.Repeats > sorted.Repeats_norandom
  awk '{if($1 !~ /_random/) print $0}' sorted.Repeats > sorted.Repeats_norandom

Revision as of 00:12, 15 February 2012

24 ChIP-seq - Stimulated and Resting Neurons

  • Analysis at repeats.
Downloaded RepeatMasker track fro mm9. Then filtered and get Non-overlapping windows:
 2072  awk '{if($12 ~ /DNA/) print $0}' RepeatMasker_mm9.txt | ./getNonOverlapping.pl > DNA_repeats.nonoverlapping &
2073  awk '{if($12 ~ /LINE/) print $0}' RepeatMasker_mm9.txt | ./getNonOverlapping.pl > LINE_repeats.nonoverlapping &
2074  awk '{if($12 ~ /SINE/) print $0}' RepeatMasker_mm9.txt | ./getNonOverlapping.pl > SINE_repeats.nonoverlapping &
2075  awk '{if($12 ~ /LTR/) print $0}' RepeatMasker_mm9.txt | ./getNonOverlapping.pl > LTR_repeats.nonoverlapping &
2076  awk '{if($12 ~ /Satellite/) print $0}' RepeatMasker_mm9.txt | ./getNonOverlapping.pl > Satellite_repeats.nonoverlapping &
File:GetNonOverlapping.txt
  • Note that some repeat segments overlaps between these difference repeat class.
cat *_repeats.nonoverlapping | sort -k1,1 -k2,3n > sorted.Repeats (3,188,316 entries)
./removeOverlapping < sorted.Repeats > tmp
mv tmp sorted.Repeats (3,187,848 entries)
File:RemoveOverlapping.txt
  • Remove chr*_random entries
awk '{if($1 !~ /_random/) print $0}' sorted.Repeats > sorted.Repeats_norandom
  • Final repeat regions file: 3,127,206 entries
Started with 3,188,316 non-overlapping entries for DNA, LINE, SINE, LTR, and Satellite repeats
Then removed entries/segments which have over-lapping labels -> 3,187,848 entries left
Then removed entries on chr*_random -> 3,127,206 entries -> Non-overlapping repeat regions