New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 97156963 (1.260)3 (0.070) 3/90 NT 4.5% coding (296/747 nt) kdsB 3‑deoxy‑manno‑octulosonate cytidylyltransferase
?AP009048 = 971584 73 (1.660)coding (311/747 nt) kdsB 3‑deoxy‑manno‑octulosonate cytidylyltransferase
Rejected: Frequency below/above cutoff threshold.

CAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTGCAGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/971494‑971569
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagggATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGACAACTT  <  AP009048/971584‑971515
                                                                                                                                                  
CAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:1197563/1‑71
 AGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:857596/1‑70
 AGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:4494911/1‑70
 AGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:4451652/1‑70
 AGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:4244636/1‑70
 AGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:2247915/1‑70
   GAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTGCAG                                                                          >  1:933738/1‑71
                GGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTG                                                                             >  1:4223344/1‑55
                     AAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTGCAG                                                                          >  1:1820921/1‑53
                            CGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTGCAGGGATCATCGGTTCATCACCCTGCA                                                  <  1:2175102/70‑1
                                AAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTGCAGG                                                                         >  1:1514142/1‑43
                                       CATTCAGCGACGACACGGTGATCGTTAATGTGCAGGGATCATCGGTTCAT                                                           <  1:3964136/50‑1
                                             GCGACGACACGGTGATCGTTAATGTGCAGGGATCATCGGTTCATCACC                                                       <  1:4426214/48‑1
                                                                       CAGGGATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGACA      <  1:86008/71‑1
                                                                        AGGGATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGAC       <  1:1902843/69‑1
                                                                          GGATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGACAACT   <  1:1174618/71‑1
                                                                          GGATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCG                              <  1:2994596/44‑1
                                                                           GATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGACAACTT  <  1:2830998/71‑1
                                                                           GATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGACAACTT  <  1:1471889/71‑1
                                                                                                                                                  
CAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGACACGGTGATCGTTAATGTGCAGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/971494‑971569
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagggATCATCGGTTCATCACCCTGCACATTAACGATCACCGTGTCGTCGCTGAATGCGCATTTTTCGACAACTT  <  AP009048/971584‑971515

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