New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3946688149 (1.000)111 (0.790) 21/184 11.9 44.2% intergenic (+81/‑12) rrlC/rrfC 23S ribosomal RNA of rrnC operon/5S ribosomal RNA of rrnC operon
?NC_000913 = 4171601 NA (NA)intergenic (+57/‑36) rrlB/rrfB 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
Rejected: Coverage evenness skew score above cutoff.

ACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAC  >  NC_000913/4171513‑4171666
                                                                                        |                                                                 
acgggTACTTATTATCCGGTAGGCGTAACCGTACATCGCCGAAGCTGTTGTGTCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGc                                                         <  1:3958869/96‑1 (MQ=1)
    gTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGACCGcttct                                                     >  2:2117778/1‑95 (MQ=0)
     tattAATGAACCGTCAGGCTGAACCTTACAACGCCGAAGCTGTTTTTGCGGATGAGCGACGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAg                                                    <  2:1667903/96‑1 (MQ=1)
                                     gCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                    >  1:1059856/1‑99 (MQ=31)
                                     gCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                    <  2:1501220/99‑1 (MQ=31)
                                     gCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                    >  2:3039290/1‑99 (MQ=31)
                                     gCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                    >  1:2202724/1‑99 (MQ=31)
                                     gCCGAAGATGTTTTGGCGGATGAGAGACGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                    <  2:2105895/99‑1 (MQ=9)
                                      ccGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                   <  1:491864/99‑1 (MQ=32)
                                      ccGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                   <  1:201448/99‑1 (MQ=32)
                                      ccGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                   <  1:2409457/99‑1 (MQ=32)
                                        gAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa                 <  1:2536901/99‑1 (MQ=33)
                                        gAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa                 <  2:1841461/99‑1 (MQ=33)
                                          aGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgc               >  1:288847/1‑99 (MQ=33)
                                           gCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcg              >  1:564746/1‑99 (MQ=34)
                                             tGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg            <  1:1092052/99‑1 (MQ=34)
                                               ttttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt          >  2:1855621/1‑99 (MQ=33)
                                               ttttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt          >  1:1805645/1‑99 (MQ=34)
                                                tttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg         <  1:2259585/99‑1 (MQ=34)
                                                tttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg         <  1:2177963/99‑1 (MQ=33)
                                                tttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg         <  1:409051/99‑1 (MQ=34)
                                                tttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg         <  1:2122494/99‑1 (MQ=34)
                                                tttGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg         <  1:1468890/99‑1 (MQ=34)
                                                  tGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt       <  2:303291/99‑1 (MQ=34)
                                                  tGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt       <  2:3506845/99‑1 (MQ=34)
                                                  tGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt       <  1:485239/99‑1 (MQ=34)
                                                  tGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt       <  1:2068778/99‑1 (MQ=34)
                                                  tGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt       <  1:1212828/99‑1 (MQ=34)
                                                    gCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAACGCGGTGGTcc     >  2:2466102/1‑99 (MQ=18)
                                                    gCGGACGAGAGAAGATTTTCAGCATGATACAGATTAAATCCGAACGAAGAAACGGTATGATAAAACAGAATTTGCCTGGCGGAAGTAGCTGGATGGTca     >  2:3325806/1‑98 (MQ=11)
                                                      ggATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCa   >  1:2538529/1‑99 (MQ=33)
                                                       gATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc  <  1:3901766/99‑1 (MQ=32)
                                                       gATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc  <  1:2718921/99‑1 (MQ=32)
                                                       gATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc  <  1:2628602/99‑1 (MQ=32)
                                                                                        |                                                                 
ACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAC  >  NC_000913/4171513‑4171666

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