New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 16888676 (1.100)3 (0.050) 3/100 NT 3.8% coding (188/726 nt) rpsB 30S ribosomal subunit protein S2
?minE = 168904 78 (1.180)coding (206/726 nt) rpsB 30S ribosomal subunit protein S2
Rejected: Frequency below/above cutoff threshold.

ATCATCAACCTTGAGAAAACTGTACCGATGTTCAACGAAGCTCTGGCTGAACTGAACAAGATTGCTTCTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/168816‑168886
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgAAAAGGATTTTACCTTTGCGAGAAGCAATCTTGTTCAGTTCAGCCAGAGCTTCGTTGAACATCGGTACA  <  minE/168904‑168836
                                                                                                                                            
ATCATCAACCTTGAGAAAACTGTACCGATGTTCAACGAAGCTCTGGCTGAACTGAACAAGATTGCTTCTCG                                                                       >  1:556023/1‑71
 TCATCAACCTTGAGAAAACTGTACCGATGTTCAACGAAGCTCTGGCTGAACTGAACAAGATTGCTTCTCG                                                                       >  1:3625146/1‑70
    TCAACCTTGAGAAAACTGTACCGATGTTCAACGAAGCTCTCGCTGAACTGAACAAGATTGCTTCTCG                                                                       >  1:2687395/1‑67
                              TTCAACGAAGCTCTGGCTGAACTGAACAAGATTGCTTCTCGAAAAGGATTTTACCTTT                                                      <  1:484054/58‑1
                                CAACGAAGCTCTGGCTGAACTGAACAAGATTGCTTCTCGAAAAGGATTTTACCTTTGC                                                    >  1:1975955/1‑58
                                             GCTGAACTGAACAAGATTGCTTCTCGAAAAGGATTTTACCTTTGCG                                                   >  1:855929/1‑46
                                                                     CGAAAAGGATTTTACCTTTGCGAGAAGCAATCTTGTTCAGTTCAGCCAGAGCTTCGTTGAACATCGGTACA  <  1:2899737/71‑1
                                                                     CGAAAAGGATTTTACCTTTGCGAGAAGCAATCTTGTTCAGTTCAGCCAGAGCTTCGTTGAACATCGGTACA  <  1:3369955/71‑1
                                                                     CGAAAAGGATTTTACCTTTGCGAGAAGCAATCTTGTTCAGTTCAGCCAGAGCTTCGTTGAACATCGGTACA  <  1:91947/71‑1
                                                                     CGAAAAGGATTTTACCTTTGCGAGAAGCAATCTTGTTCAGTTCAGCC                          <  1:2994815/47‑1
                                                                                                                                            
ATCATCAACCTTGAGAAAACTGTACCGATGTTCAACGAAGCTCTGGCTGAACTGAACAAGATTGCTTCTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/168816‑168886
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgAAAAGGATTTTACCTTTGCGAGAAGCAATCTTGTTCAGTTCAGCCAGAGCTTCGTTGAACATCGGTACA  <  minE/168904‑168836

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