Examples of SAMRecord

Examples of net.sf.samtools.SAMRecord

68
69
70
71
72
73
74
75
76
77
78
79
80
            Cigar newCigar = shiftCigarLeft(read.getCigar(), i);
           
            String shiftedReadAltRef = c2r.getAlternateReference(read, newCigar);
           
            if ((shiftedReadAltRef != null) && (origReadAltRef.equals(shiftedReadAltRef))) {         
              SAMRecord newRead = cloneRead(read);
              newRead.setCigar(newCigar);
              newRead.setAttribute("IS", read.getCigarString());
              return newRead;
            }
          }
        }
      }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

198
199
200
201
202
203
204
205
206
207
208
    c2r.init(ref);
   
    IndelShifter indelShifter = new IndelShifter();

    for (SAMRecord read : reader) {
      SAMRecord shiftedRead = indelShifter.shiftIndelsLeft(read, c2r);
      writer.addAlignment(shiftedRead);
    }
   
    writer.close();
    reader.close();
View Full Code Here

Examples of net.sf.samtools.SAMRecord

42
43
44
45
46
47
48
49
50
51
52
53
*/
   
//    String readStr = "chr11_118380033_118380233_2__numedges:341_totaledgecounts:133107_medianedgecount:417_minedgecount:4_terminatedatrepeat:false  0  chr11  118379934  239  235M3I164M  *  0  0  GGTTAGAATCAGAGAATATCAATGCTAAAAGGATTATGAGAATCACCCACTTTACCTACTTATTTTCTACATTTAAAAAAAAAAATCTAAGCTCCAAAGAAGTTAAGTGATTTGGCCCACATTGGACTGAAACTTGGCGCACCTGTCTCTCGGTGCAGTGTTCTTCCAGTACATATTGTGTGATCACCTGTCAGCTAAGGACTCAAGACCATACCCATACTCTTCTGCTGTACTGGCTGTTTTACCAGCACTGAGGCTTAAATAGCTAGTAATAACCTGACTTCACTTTTTAGTTGTTACTAAAGAAAACTAAGAACCATTTTTATTAGATAGTCAGATTTTGGTTACAATACCAGATACATCTCCATGGCATTTTCCATCAGTTCTAATGAATTTGATTAG  *  XE:i:9  XF:i:0  XN:i:0  AS:i:388";
    String readStr = "chr11_118380033_118380233_2__numedges:341_totaledgecounts:133107_medianedgecount:417_minedgecount:4_terminatedatrepeat:false  0  chr11  110000001  239  100M800D135M3I164M  *  0  0  GGTTAGAATCAGAGAATATCAATGCTAAAAGGATTATGAGAATCACCCACTTTACCTACTTATTTTCTACATTTAAAAAAAAAAATCTAAGCTCCAAAGAAGTTAAGTGATTTGGCCCACATTGGACTGAAACTTGGCGCACCTGTCTCTCGGTGCAGTGTTCTTCCAGTACATATTGTGTGATCACCTGTCAGCTAAGGACTCAAGACCATACCCATACTCTTCTGCTGTACTGGCTGTTTTACCAGCACTGAGGCTTAAATAGCTAGTAATAACCTGACTTCACTTTTTAGTTGTTACTAAAGAAAACTAAGAACCATTTTTATTAGATAGTCAGATTTTGGTTACAATACCAGATACATCTCCATGGCATTTTCCATCAGTTCTAATGAATTTGATTAG  *  XE:i:9  XF:i:0  XN:i:0  AS:i:388";
    SamStringReader rdr = new SamStringReader(new SAMFileHeader());
    SAMRecord read = rdr.getRead(readStr);
    System.out.println("read: " + read);
    System.out.println(read.getSAMString());
    System.out.println(read.getAlignmentStart());
    System.out.println(read.getAlignmentEnd());
    System.out.println(read.getReadLength());
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

29
30
31
32
33
34
35
36
37
    reads.add(newRecord(6162064, "129S120M", 249));
    reads.add(newRecord(6161985, "79M170S", 249));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 1);
    SAMRecord combined = outputList.get(0);
    assertEquals(combined.getAlignmentStart(), 6161985);
    assertEquals(combined.getCigarString(), "79M50I120M");
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

42
43
44
45
46
47
48
49
50
    reads.add(newRecord(6161985, "129M120S", 249));
    reads.add(newRecord(6162064, "129S120M", 249));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 1);
    SAMRecord combined = outputList.get(0);
    assertEquals(combined.getAlignmentStart(), 6161985);
    assertEquals(combined.getCigarString(), "79M50I120M");
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

55
56
57
58
59
60
61
62
63
64
    reads.add(newRecord(16085602, "60S65M", 125));
    reads.add(newRecord(16085517, "60M65S", 125));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 1);
    SAMRecord combined = outputList.get(0);
    assertEquals(combined.getAlignmentStart(), 16085517);
    assertEquals(combined.getCigarString(), "60M25D65M");
    assertEquals(combined.getReadLength(), 125);
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

69
70
71
72
73
74
75
76
77
78
    reads.add(newRecord(16085605, "49S62M", 111));
    reads.add(newRecord(16085517, "53M58S", 111));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 1);
    SAMRecord combined = outputList.get(0);
    assertEquals(combined.getAlignmentStart(), 16085517);
    assertEquals(combined.getCigarString(), "49M39D62M");
    assertEquals(combined.getReadLength(), 111);
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

84
85
86
87
88
89
90
91
92
    reads.add(newRecord(16085717, "300S200M", 500));
    reads.add(newRecord(33141553, "199S102M199S", 500));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 1);
    SAMRecord combined = outputList.get(0);
    assertEquals(combined.getAlignmentStart(), 16085517);
    assertEquals(combined.getCigarString(), "200M100I200M");
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

98
99
100
101
102
103
104
105
106
    reads.add(newRecord(30826191, "90S51M43S", 184));
    reads.add(newRecord(30726192, "141S43M", 184));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 1);
    SAMRecord combined = outputList.get(0);
    assertEquals(combined.getAlignmentStart(), 30726101);
    assertEquals(combined.getCigarString(), "90M50I44M");
  }
View Full Code Here

Examples of net.sf.samtools.SAMRecord

111
112
113
114
115
116
117
118
119
120
    reads.add(newRecord(11032656, "63M31S", 94));
    reads.add(newRecord(214990735, "34S60M", 94));
   
    List<SAMRecord> outputList = cc.processRead(reads);
    assertEquals(outputList.size(), 2);
    SAMRecord read1 = outputList.get(0);
    SAMRecord read2 = outputList.get(1);
    assertEquals(read1.getAlignmentStart(), 11032656);
    assertEquals(read2.getAlignmentStart(), 214990735);
  }
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.