New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 442673 =25 (0.710)5 (0.160) 4/146 2.4 18.3% coding (489/891 nt) cyoE heme o synthase
?NZ_CP009273 442734 = 23 (0.750)coding (428/891 nt) cyoE heme o synthase
Rejected: Frequency below/above cutoff threshold.

TCTGTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/442768‑442673
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcgcagcgCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATC  >  NZ_CP009273/442734‑442820
                                                                                                                                                                                       
TCTGTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGC                                                                                               <  2:1082801/90‑1
TCTGTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGC                                                                                               <  2:237893/90‑1
  TGTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGC                                                                                             <  1:149717/90‑1
   GTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCA                                                                                            >  2:942817/1‑90
   GTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCA                                                                                            >  1:796343/1‑90
     CTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGC                                                                                          <  1:93201/90‑1
      TACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCG                                                                                         >  1:829132/1‑90
      TACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCG                                                                                         >  2:113977/1‑90
             CGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAG                                                                                           >  1:818054/1‑81
             CGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAG                                                                                           <  2:818054/81‑1
                            TCTCCGGCGCTGCGCCTCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCGCCGGAGAGCGAACCAATCAACG                                                                   <  2:847836/90‑1
                            TCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCGCCGGAGAGCGAACCAATCAACG                                                                   <  2:530111/90‑1
                                   CGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTA                                                            <  2:130873/90‑1
                                     CTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGA                                                          <  1:1011523/90‑1
                                                                      CCGGTGAGTTCGATAGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGC                         <  2:854691/90‑1
                                                                                        GCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGC       <  2:720616/90‑1
                                                                                          GCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCC     >  2:858417/1‑90
                                                                                          GCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCC     >  2:999841/1‑90
                                                                                             GCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATC  <  2:128206/90‑1
                                                                                             GCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATC  <  1:514445/90‑1
                                                                                                                                                                                       
TCTGTCTACGGCACGTTGATTGGTTCGCTCTCCGGCGCTGCGCCGCCGGTGATCGGCTACTGTGCGGTAACCGGTGAGTTCGATAGCGGCGCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/442768‑442673
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcgcagcgCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATC  >  NZ_CP009273/442734‑442820

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