• Intro

  • step_01
  • step_02
  • step_03
  • step_04
  • step_05
  • step_06
  • step_07
  • step_08
  • step_09
  • step_10
  • step_11
  • step_12
  • step_13
  • step_14
  • step_15
  • step_16
  • step_17
  • step_18
  • step_19
  • step_20
  • step_21
  • step_22
  • step_23
  • step_24

  • Credits



  • CGPDB
    UC Davis


    Step 12: SNPs/INDELs finding
    Genotype "A" versus "B" and "C"


    Finally, we are ready to run Python_CAP3_contig_poly_DIS_Feb_27_2004.py script to find all SNP/INDEL candidates in our assembly.

    Copy L_ABC_Contigs_All.aln.MM file into one directory up (where tomato_ABC.cap3.out.Info is located)

    $ cp L_ABC_Contigs_All.aln.MM ../

    $ cd ../


    We have run Python_CAP3_contig_poly_DIS_Feb_27_2004.py on two input files tomato_ABC.cap3.out.Info and L_ABC_Contigs_All.aln.MM. Example dialog is shown below:

    $ python Python_CAP3_contig_poly_DIS_Feb_27_2004.py
    Enter the mismatch file name : L_ABC_Contigs_All.aln.MM
    Enter the contig info file name : tomato_ABC.cap3.out.Info
    Distinguish by (p)refix or (s)uffix? : p
    Please input the size of the prefix : 1
    Specify Group 1
    Group 1 has []
    Please input element : A
    Do you want to add more? (y/n) : n
    Specify Group 2
    Group 2 has []
    Please input element : B
    Do you want to add more? (y/n) : y
    Group 2 has ['B']
    Please input element : C
    Do you want to add more? (y/n) : n

    ........ [ after you enter last "n" program will start SNP/INDEL search ]

    See description of output files on the next page