New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 349144039 (0.970)4 (0.110) 4/152 2.5 9.5% coding (254/1374 nt) cysG siroheme synthase CysG
?NZ_CP009273 = 3491445 41 (1.140)coding (259/1374 nt) cysG siroheme synthase CysG
Rejected: Frequency below/above cutoff threshold.

AGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3491345‑3491440
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgacgcgctGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAGAAGGCTTTCATCAAATGGCCCTTCGACGAGGGTTAAC  <  NZ_CP009273/3491445‑3491357
                                                                                                                                                                                         
AGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGA                                                                                                 >  1:483870/1‑90
   TGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGC                                                                                              <  1:655217/90‑1
                             GGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCCACGGATGATGACGCGCTGGGTAAGGGCGGCATCATCCGTC                                                                    >  2:191460/1‑90
                                          AAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGGTGATGACGCGCTGGGTAAGCGCGTCATCATCCGTCGCTGCAATCGCCA                                                       >  2:239716/1‑90
                                                          TGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAG                                       >  2:510604/1‑90
                                                             TGGCTGGCGATTGCAGCGACGGATGATGACGCGCTGGGTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGGGTCGAGAAG                                    >  1:316920/1‑90
                                                                                       TGACGCGCTGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAGAAGGCTTTCATCAAATGGCCCTTCGACGA          <  1:954568/90‑1
                                                                                        GACGCGCTGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAGAAGGCTTTCATCAAATGGCCCTTCGACGAG         <  2:207923/90‑1
                                                                                         ACGCGCTGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAGAAGGCTTTCATCAAATGGCCCTTCGACGAGG        >  1:457490/1‑90
                                                                                          CGCGCTGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGG                                              <  1:325145/51‑1
                                                                                          CGCGCTGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGG                                              >  2:325145/1‑51
                                                                                               TGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAGAAGGCTTTCATCAAATGGCCCTTCGACGAGGGTTAAC  <  1:630495/90‑1
                                                                                                                                                                                         
AGATGCAGGCATGTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3491345‑3491440
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgacgcgctGGTTAAGCGCGTCATCATCCGTCGCTGCAATCGCCAGCCAGCAGGTGTCGAGAAGGCTTTCATCAAATGGCCCTTCGACGAGGGTTAAC  <  NZ_CP009273/3491445‑3491357

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