New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1362001 =20 (1.000)6 (0.350) 4/114 1.4 28.4% coding (91/558 nt) puuR repressor for the divergent puu operons, putrescine inducible
?NC_000913 1362051 = 13 (0.760)coding (141/558 nt) puuR repressor for the divergent puu operons, putrescine inducible

GCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCG  >  NC_000913/1361986‑1362115
                                                                 |                                                                
gccgccGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGc                                                                <  1:1678/68‑1 (MQ=255)
agcgtactgatggcaGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGc                                                                <  1:465013‑M1/53‑1 (MQ=255)
  cgtactgatggcagggctgactttatcttgttctatcgtactgatagcactgtgagccagcccTGCCa                                                              <  2:284174‑M2/5‑1 (MQ=255)
  cgccgAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCa                                                              >  1:137799/1‑68 (MQ=255)
       aaCTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGt                                                         <  2:373550/68‑1 (MQ=255)
          tggcaGGGCTGGCTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACg                                                      <  2:284174‑M1/63‑1 (MQ=255)
            tcCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCt                                                    <  1:622838/68‑1 (MQ=255)
            gcagggctgactttatcttgttctatcgtactgatagcactgtgagtcagcccTGCCATCAGTACGCt                                                    <  1:465013‑M2/15‑1 (MQ=255)
                   TGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAg                                             >  1:234326‑M1/1‑68 (MQ=255)
                   TGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAg                                             >  1:183326‑M1/1‑68 (MQ=255)
                            gTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAg                                    >  2:17687/1‑68 (MQ=255)
                             ttgttctatcgtactgatagcactgtgagtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAgg                                   >  1:526832‑M2/37‑68 (MQ=255)
                             ttgttctatcgtactgatagcactgtgagtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAgg                                   >  1:527499‑M2/37‑68 (MQ=255)
                               cTATCAGTACGATAGATCAAGATAAAGTCAGCCCTGCCATCAGTTCGCTGCTAAAGCTGCTGAAGgtg                                 <  1:33217/68‑1 (MQ=255)
                                 aTCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTa                               >  1:332555/1‑68 (MQ=255)
                                 aTCAGTACGATAGAACAAGATAAAGTCAGAACTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGtata                               >  2:137457/1‑68 (MQ=255)
                                   cAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATg                             >  2:384030/1‑68 (MQ=255)
                                          gtAGGACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTc                      >  1:511534/2‑68 (MQ=255)
                                          aTAGAACAAGATAAAGTCAGCACTGCAATAAGTACGCTGCAAAAGCAGCAGAAGGTGTATGGTCTGTc                      >  2:528532/1‑68 (MQ=255)
                                               aCAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCActct                 >  2:618315/1‑68 (MQ=255)
                                               aCAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCActct                 >  2:188077/1‑68 (MQ=255)
                                                     gtgagtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAAt           <  2:526833‑M2/56‑1 (MQ=255)
                                                     gtgagtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAAt           <  2:527500‑M2/56‑1 (MQ=255)
                                                         gtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCtt       >  1:727200‑M2/9‑68 (MQ=255)
                                                          tcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCttt      >  2:140023‑M2/8‑68 (MQ=255)
                                                          tcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCttt      >  1:476829‑M2/8‑68 (MQ=255)
                                                              cccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCg  >  2:461918‑M2/4‑68 (MQ=255)
                                                                 |                                                                
GCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCG  >  NC_000913/1361986‑1362115

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.