New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 225281NA (NA)52 (0.250) 19/186 19.3 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.

GTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTG  >  NC_000913/225267‑225387
                          |                                                                                              
gtcaCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCTGTGCTCACACAGATTGTGTGATGAAAATGAg                        >  2:4327289/4‑99 (MQ=0)
                     aatcTGCGGTGGGATCACCTCCTTACCTTAAAGAAGCTTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCtt   <  2:930228/96‑1 (MQ=0)
                      gccggCGGTTGGATCATCTCCGTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTg  <  2:1458419/95‑1 (MQ=1)
                          |                                                                                              
GTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTG  >  NC_000913/225267‑225387

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