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.

ACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038652‑4038867
                                                                                                                                             |                                                                           
aCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAg                                                                           <  1:259885/143‑1 (MQ=255)
 cGGGGATAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAAGAGAGCGGTTTGGAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAGAAGTAACGATAAAGc                                                                          >  1:36646/1‑143 (MQ=14)
  ggggCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCg                                                                         >  1:150335/1‑143 (MQ=255)
   ggggTAAACCATGCACCGAGGCTGCGGGAGCGGCACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  1:888431/143‑1 (MQ=255)
   gggCTAAACGATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  2:952801/143‑1 (MQ=255)
   gggCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  1:371363/143‑1 (MQ=255)
     gCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGt                                                                      <  1:716605/143‑1 (MQ=255)
      cTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAGGTAACGATAAAGCGGGTg                                                                     >  1:122758/1‑143 (MQ=255)
      cTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTg                                                                     >  1:295163/1‑143 (MQ=255)
         aaCCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACTTAAGTAACGATAAAGCGGGTGaaa                                                                  >  2:637643/1‑143 (MQ=255)
           ccATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGCTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAg                                                                >  1:403546/1‑143 (MQ=255)
             aTGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCGTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAGGcc                                                              >  1:888681/1‑143 (MQ=255)
              ggCACGGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGccc                                                             <  1:250943/142‑1 (MQ=255)
               gCACCGAAGCTGCGGCAGCGACACTGTGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCg                                                            <  1:531112/143‑1 (MQ=255)
               gCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGACTGCTTACATGAGTAACGATAAAGCGGGTGAAAAGCCCg                                                            >  1:274050/1‑143 (MQ=255)
                  ccGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTc                                                         >  2:232298/1‑143 (MQ=255)
                        cttcGGTAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                   <  2:753486/140‑1 (MQ=255)
                        cTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                   <  2:574529/143‑1 (MQ=255)
                         tGCGGGAGCGACATTTGTGGTGGTTGGGTGGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                  <  2:731603/143‑1 (MQ=11)
                         tGCGGCAGCGAGACGCGGTGTTGTTGGGTAGGGGAGCGTGCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                  <  1:164127/143‑1 (MQ=11)
                              cAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATCATGGCACCGCCATATGTGCGCCTCGATACGTAAGTAACGATAAAGCGGGTGGAGAGCCCGCTCGCCGGAAGACCa                                             >  2:423184/1‑143 (MQ=9)
                               aGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCaa                                            <  2:750844/143‑1 (MQ=255)
                                 cGACGCTATGGGGTGTTTGGTAGGGGAGCGTTCTGGAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                          <  2:930706/143‑1 (MQ=11)
                                 cGACACTATGGGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                          <  1:328327/143‑1 (MQ=38)
                                  gACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGCTAAAGCGGGTGAAAAGCCCGCTCGCCGCCAGACCAAggg                                         >  1:897296/1‑143 (MQ=255)
                                   acacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGt                                        >  1:925933/1‑143 (MQ=38)
                                     acTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTCTGAACGTGTGCTGTGAGGCATGCTGGAGGCACCAGAAGCGCGAATCCCTACATAAGTAACGATAAAGCGCGCAAGAAGCCCCCCCCCCGCAACACCAAgcgccc                                      >  2:109310/1‑138 (MQ=2)
                                       tATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCACGGGTGCCt                                    >  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)
                                                                                                                                             |                                                                           
ACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  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.