New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 22705220 (0.380)10 (0.200) 4/272 7.3 25.9% noncoding (1294/2904 nt) rrlH 23S ribosomal RNA of rrnH operon
?NC_000913 = 4038792 38 (0.750)noncoding (1274/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

ACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038652‑4038867
                                                                                                                                             |                                                                           
aCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAg                                                                           <  1:259885/143‑1 (MQ=255)
 cGGGGATAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAAGAGAGCGGTTTGGAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAGAAGTAACGATAAAGc                                                                          >  1:36646/1‑143 (MQ=14)
  ggggCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCg                                                                         >  1:150335/1‑143 (MQ=255)
   ggggTAAACCATGCACCGAGGCTGCGGGAGCGGCACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  1:888431/143‑1 (MQ=255)
   gggCTAAACGATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  2:952801/143‑1 (MQ=255)
   gggCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  1:371363/143‑1 (MQ=255)
     gCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGt                                                                      <  1:716605/143‑1 (MQ=255)
      cTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTg                                                                     >  1:122758/1‑143 (MQ=255)
      cTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTg                                                                     >  1:295163/1‑143 (MQ=255)
         aaCCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACTTAAGTAACGATAAAGCGGGTGaaa                                                                  >  2:637643/1‑143 (MQ=255)
           ccATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGCTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAg                                                                >  1:403546/1‑143 (MQ=255)
             aTGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCGTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAGGcc                                                              >  1:888681/1‑143 (MQ=255)
              ggCACGGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGccc                                                             <  1:250943/142‑1 (MQ=255)
               gCACCGAAGCTGCGGCAGCGACACTGTGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCg                                                            <  1:531112/143‑1 (MQ=255)
               gCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGACTGCTTACATGAGTAACGATAAAGCGGGTGAAAAGCCCg                                                            >  1:274050/1‑143 (MQ=255)
                  ccGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTc                                                         >  2:232298/1‑143 (MQ=255)
                        cttcGGTAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                   <  2:753486/140‑1 (MQ=255)
                        cTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                   <  2:574529/143‑1 (MQ=255)
                         tGCGGGAGCGACA‑TTTGTGGTGGTTGGGTGGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                  <  2:731603/143‑1 (MQ=11)
                         tGCGGCAGCGAGACGCGGT‑GTTGTTGGGTAGGGGAGCGTGCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                  <  1:164127/143‑1 (MQ=11)
                              cAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATCATGGCACCGCCATATGTGCGCCTCGATACGTAAGTAACGATAAAGCGGGTGGAGAGCCCGCTCGCCGGAAGACCa                                             >  2:423184/1‑143 (MQ=9)
                               aGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa                                            <  2:750844/143‑1 (MQ=255)
                                 cGACGCTATGG‑GGTGTTTGGTAGGGGAGCGTTCTGGAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                          <  2:930706/143‑1 (MQ=11)
                                 cGACACTATGG‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                          <  1:328327/143‑1 (MQ=38)
                                  gACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGCTAAAGCGGGTGAAAAGCCCGCTCGCCGCCAGACCAAggg                                         >  1:897296/1‑143 (MQ=255)
                                   acacTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGt                                        >  1:925933/1‑143 (MQ=38)
                                     acTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTCTGAACGTGTGCTGTGAGGCATGCTGGAGGCACCAGAAGCGCGAATCCCTACATAAGTAACGATAAAGCGCGCAAGAAGCCCCCCCCCCGCAACACCAAgcgccc                                      >  2:109310/1‑138 (MQ=2)
                                       tATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCACGGGTGCCt                                    >  2:206065/1‑143 (MQ=25)
                                             gttgttGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGCGGTCTCCGAACTCCGAATGCTGACATAAGTAACGATAAAGCGGGTGACAAGCCCGCTCGCCGGAAGCCCAAGGGATCCTGTCCa                               >  1:342071/1‑143 (MQ=11)
                                               tgttgGTGAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACTATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAAc                             <  1:107293/143‑1 (MQ=12)
                                               tggtgGGTAGGGGGGCGTTCTGTAGGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAAc                             <  2:605705/140‑1 (MQ=12)
                                                gtgtggTAGGGGAGCGGTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACg                            <  2:31980/139‑1 (MQ=17)
                                                 ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAACGGTTCCTGTCCAACGt                           >  2:784254/1‑143 (MQ=17)
                                                       aGGGGAGCGTTCTGTAAGCTCGTGAAGGTGTGCTGTGAAGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGACAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCg                     <  1:660471/143‑1 (MQ=255)
                                                       aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCg                     >  1:928085/1‑143 (MQ=32)
                                                         gggAGCGTTCTGTGAGCGTGTGATGGTGTGCTGTGAGTCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                   <  2:567319/143‑1 (MQ=255)
                                                         gggAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                   <  1:609341/143‑1 (MQ=32)
                                                              cggtCTGGAAGCTTGTGACGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAgg              <  2:31806/140‑1 (MQ=255)
                                                              cGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAgg              >  2:318567/1‑143 (MQ=32)
                                                               gtgcTGTAAGCCTGTGAACGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAggg             <  2:667624/140‑1 (MQ=17)
                                                                 taggTAAGGCTGTGAAGGGGTGCGGTGAGGCATGCGGGAGGGAGCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg           <  1:775359/140‑1 (MQ=11)
                                                                        gCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCCACGTGAATCGGGGCAGGGTGAGTCGAc    >  2:915539/1‑143 (MQ=17)
                                                                        gCCTGTGAAGGTGAGCTCTTTTGCATTCTGGAGGTAGCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAc    <  2:150335/143‑1 (MQ=11)
                                                                          cTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAccc  <  1:318567/143‑1 (MQ=32)
                                                                                                                                             |                                                                           
ACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGT‑GTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038652‑4038867

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