New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 34860233 (1.600)3 (0.160) 3/156 1.6 9.1% coding (186/1887 nt) prpE propionate‑‑CoA ligase
?NZ_CP009273 = 348608 30 (1.590)coding (192/1887 nt) prpE propionate‑‑CoA ligase
Rejected: Frequency below/above cutoff threshold.

CAGCAATCCGCCGTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGAC  >  NZ_CP009273/348527‑348689
                                                                                 |                                                                                 
caGCAATCCGCCGTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCaga                                                                           >  1:177160/1‑90 (MQ=255)
caGCAATCCGCCGTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCaga                                                                           >  2:63958/1‑90 (MQ=255)
  gCAATCCGCCGTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAgg                                                                         >  1:275430/1‑90 (MQ=255)
     aTCCGCCGTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAggcgc                                                                      <  1:470832/90‑1 (MQ=255)
        cgccgTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAggcgctgg                                                                   <  1:343044/90‑1 (MQ=255)
         gccgTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAggcgctggc                                                                  <  2:483373/90‑1 (MQ=255)
            gTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAggcgctggcgct                                                               >  1:443664/1‑90 (MQ=255)
             tttGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGgcgctggcgctg                                                              <  2:425313/90‑1 (MQ=255)
              ttGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAAcagccagcggtcgatggcgttg                                                             >  2:307106‑M2/1‑68 (MQ=255)
                               gAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGaa                                            <  2:98658/90‑1 (MQ=255)
                               gAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGaa                                            <  2:444569/90‑1 (MQ=255)
                               gAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGaa                                            <  1:535897/90‑1 (MQ=255)
                                 aGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAAc                                          <  2:182980/90‑1 (MQ=255)
                                 aGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAAc                                          <  2:110156/90‑1 (MQ=255)
                                  ggCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACa                                         >  1:93815/1‑90 (MQ=255)
                                     cGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAAcagccagcggtcgatggcgttgtggcacaagttggttcggccttc                                      >  1:211434‑M2/1‑45 (MQ=255)
                                         ccAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAAcagccagcggtcgatggcgttgtggcacaagttggttcggccttcacaa                                  >  2:549474‑M2/1‑41 (MQ=255)
                                            aCTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgc                               <  1:467907/90‑1 (MQ=255)
                                             cTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcg                              <  2:265287/90‑1 (MQ=255)
                                             cTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcg                              <  2:425868/90‑1 (MQ=255)
                                              ttGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcgc                             <  1:254728/90‑1 (MQ=255)
                                              ttGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcgc                             <  1:348604/90‑1 (MQ=255)
                                                    cacaACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCaccttt                       <  2:177160/90‑1 (MQ=255)
                                                      caACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcgc                             >  1:533200/1‑82 (MQ=255)
                                                      caACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcgc                             <  2:533200/82‑1 (MQ=255)
                                                      caACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCacctttac                     >  2:449306/1‑90 (MQ=255)
                                                      caACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCacctttac                     >  2:418468/1‑90 (MQ=255)
                                                                gACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCg                                               <  2:211546/54‑1 (MQ=255)
                                                                gACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCg                                               >  1:211546/1‑54 (MQ=255)
                                                                gACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAg           >  2:5077/1‑90 (MQ=255)
                                                                       ggctggAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcgc                             >  2:522794/1‑65 (MQ=255)
                                                                       ggctggAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAgcgc                             <  1:522794/65‑1 (MQ=255)
                                                                         ctggAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGAc  <  1:305575/90‑1 (MQ=255)
                                                                                 |                                                                                 
CAGCAATCCGCCGTTTGCCCGTTGGTTTTGTGAAGGCCGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGAC  >  NZ_CP009273/348527‑348689

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