New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4596273 =119 (0.880)7 (0.050) 5/168 NT 5.9% coding (284/1023 nt) lgoD L‑galactonate oxidoreductase; L‑gulonate oxidoreductase
?NC_000913 4596294 = 110 (0.860)coding (305/1023 nt) lgoD L‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Rejected: Frequency below cutoff threshold.

CGCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4596376‑4596273
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggcGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCG  >  NC_000913/4596294‑4596389
                                                                                                                                                                                                        
CGCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCC                                                                                                      <  2:3726685/100‑1
 GCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCG                                                                                                     >  1:977776/1‑100
 GCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCG                                                                                                     <  2:435007/100‑1
 GCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCG                                                                                                     >  2:3787236/1‑100
 GCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCG                                                                                                     >  1:346868/1‑100
  CCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGG                                                                                                    >  2:18561/1‑100
  CCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGG                                                                                                    >  1:1280644/1‑100
   CACCGGCACGCTCAAATACTCACTTAAACCTCCATCCTGATGCACGCCAATGACTTAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGGt                                                                                                   <  2:2830069/100‑2
   CACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGGG                                                                                                   <  2:2494616/100‑1
   CACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGGG                                                                                                   <  2:551759/100‑1
   CACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGGG                                                                                                   >  1:1027566/1‑100
                                                                    CACAGCAATTGGTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTA                                  >  1:2620487/1‑100
                                                                           ATTGGTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGC                           >  1:3601863/1‑100
                                                                           ATTGGTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGC                           >  1:924344/1‑100
                                                                             TGGTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCCGGATGGCGGTTTGAGTGAGTATTTGAGCGT                         >  1:1023912/1‑100
                                                                              GGTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTG                        >  1:39670/1‑100
                                                                               GTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGC                       >  2:661273/1‑100
                                                                               GTACGCCCGCTTTTACACGCCGGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGC                       >  2:2695745/1‑100
                                                                                                    GGGCGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCG  <  2:1428796/100‑1
                                                                                                                                                                                                        
CGCCACCGGCACGCTCAAATACTCACTAAAACCGCCATCCTGATGCACGCCAATGACTGAAATTTTTTCACAGCAATTGGTACGCCCGCTTTTACACGCCGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4596376‑4596273
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggcGTACCAATTGCTGTGAAAAAATTTCAGTCATTGGCGTGCATCAGGATGGCGGTTTTAGTGAGTATTTGAGCGTGCCGGTGGCGAACATTTTGCCCG  >  NC_000913/4596294‑4596389

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