New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1362001 =19 (0.860)3 (0.160) 3/114 1.8 14.2% coding (91/558 nt) puuR repressor for the divergent puu operons, putrescine inducible
?NC_000913 1362051 = 20 (1.050)coding (141/558 nt) puuR repressor for the divergent puu operons, putrescine inducible
Rejected: Frequency below/above cutoff threshold.

TGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTT  >  NC_000913/1361985‑1362110
                                                                  |                                                           
cagcgtactgctggcaGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTg                                                            <  1:363307‑M1/52‑1 (MQ=255)
   cgtactgatggcagtgctgactttatcttgttctatcgtactgatagcactgtgagtcagcccTGCCa                                                         <  2:645979‑M2/5‑1 (MQ=255)
           tggcaGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCACTGCCATCAGTACg                                                 <  2:645979‑M1/63‑1 (MQ=255)
             tcCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCt                                               <  2:777370/68‑1 (MQ=255)
              cagggctgactttatcttgttctatcgtactgatagcactgtgagtcagcccTGCCAGCAGTACGCTg                                              <  1:363307‑M2/16‑1 (MQ=255)
               cGGTCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGc                                             <  1:644636/68‑1 (MQ=255)
                   CTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGcaaag                                         <  2:520147‑M1/68‑2 (MQ=255)
                       ctttatcttgttctatcgtactgatagcactgtgagtcagcccTGCCATCAGTACGCTGCAAAAgctg                                     <  1:770795‑M2/25‑1 (MQ=255)
                            aGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGaa                                >  2:523430/1‑68 (MQ=255)
                            aGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGaa                                >  2:483306/1‑68 (MQ=255)
                            aGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGaa                                >  2:1813/1‑68 (MQ=255)
                            aGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGaa                                >  2:785/1‑68 (MQ=255)
                            aGTGCTATCAGAACCATAGAACAAGATAAAGTCAGCCCTGCCATCAGTCCGCTGCAAAAGCTGCTGaa                                >  1:39128/1‑68 (MQ=255)
                              tGCTATCAGTACGATAGAACAAGTTAAAGTCAGCCCTGCCATCAGTACGCTGCTATCGCTGCTGAAgg                              <  1:576901/68‑1 (MQ=255)
                              tGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAgg                              <  1:804599/68‑1 (MQ=255)
                              tGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAgg                              <  1:576899/68‑1 (MQ=255)
                              tGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAgg                              <  1:308574/68‑1 (MQ=255)
                                  aTCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTa                          >  2:214579/1‑68 (MQ=255)
                                  aTCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTa                          >  2:638245/1‑68 (MQ=255)
                                  aTCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTa                          >  2:387397/1‑68 (MQ=255)
                                    cAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATg                        >  1:690400/1‑68 (MQ=255)
                                      gTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGt                      <  2:96953/68‑1 (MQ=255)
                                      gTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGt                      <  2:367667/68‑1 (MQ=255)
                                             aGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCAc               >  1:554980/1‑68 (MQ=255)
                                                     aTAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGaa       >  2:244763/1‑68 (MQ=255)
                                                        aaGTCAGCCCGGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTc    <  2:244875/68‑1 (MQ=255)
                                                          gtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGATGGTGTATGGTCTGTCACTCTCGGAATTCtt  >  2:526783‑M2/9‑68 (MQ=255)
                                                          gtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCtt  >  2:264662‑M2/9‑68 (MQ=255)
                                                          gtcagcccTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCtt  >  1:476450‑M2/9‑68 (MQ=255)
                                                                  |                                                           
TGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCTGCAAAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTT  >  NC_000913/1361985‑1362110

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