New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1407316 =42 (0.850)3 (0.060) 3/260 4.5 7.8% coding (480/1233 nt) dgcM diguanylate cyclase DgcM
?NC_000913 1416631 = 31 (0.670)coding (756/2601 nt) recE Rac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTA  >  NC_000913/1407182‑1407452
                                                                                                                                      |                                                                                                                                        
aaaTTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCgcggc                                                                                                                                     <  2:825419/140‑1 (MQ=255)
aaaTTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCgcggc                                                                                                                                     <  2:43901/140‑1 (MQ=255)
aaaTTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCgcggc                                                                                                                                     <  1:630558/140‑1 (MQ=255)
     cAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGt                                                                                                                                <  2:535766/140‑1 (MQ=255)
                   caccGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGAtt                                                                                                                  <  2:333859/140‑1 (MQ=255)
                      cGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGCTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGt                                                                                                               <  1:474082/140‑1 (MQ=255)
                       gggTAAATGCGAGATCTCATGCCTGCTAGACATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTc                                                                                                              <  1:262103/140‑1 (MQ=255)
                          tAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCAt                                                                                                           <  2:8743/140‑1 (MQ=255)
                                cGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTg                                                                                                     >  2:481600/1‑140 (MQ=255)
                                         atgcatgATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaa                                                                                            <  1:525214/140‑1 (MQ=255)
                                          tgcatgATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaat                                                                                           >  2:761574/1‑140 (MQ=255)
                                              tgATAGGCATGACGCGACGCCCCAGCATATTTATTTCCCCGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaac                                                                                       <  2:250654/140‑1 (MQ=255)
                                                 tAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCaacaac                                                                                    >  2:654326/1‑140 (MQ=255)
                                                   ggCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAacaaca                                                                                   <  1:189045/139‑1 (MQ=255)
                                                   ggCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAacaaca                                                                                   >  2:189045/1‑139 (MQ=255)
                                                   ggCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACAt                                                                                  <  1:781289/140‑1 (MQ=255)
                                                     cATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTg                                                                                >  1:724187/1‑140 (MQ=255)
                                                        gacgcgacgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTg                                                                             <  1:480969/140‑1 (MQ=255)
                                                        gacgcgacgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTg                                                                             <  2:504906/140‑1 (MQ=255)
                                                              acgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCagag                                                                        >  2:708734/1‑139 (MQ=255)
                                                              acgcCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCagag                                                                        <  1:708734/139‑1 (MQ=255)
                                                                 gacaggCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAg                                                                    <  2:86841/135‑1 (MQ=255)
                                                                 cccGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAg                                                                    <  1:481600/140‑1 (MQ=255)
                                                                     aGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACAt                                                                                  <  2:803553/122‑1 (MQ=255)
                                                                     aGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACAt                                                                                  >  1:803553/1‑122 (MQ=255)
                                                                              atttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcg                                                       <  1:654326/140‑1 (MQ=255)
                                                                               tttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcgc                                                      <  1:761574/140‑1 (MQ=255)
                                                                                ttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATAAACAACATTGGTGCGGAGTTATTCAGAAAAAAgcgcg                                                     >  1:315904/1‑140 (MQ=255)
                                                                                ttCCCAGGTATGTTTCTGGCACATAGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATTAACAACAATGGGGCAGAGTTATTCAGAAAAAAgcgcg                                                     >  1:288698/1‑140 (MQ=255)
                                                                                  cccAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcgcg                                                     >  2:282544/1‑138 (MQ=255)
                                                                                  cccAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcgcg                                                     <  1:282544/138‑1 (MQ=255)
                                                                                   ccAGGTATGTTTCTAGCACATCGTTTCATGCTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCAt                                                  >  2:504234/1‑140 (MQ=255)
                                                                                      ggTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATaaa                                               >  2:424866/1‑140 (MQ=255)
                                                                                      ggTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGGGCAGAGTTAGTCAAAAAAAAGCGCGCATaaa                                               >  2:263226/1‑140 (MQ=255)
                                                                                                 ggCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAATATTGGTGCGGAGTTAGGCAGAAAAAAGCGCGCATAAAAACCTTGtttt                                    >  2:237078/1‑140 (MQ=255)
                                                                                                       atcgttCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTgcgc                              >  2:794276/5‑140 (MQ=255)
                                                                                                                  ttATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGc                   >  1:79627/1‑140 (MQ=255)
                                                                                                                    ataACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCg                 >  2:649529/1‑140 (MQ=255)
                                                                                                                           tAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACt          >  1:665348/1‑140 (MQ=255)
                                                                                                                            aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAACGCGCATAGAAAGCTTGTTTTTTGCACTGATAATTTGCCGAGCGACTg         >  2:219117/1‑140 (MQ=255)
                                                                                                                             gAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTgg        <  1:794276/140‑1 (MQ=255)
                                                                                                                                aTTGAGCGCGGCGAGGTTAGCATCGACGATTCGTCTATCCCGTGACGGTTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCtt     <  2:315904/140‑1 (MQ=255)
                                                                                                                                   gatcgcgGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTa  >  1:467342/4‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
AAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTA  >  NC_000913/1407182‑1407452

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