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

GAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCT  >  NC_000913/1361936‑1362065
                                                                 |                                                                
gAAATCGCTTGGCGGAAATCCTCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCTAACTCTCCggg                                                                <  2:208373/68‑1 (MQ=255)
gAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCggg                                                                <  2:481530/68‑1 (MQ=255)
     aagaattccgagagtgacagaccatacaccttctgcagcttttgcagcgtactgatggcaGGGCTGAc                                                           >  1:264671‑M1/61‑68 (MQ=255)
     aagaattccgagagtgacagaccatacaccttcagcagcttttgcagcgtactgatggcaGGGCTGAc                                                           >  2:802108‑M1/61‑68 (MQ=255)
         tGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTcac                                                       <  1:163288/68‑1 (MQ=255)
           tCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAg                                                     <  2:708552/68‑1 (MQ=255)
           tCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAg                                                     <  2:451046/68‑1 (MQ=255)
           tCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAg                                                     <  2:216461/68‑1 (MQ=255)
            cGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGt                                                    <  1:14768/68‑1 (MQ=255)
            cGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGt                                                    <  2:711581/68‑1 (MQ=255)
                        atcagcagGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACg                                        <  1:576819/66‑1 (MQ=255)
                            gcagGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGGCTCACAGTGCTATCAGTACGATAg                                    >  1:774287/1‑68 (MQ=255)
                            atacaccttcagcagcttttgcagcgtactgatggcaGGGCTGACTCACAGTGCTATCAGTACGATAg                                    <  2:97965‑M1/31‑1 (MQ=255)
                             cagGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGa                                   <  1:44141/68‑1 (MQ=255)
                                accttcagcagcttttgcagcgtactgctggccGGGCAGACTCGGAGTGCTATCAGTACGAAAGAACa                                <  1:660106‑M1/35‑1 (MQ=255)
                                   cTTTCACAACGTCGTGCCGCCTAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGa                             <  2:100870/68‑1 (MQ=255)
                                   cTTTCACAACGTCGTGCCGCCGAACTCTCCGTGCTGACGCACAGTGCTATCAGTACGATAGAACTAGa                             <  2:291687/68‑1 (MQ=255)
                                   cTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGa                             <  2:386207/68‑1 (MQ=255)
                                             gtcgtGCCGCCGAACTCTCCGGGCTGCCTCACAGTGCTATCACTACGATAGAACAAGATAAAGTCAGc                   >  1:227683/1‑68 (MQ=255)
                                             gtcgtGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGc                   <  2:728983/68‑1 (MQ=255)
                                                  gccgccGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGc              <  2:277318/68‑1 (MQ=255)
                                                  gccgccGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGc              <  2:569411/68‑1 (MQ=255)
                                                  agcgtactgatggcaGGGCTGACTCGCAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGc              <  2:489814‑M1/53‑1 (MQ=255)
                                                    cgccgAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCa            <  1:138833/68‑1 (MQ=255)
                                                              gcaGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCAGCCATCAGTACGCt  <  1:63997‑M1/65‑1 (MQ=255)
                                                                 |                                                                
GAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTCACAGTGCTATCAGTACGATAGAACAAGATAAAGTCAGCCCTGCCATCAGTACGCT  >  NC_000913/1361936‑1362065

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