New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3704467 =29 (0.950)7 (0.260) 3/116 2.6 20.2% coding (295/903 nt) dppC dipeptide/heme ABC transporter permease
?NC_000913 3704731 = 30 (1.120)coding (31/903 nt) dppC dipeptide/heme ABC transporter permease

ACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGG  >  NC_000913/3704664‑3704795
                                                                   |                                                                
acCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCt                                                                  >  2:133602/1‑68 (MQ=255)
      ccTTTGTCGCGTCTAAAATAGTGCCAGAACCCCTGTAACGGGGGCATCGGCACCGGCGCGCTAATCAc                                                            <  2:264658/68‑1 (MQ=255)
        tttGTTGCGTTTAAAATAGCGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACtt                                                          <  2:1107402/68‑1 (MQ=255)
           gttgCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTAt                                                       <  2:1024425/68‑1 (MQ=255)
               cGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTc                                                   <  1:42736/68‑1 (MQ=255)
               cGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTc                                                   <  1:404403/68‑1 (MQ=255)
                    aaaaTTTTGCCATATCTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCTCTTTATTTTCAGTat                                              <  2:466466/68‑2 (MQ=255)
                     gatgacgtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTAATCACTTTATTTTCAGTAAc                                             >  2:206175‑M2/47‑68 (MQ=255)
                      aaTAGTGCGATAACTCCTGTTACGGGGTCATCGTCACTGGTGCTATAATCACTTTATTTTCAGTTAcc                                            <  2:129910/68‑1 (MQ=255)
                      aaTAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAAcc                                            <  1:871076/68‑1 (MQ=255)
                       aTAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCt                                           >  1:796388/1‑68 (MQ=255)
                         aGTGCCAGTTCTCCTGTTTCGGGGTCATCGGCTCCGGTTCGCTTTTCACTTTATTTTCAGTAACCtgt                                         <  2:365399/68‑1 (MQ=255)
                         aGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCtgt                                         <  1:588298/68‑1 (MQ=255)
                         aGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCtgt                                         <  2:1126637/68‑1 (MQ=255)
                          cgtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTAATCACTTTATTTTCAGTAACCtgtg                                        <  1:206175‑M2/27‑1 (MQ=255)
                           tGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGa                                       <  1:1090584/68‑1 (MQ=255)
                           tGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGa                                       <  2:885309/68‑1 (MQ=255)
                           tGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGa                                       <  1:175694/68‑1 (MQ=255)
                           gtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTTATCACTTTATTTTCAGTAACctgtgt                                       <  1:405368‑M2/28‑2 (MQ=255)
                           gtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTAATCACTTTATTTTCAGTAACCTGTGa                                       <  2:709964‑M2/28‑1 (MQ=255)
                           gtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTAATCACTTTATTTTCAGTAACCTGTGa                                       <  1:178855‑M2/28‑1 (MQ=255)
                           gtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTAATCACTTTATTTTCAGTAACCTGTGa                                       <  2:40556‑M2/28‑1 (MQ=255)
                           gtaggccgtgatgtgctgtcgcgcctgatgtacggtgcgcTAATCACTTTATTTTCAGTAACCTGCGa                                       <  2:689723‑M2/28‑1 (MQ=255)
                           cgCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGa                                       <  2:1096689/67‑1 (MQ=255)
                              cAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACat                                    <  1:481963/68‑1 (MQ=255)
                                    tCCTGTAACGGGGTCAGCGGCACCTGTGCGCTAATCACTTTATTTTCAGTAACCTGCGACATGATGTc                              <  2:913699/68‑1 (MQ=255)
                                    tCCTGTAAATGGGTCATAGGAACCGGTGCGCTAATCACTTTATTATCAGTAAACTGTGACATGAAGTc                              >  1:250939/1‑68 (MQ=255)
                                     ccTGTAACGGGGCCATCGGCACCGGCGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTcc                             <  2:159558/68‑1 (MQ=255)
                                          aaCGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTa                        >  1:972122/1‑68 (MQ=255)
                                          aaCGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTa                        >  1:218922/1‑68 (MQ=255)
                                            cGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTAct                      >  1:876459/1‑68 (MQ=255)
                                            cGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTAct                      >  2:83105/1‑68 (MQ=255)
                                               ggTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAGCCTGTGACATGATGTCCCCTTActtct                   >  2:297644/1‑68 (MQ=255)
                                               ggTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTActtct                   >  1:843084/1‑68 (MQ=255)
                                               ggTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTActtct                   >  1:651496/1‑68 (MQ=255)
                                                     cGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGAc             <  1:149914/68‑1 (MQ=255)
                                                       gCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGa           <  2:426831/68‑1 (MQ=255)
                                                              tgcgcTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATAcgcg    >  1:520340‑M2/6‑68 (MQ=255)
                                                               gcgcTAATCACTTTATTTTCAGTATCCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCgg   <  2:837874‑M2/64‑1 (MQ=255)
                                                                cgcTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCggg  <  2:220340‑M2/65‑1 (MQ=255)
                                                                   |                                                                
ACCGCGCCTTTGTTGCGTTTAAAATAGTGCCAGAACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGG  >  NC_000913/3704664‑3704795

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.