New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 228859NA (NA)9 (0.310)
+AG
5/94 2.5 27.3% noncoding (104/120 nt) rrfH 5S ribosomal RNA
?NC_000913 = 1363199 25 (0.820)coding (457/1488 nt) puuC gamma‑glutamyl‑gamma‑aminobutyraldehyde dehydrogenase

AGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/228821‑228859
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TCATCGCCAGCTCAT  <  NC_000913/1363199‑1363185
                                       ||               
AGCGCTGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAG       >  1:226475/1‑50
AGCGCTGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAG       >  1:98314/1‑50
AGCGCTGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCA        <  1:66229/49‑1
AGCGCTGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGC          <  1:2640942/47‑1
    ctGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAGCTCA   >  1:636206/3‑50
    ctGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAGCTCA   >  1:692370/3‑50
     tGATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAGCTCAT  >  1:951916/2‑50
      GATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAGCTCAT  >  1:2443704/1‑49
       ATGGTAGTGTGGGGTTACCCATGTGAGAGTAAGTCATCGCCAGCTCAT  <  1:536444/48‑1
                                       ||               
AGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/228821‑228859
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TCATCGCCAGCTCAT  <  NC_000913/1363199‑1363185

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