New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1255973 =72 (0.700)12 (0.140) 4/40 NT 19.0% coding (130/2868 nt) ycgV predicted adhesin
?REL606 = 3223859 42 (0.490)coding (537/696 nt) yraI predicted periplasmic pilin chaperone

ACCCAACTGAATATTACGTCGTATTTGCCGGATTGGCAGTCGATAAAATCGATCTCACAAGCAAAATT  >  REL606/3223825‑3223892
                                  |                                 
gcccAACTGAATATTACGTCGTATTTGCCGGATTg                                   <  1:4839643/34‑1 (MQ=255)
aCCCAACTGAATATTACGTCGTATTTGCCGGATTg                                   >  1:7327935/1‑35 (MQ=255)
aCCCAACTGAATATTACGTCGTATTTGCCGGATTg                                   <  1:14230321/35‑1 (MQ=255)
 cccAACTGAATATTACGTCGTATTTGCCGGATTGg                                  <  1:17789208‑M2/35‑2 (MQ=255)
 cccAACTGAATATTACGTCGTATTTGCCGGATTGg                                  <  1:18700468‑M2/35‑2 (MQ=255)
 cccAACTGAATATTACGTCGTATTTGCCGGATTGg                                  >  1:1324206‑M2/1‑34 (MQ=255)
 cccAACTGAATATTACGTCGTATTTGCCGGATTGg                                  >  1:13410612‑M2/1‑34 (MQ=255)
  ccAACTGAATATTACGTCGTATTTGCCGGATTGgc                                 <  1:16487033‑M2/35‑3 (MQ=255)
  acaACTGAATATTACGTCGTATTTGCCGGATTGgc                                 >  1:1428666‑M2/2‑33 (MQ=255)
    aaCTGAATATTACGTCGTATTTGCCGGATTGgcag                               >  1:189368‑M2/1‑31 (MQ=255)
    aaCTGAATATTACGTCGTATTTGCCGGATTGgcag                               >  1:14751169‑M2/1‑31 (MQ=255)
     aCTGAATATTACGTCGTATTTGCCGGATTGGCAGt                              <  1:7008532/35‑1 (MQ=255)
     aCTGAATATTACGTCGTATTTGCCGGATTGGCAGt                              <  1:11031648/35‑1 (MQ=255)
      cTGAATATTACGTCGTATTTGCCGGATTGgcaggg                             <  1:4100215‑M2/35‑7 (MQ=255)
      cTGAATATTACGTCGTATTTGCCGGATTGgcaggg                             <  1:9596495‑M2/35‑7 (MQ=255)
       tGAATATTACGTCGTATTTGCCGGGTTGGCAGTCg                            <  1:13640949/35‑1 (MQ=255)
       tGAATATTACGTCGTATTTGCCGGATTGgcagggg                            <  1:20539765‑M2/35‑8 (MQ=255)
       tGAATATTACGTCGTATTTGCCGGATTGGCAGTCg                            <  1:14133670/35‑1 (MQ=255)
       tGAATATTACGTCGTATTTGCCGGATTGGCAGTCg                            <  1:16501593/35‑1 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGgcaggggg                           <  1:17261712‑M2/35‑9 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGgcaggggg                           <  1:7830448‑M2/35‑9 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGgcaggggg                           <  1:4602869‑M2/35‑9 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGgcaggggg                           <  1:11375488‑M2/35‑9 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGgcaggggg                           <  1:27625‑M2/35‑9 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGgcaggggg                           <  1:3485255‑M2/35‑9 (MQ=255)
        gAATATTACGTCGTATTTGCCGGATTGGCAGTCGa                           <  1:7725751/35‑1 (MQ=255)
         aaTATTACGTCGTATTTGCCGGATTGgcagggggg                          <  1:6567349‑M2/35‑10 (MQ=255)
         aaTATTACGTCGTATTTGCCGGATTGGCAGTCGAt                          <  1:16140491/35‑1 (MQ=255)
         aaTATTACGTCGTATTTGCCGGATTGGCAGTCGAt                          <  1:10375571/35‑1 (MQ=255)
          atatTACGTCGTATTTGCCGGATTGgcaggggggg                         <  1:15299865‑M2/35‑11 (MQ=255)
          atatTACGTCGTATTTGCCGGATTGgcaggggggg                         <  1:3945430‑M2/35‑11 (MQ=255)
          atatTACGTCGTATTTGCCGGATTGgcaggggggg                         <  1:10400319‑M2/35‑11 (MQ=255)
          atatTACGTCGTATTTGCCGGATTGGCAGTCGATa                         <  1:16501002/35‑1 (MQ=255)
          atatTACGTCGTATTTGCCGGATTGGCAGTCGATa                         <  1:7856784/35‑1 (MQ=255)
           tatTACGTCGTATTTGCCGGATTGgcaggggggga                        <  1:3219985‑M2/35‑12 (MQ=255)
           tatTACGTCGTATTTGCCGGATTGGCAGTCGATaa                        <  1:17822723/35‑1 (MQ=255)
           tatTACGTCGTATTTGCCGGATTGGCAGTCGATaa                        <  1:18925069/35‑1 (MQ=255)
            atTACGTCGTATTTGCCGGATTGgcagggggggaa                       <  1:9875090‑M2/35‑13 (MQ=255)
            atTACGTCGTATTTGCCGGATTGGCAGTCGATaaa                       <  1:10477216/35‑1 (MQ=255)
             ttACGTCGTATTTGCCGGATTGgcagggggggaaa                      <  1:6505031‑M2/35‑14 (MQ=255)
             ttACGTCGTATTTGCCGGATTGGCAGTCGATaaaa                      <  1:2580426/35‑1 (MQ=255)
             ttACGTCGTATTTGCCGGATTGGCAGTCGATaaaa                      <  1:8979145/35‑1 (MQ=255)
             ttACGTAGTATTTGCCGGATTGGCAGTCGATaaaa                      <  1:13940513/35‑1 (MQ=255)
              tACGTCGTATTTGCCGGATTGGCAGTCGATAAAAt                     <  1:2587131/35‑1 (MQ=255)
              tACGTCGTATTTGCCGGATTGGCAGTCGATAAAAt                     <  1:18785772/35‑1 (MQ=255)
              tACGTCGTATTTGCCGGATTGGCAGTCGATAAAAt                     <  1:16952846/35‑1 (MQ=255)
               aCGTCGTATTTGCCGGATTGgcagggggggaaaaa                    <  1:19415412‑M2/35‑16 (MQ=255)
               aCGTCGTATTTGCCGGATTGgcagggggggaaaaa                    <  1:5174791‑M2/35‑16 (MQ=255)
               aCGTCGTATTTGCCGGATTGgcagggggggaaaaa                    <  1:2173898‑M2/35‑16 (MQ=255)
               aCGTCGTATTTGCCGGATTGgcagggggggaaaaa                    <  1:11405117‑M2/35‑16 (MQ=255)
               aCGTCGTATTTGCCGGATTGgcagggggggaaaaa                    <  1:16270749‑M2/35‑16 (MQ=255)
               aCGTCGTATTTGCCGGATTGgcagggggggaaaaa                    <  1:16706786‑M2/35‑16 (MQ=255)
               aCGTCGTATTTGCCGGATTGGCAGTCGATAAAATc                    <  1:14115574/35‑1 (MQ=255)
               aCGTCGTATTTGCCGGATTGGCAGTCGATAAAATc                    <  1:17658761/35‑1 (MQ=255)
                cgtcgtATTTGCCGGATTGGCAGTCGATAAAATCg                   <  1:14550250/35‑1 (MQ=255)
                cgtcgtATTTGCCGGATTGGCAGTCGATAAAATCg                   <  1:14577969/35‑1 (MQ=255)
                cgtcgtATTTGCCGGATTGGCAGTCGATAAAAACg                   <  1:7619334/35‑1 (MQ=255)
                 gtcgtATTTGCCGGATTGGCAGTCGATAAAATCGa                  <  1:15987406/35‑1 (MQ=255)
                 gtcgtATTTGCCGGATTGGCAGTCGATAAAATCGa                  <  1:19336661/35‑1 (MQ=255)
                  tcgtATTTGCCGGATTGGCAGTCGATAAAATCGAt                 <  1:18249637/35‑1 (MQ=255)
                   cgtATTTGCCGGATTGGCAGTCGATAAAATCGAtc                <  1:8910829/35‑1 (MQ=255)
                    gtATTTGCCGGATTGGCAGTCGATAAAATCGAtct               <  1:4485180/35‑1 (MQ=255)
                       tttGCCGGATTGGCAGTCGATAAAATCGATCTcac            <  1:14920974/35‑1 (MQ=255)
                       tttGCCGGATTGGCAGTCGATAAAATCGATCTcac            <  1:6720511/35‑1 (MQ=255)
                       tttGCCGGATTGGCAGTCGATAAAATCGATCTcac            <  1:13569293/35‑1 (MQ=255)
                        ttGCCGGATTGGCAGTCGATAAAATCGATCTcaca           <  1:4482145/35‑1 (MQ=255)
                        ttGCCGGATTGGCAGTCGATAAAATCGATCTcaca           <  1:20106767/35‑1 (MQ=255)
                         tGCCGGGTTGGCAGTCGATAAAATCGATCTCACaa          <  1:7932022/35‑1 (MQ=255)
                         tGCCGGATTGGCAGTCGATAAAATCGATCTCACaa          <  1:8281748/35‑1 (MQ=255)
                         tGCCGGATTGGCAGTCGATAAAATCGATCTCACaa          <  1:2214731/35‑1 (MQ=255)
                         tGCCGGATTGGCAGTCGATAAAATCGATCTCACaa          <  1:3516182/35‑1 (MQ=255)
                         tGCCGGATTGGCAGTCGATAAAATCGATCTCACaa          <  1:11012821/35‑1 (MQ=255)
                          gCCGGATTGGCAGTCGATAAAATCGATTTCACAAg         <  1:7218758/35‑1 (MQ=255)
                          gCCGGATTGGCAGTCGATAAAATCGATCTCACAAg         <  1:8556142/35‑1 (MQ=255)
                           ccGGATTGGCAGTCGATAAAATCGATCTCACAAGc        <  1:9879471/35‑1 (MQ=255)
                           ccGGATTGGCAGTCGATAAAATCGATCTCACAAGc        <  1:7653488/35‑1 (MQ=255)
                           ccGGATTGGCAGTCGATAAAATCGATCTCACAAGc        <  1:16572052/35‑1 (MQ=255)
                           ccGGATTGGCAGTCGATAAAATCGATCTCACAAGc        <  1:19545704/35‑1 (MQ=255)
                            cGGATTGGCAGTCGATAAAATCGATCTCACAAGCa       <  1:19805563/35‑1 (MQ=255)
                            cGGATTGGCAGTCGATAAAATCGATCTCACAAGCa       <  1:2118245/35‑1 (MQ=255)
                            cGGATTGGCAGTCGATAAAATCGATCTCACAAGCa       <  1:2943681/35‑1 (MQ=255)
                             ggATTGGCAGTCGATAAAATCGATCTCACAAGCaa      <  1:14593722/35‑1 (MQ=255)
                             ggATTGGCAGTCGATAAAATCGATCTCACAAGCaa      <  1:13079995/35‑1 (MQ=255)
                              gATTGGCAGTCGATAAAATCGATCTCACAAGCaaa     <  1:11500441/35‑1 (MQ=255)
                               aTTGGCAGTCGATAAAATCGATCTCACAAGCaaaa    <  1:19099620/35‑1 (MQ=255)
                                 tGGCAGTCGATAAAATCGATCTCACAAGCAAAAtt  <  1:4302254/35‑1 (MQ=255)
                                 tGGCAGTCGATAAAATCGATCTCACAAGCAAAAtt  <  1:3084812/35‑1 (MQ=255)
                                  |                                 
ACCCAACTGAATATTACGTCGTATTTGCCGGATTGGCAGTCGATAAAATCGATCTCACAAGCAAAATT  >  REL606/3223825‑3223892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.