New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 1728325248 (1.010)5 (0.020) 4/180 26.2 0.02 coding (711/792 nt) cobA1 uroporphyrinogen‑III C‑methyltransferase
?NC_014034 1728329 = 209 (0.850)coding (707/792 nt) cobA1 uroporphyrinogen‑III C‑methyltransferase
Rejected: Position hash score below cutoff.

TCGGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGC.......................................  >  NC_014034/1728226‑1728325
....................................................................................................GCGGCGGCGGGCAGATCGGCCAGGGTGGTGGCGATCACG  >  NC_014034/1728329‑1728367
                                                                                                                                           
TCGGAACGAAGGAACGCCTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAAAAGGATGACGGCGGGACCTTCCAGACGCG                                        <  1:1560903/101‑1 (MQ=255)
TCGGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCG                                        >  1:2706521/1‑101 (MQ=255)
 CGGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGC                                       >  1:228783/1‑101 (MQ=255)
 CGGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGC                                       >  1:4581067/1‑101 (MQ=255)
 CGGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGAGGGACCTTCCAGAAGCGC                                       >  2:3851661/1‑101 (MQ=255)
  GGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCG                                      >  2:1463917/1‑101 (MQ=255)
   GAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGG                                     <  2:458313/101‑1 (MQ=255)
   GAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGG                                     >  2:2169994/1‑101 (MQ=255)
   GAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGG                                     >  2:1690857/1‑101 (MQ=255)
   GAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGG                                     <  2:1378303/101‑1 (MQ=255)
   GAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGG                                     >  1:3164355/1‑101 (MQ=255)
    AACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGC                                    <  1:5041584/101‑1 (MQ=255)
     ACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCG                                   >  2:2620973/1‑101 (MQ=255)
     ACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCCTCCAGACGGGCGGCG                                   >  2:589609/1‑101 (MQ=255)
      CGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGG                                  >  2:740099/1‑101 (MQ=255)
      CGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGG                                  >  1:2146539/1‑101 (MQ=255)
      CGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGG                                  >  2:4110010/1‑101 (MQ=255)
       GAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGGGGCGGC                                 >  1:832543/1‑101 (MQ=255)
       GAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGC                                 >  2:999431/1‑101 (MQ=255)
       GAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGC                                 >  1:3189081/1‑101 (MQ=255)
       GAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGC                                 >  2:1181799/1‑101 (MQ=255)
       GAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGC                                 >  1:4659960/1‑101 (MQ=255)
        AAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCG                                >  1:3260188/1‑101 (MQ=255)
        AAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCG                                >  1:3654389/1‑101 (MQ=255)
        AAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCG                                <  2:3583643/101‑1 (MQ=255)
        AAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCG                                >  2:4174500/1‑101 (MQ=255)
        AAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCG                                <  2:4325750/101‑1 (MQ=255)
         AGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGG                               <  1:2652247/101‑1 (MQ=255)
         AGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGG                               >  2:4044867/1‑101 (MQ=255)
         AGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGG                               >  1:636773/1‑101 (MQ=255)
          GGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGGGGG                              >  2:3120710/1‑101 (MQ=255)
          GGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGGG                              >  1:684532/1‑101 (MQ=255)
           GAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGGGg                             >  2:3642347/1‑100 (MQ=255)
                     TATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGGGCAGATCGGCCA                   >  1:737719/1‑101 (MQ=255)
                                      ATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGCGCGGCGGCGGGCAGATCGGCCAGGGTGGTGGCGATCACG  <  2:737630/101‑1 (MQ=255)
                                                                                                                                           
TCGGAACGAAGGAACGACTCATATCACGGCCTCGAGCGATTGGGCGGCGGCCGCGCGCGGGGCGAGACCGAACAGGATGACGGCGGGACCTTCCAGACGC.......................................  >  NC_014034/1728226‑1728325
....................................................................................................GCGGCGGCGGGCAGATCGGCCAGGGTGGTGGCGATCACG  >  NC_014034/1728329‑1728367

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG