New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 214416442 (0.880)4 (0.090) 3/264 7.2 9.2% coding (899/3318 nt) yegE putative diguanylate cyclase
?NC_000913 = 2144173 40 (0.890)coding (908/3318 nt) yegE putative diguanylate cyclase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAAC  >  NC_000913/2144031‑2144313
                                                                                                                                              |                                                                                                                                            
ctGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAAAGCAAACACATATCCGAAAGCGAAACCCa                                                                                                                                              >  2:7072/1‑143 (MQ=255)
     tgctACGCGGGGTATGTACCTGAAGAGTCATGTGCCGTGGTTGCTGTTTTCGCCAAACCAACAAACCATTAATTTTATTACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACTTTTCCGAAAGCGAAACCCACTTTc                                                                                                                                         <  2:185529/143‑1 (MQ=255)
        tACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAGGCGAGACCCGCTTTCGGa                                                                                                                                      >  1:18095/1‑143 (MQ=255)
         aggccgcGTGCGGACCGGCTGAGCCATATGCCGGGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACAGCATGACCATGGGGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGaa                                                                                                                                     <  2:154169/141‑1 (MQ=255)
               cgtacgtaCCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGAt                                                                                                                               <  1:3185/143‑1 (MQ=255)
                  aggtacCTGACGAGCCCTATGCCGTGGCTCCCGTTTCTCCGTATTCTGCTGCCCGCCAACATCATGCCCATGGTGATGTATGCCTGTCGTGCGGAACGCACACACATTTCCGAAAGCGAAATCCACTTTCGGAACGCGATGGa                                                                                                                            <  2:59311/141‑1 (MQ=255)
                       ccTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATAtt                                                                                                                       <  1:593087/143‑1 (MQ=255)
                          gatgaGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCg                                                                                                                    >  1:512817/1‑143 (MQ=255)
                                   tatGCCGTCGATACCGTATTCACTTATCCTTTTGCCCGCCACCATCATCACAATATTGTTTTATGCCTCTCGCTTAGCACGCACACACATTTCCGAAAGCGAAATCCTGTTTCGGACCTCGATGTAATGTTCATCTAttgctt                                                                                                           >  2:156430/1‑138 (MQ=255)
                                       ccGTGGCTACCGTTTTTGCTGATCCTGCTGCCTGCCAACATCATGACCATGGTGATGTATACCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACACACTTTCGGAACGCAGTGGCATATTCCGCTAACggactcgc                                                                                                       >  2:358627/1‑137 (MQ=255)
                                         ggggCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTCTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGt                                                                                                     <  1:440869/141‑1 (MQ=255)
                                              tACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACGCACTTTCGGAACGCGATGGAATATTCCGCTAGCGGTATGGCGGTAGTg                                                                                                >  1:130663/1‑143 (MQ=255)
                                                          tGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGa                                                                                    <  1:366246/143‑1 (MQ=255)
                                                             tCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACaa                                                                                 <  2:742755/143‑1 (MQ=255)
                                                             tCCTGCATCCCGCCAACACCCTGACCATGGGGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAATCCAAACCCACGTTCCGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACaa                                                                                 <  1:104382/143‑1 (MQ=255)
                                                              ccTGCTGCCCGCCAACCTCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAAt                                                                                <  1:88898/143‑1 (MQ=255)
                                                               ctgctgCCCGCCAACATCATGACCATGGTGATGTATGGCTTTCGTGCGGAACGCAAACACATTGCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATg                                                                               <  1:858611/143‑1 (MQ=255)
                                                                    gcCCGCCAACTGCATGCCCCTCGGGACGCGTGCCTTTCGTACGGATCGCTAACTCATTTCCGAAAGCGAGACCCTCTTTCGGATCGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGc                                                                          <  2:97792/143‑1 (MQ=255)
                                                                    gcCCGCCAACGTCAGGACCGTGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGc                                                                          <  1:300149/143‑1 (MQ=255)
                                                                    gcCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCAATCGGTATGGCGTCAGTGGGCCCCGAGGGACAATGGCTGc                                                                          >  2:619465/1‑143 (MQ=255)
                                                                           aaCATCATGACCATGGTGATGTATGCCTGTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACcaa                                                                   <  1:174902/143‑1 (MQ=255)
                                                                                cacgACCATGGTGATGTACGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGGTGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAg                                                              <  2:493865/140‑1 (MQ=255)
                                                                                 atGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACAAACAAAgc                                                             >  2:412109/1‑143 (MQ=255)
                                                                                  tGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACCATGGCTGCAAACCAACAAAgcg                                                            >  1:68484/1‑143 (MQ=255)
                                                                                   gACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAgcgc                                                           >  1:669132/1‑143 (MQ=255)
                                                                                      cATGGGTATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGctct                                                        <  2:694228/143‑1 (MQ=255)
                                                                                                tATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGGGGGCCCCTGCCTTCTGCCCAACTAAGCGCTCTTCCAATTtct                                              >  1:236599/1‑143 (MQ=255)
                                                                                                   gCCTTTCGGGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCgg                                           <  1:619465/143‑1 (MQ=255)
                                                                                                      tttCGTGCGGAACGCAAACACATTTCCGAAGGCGAAACCCACTTTCGGAACGCTCTGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGGGCTCTGCCAGTTTCTCgggca                                        >  2:342933/1‑141 (MQ=255)
                                                                                                      tttCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTa                                        >  1:346158/1‑143 (MQ=255)
                                                                                                      gcgcGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTa                                        <  1:118925/140‑1 (MQ=255)
                                                                                                               gAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGa                               >  1:760810/1‑143 (MQ=255)
                                                                                                                aaCGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGCTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTATGCCAGTGTCTCGGGGACAGTCAGGaa                              >  2:521057/1‑143 (MQ=255)
                                                                                                                    cgagcacaTGTCCGACAGCGGAACCCACGTTGGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGc                          <  2:883269/139‑1 (MQ=255)
                                                                                                                     aaaCACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCCGTTTCTCGGGTACAGTCAGGACGAGCt                         >  1:769504/1‑143 (MQ=255)
                                                                                                                     aaaCACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGGTTCTCGGGTGCAGGCAGGAAGGGCt                         >  2:512311/1‑143 (MQ=255)
                                                                                                                         acaTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgcgc                     <  1:671859/143‑1 (MQ=255)
                                                                                                                         acaTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTCCCAGTGTCGCGGGTACAGTCAGGAGGACCTgcgc                     >  2:265283/1‑143 (MQ=255)
                                                                                                                           aTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAGGGGCTGCGCgg                   >  1:291470/1‑143 (MQ=255)
                                                                                                                            tttCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGAATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAACCCAACAAAGCGCCCTGCCAGTTTCTCGGGTACAGTCAGGCAGAGCAGCCCGGa                  >  2:621009/1‑143 (MQ=255)
                                                                                                                             ttgccaaaGCGAACCCCACTTTCGGAACGCGATGGAATCTTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGAAAGAGCTGCGCGGAc                 <  1:7072/138‑1 (MQ=255)
                                                                                                                                           cccACTGTCGCAGCGCGCTGGAATATCCCGCTATCGGTATGGCGGTAGTGGGCACCGAGGGACAATGGCGGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTAAGGAAGAGCTGGGCGGACTCACCTTTCAGCaa   <  2:458831/143‑1 (MQ=255)
                                                                                                                                           cccACTCTCCGAACGCGATAGAAGATTCCGGTGTCGGTATGTCGTTAGTGTGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCaa   <  1:119558/143‑1 (MQ=255)
                                                                                                                                            ccACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGGACAGTCAGAAAGAGCTGCGCGGACTGACCTTTCAGCAAc  >  2:606784/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CTGCTTGCTACGCCGCGTACGTACCTGATGAGCCATATGCCGTGGCTACCGTTTTTGCTGATCCTGCTGCCCGCCAACATCATGACCATGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAAC  >  NC_000913/2144031‑2144313

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