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

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

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.