Perl script

[bioenvisage]

Can any one help me with a perl script for removing the repeats in the reads , for eg i will paste the format of the seq below

HWI-EAS373:2:100:1792:1509#0/1 AAAAAAAAAAAAAAAAAACAACAAAAAAACAAAACAAAAACAAAACCAACACC ]_a`_Z_IT`b_\[_Ya\[\[]S\[RHUR^a^YY_V]aa^[TaW\Y\W_`^][aYR_BBBBBBBBBBBBBBBBBBB NF
HWI-EAS373:2:100:1792:1509#0/2 ACACACACATGGTCCACCATATTTTTTTACTTGGTTGTA aRaPZ__\__]VG[]RMGX\_Z_aa_P_NQ[_\VTFZTOa`R_[Q]ZZZXaBBBBBBBBBBBBBBBBBBBBBBBBB NF
HWI-EAS373:2:100:1792:1691#0/1 ACACACACAGTGTAGCTGGGGAGCAGGGATCCATTGATC abaa^]Waa]b_`Vb_b`aa[^`aa_aaXD^H]`]QWYa`ZaZaH]`TMS]`^BBBBBBBBBBBBBBBBBBBBBBB NF
HWI-EAS373:2:100:1792:1691#0/2 GGCTTTTTTGGTATCCTTTTCTCATGTTAGATGATGGGAGCATTTTTCTTCAGTgggatggatggtctggtagggc a^aY`_aaVa`UUabWaWa_bab_`a`b`aaOb``YN[a]GR`a`a`ba]_[J[XYBBBBBBBBBBBBBBBBBBBB NF
HWI-EAS373:2:100:1792:198#0/1 CGGCATTCCTTTTATTATAGCCCCTCTAGCTAGTTACAGTAGATAGGAACGtgcatgaatctntaaatggntgnan aZ]`]ab``aaab`a`]`YT`a^`aa`UZ\^X_Y]^Z^aYY[TYV[\XVLYBBBBBBBBBBBBBBBBBBBBBBBBB NF
HWI-EAS373:2:100:1792:198#0/2 agCTGATCTAGCGTCGTCTGCAACAACAACCGCGGGGGCGTCatcaacggcaagtgcggctcagcctcgggtgttg HOT_TTGYZGV_]GUQ_XNGSQZ\QIYTXT\_RKQGGL]O\ZBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NF
HWI-EAS373:2:100:1793:876#0/1 CCNCTGCCTCTACCTCCACGCCCTCGGCCTCTGCCACGCCCGCGGCCTGTATCTccagtgctctactcgcacanan `WDV^`a``a`aa^a`_aa[a]aY[`a\][a`\^``\`\]\^^S]Z[ZXW]ZP\SQBBBBBBBBBBBBBBBBBBBB NF

[maubp]

I'm having trouble making out what the data should look like - try wrapping the example with code tags,
[ code ] sequence data [ /code ]

[krobison]

By repeats in the reads do you mean
(a) Within a read, the repetition of a single nucleotide or simple sequence
(b) For sets of reads which are identical (presumed PCR duplicates), report only one
(c) Something else

SAMTools will do (b)

[Dave S.]

If you just want to remove homopolymers of DNA of some arbitrary length, use something like:

$min = 4;

while (<>)
{
s/(G){$min,}|(A){$min,}|(T){$min,}|(C){$min,}/$1$2$3$4/g;
print;
}

You are probably better off determining where and when they occur before wiping them out, e.g. see http://www.bioperl.org/wiki/Finding_...hes_in_contigs

If your sequencing method is generating spurious homopolymers you will need a much more sophisticated approach to determining which ones are real.

[bioenvisage]

hi krobison ...iam telling about with in the read the repetation of single nucleotide and also simple repeats.