New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 225281NA (NA)10 (0.150) 8/272 7.6 NA noncoding (1511/1542 nt) rrsH 16S ribosomal RNA of rrnH operon
?NC_000913 = 225293 NA (NA)noncoding (1523/1542 nt) rrsH 16S ribosomal RNA of rrnH operon
Rejected: Coverage evenness skew score above cutoff.

GTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCC  >  NC_000913/225205‑225412
                                                                            |                                                                                                                                   
gtcggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGa                                                                   <  1:254804/140‑1 (MQ=1)
        cgtatccTTAGGGAGGTCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCtga                                                           <  2:630836/138‑1 (MQ=0)
               ttggggAGGGCGCTTACCACTGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAAt                                                    <  2:1063947/140‑1 (MQ=1)
                           ctgaCCACTTTGTGATTCATGCCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTGCTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAAc                                        <  1:1133268/140‑1 (MQ=0)
                                      atgcgtcgtGACTGCGGTCAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCCGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAACCCTCTACAGGCt                             <  2:1060593/135‑1 (MQ=1)
                                                     ggggAAGGCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAggtggt              <  1:32112/140‑1 (MQ=1)
                                                                 gcagggTAACCGGCGGGGCACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACCGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCAccc  <  2:255371/139‑1 (MQ=1)
                                                                            |                                                                                                                                   
GTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCC  >  NC_000913/225205‑225412

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.