New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4596219123 (0.780)4 (0.030) 4/164 NT 3.2% coding (230/1023 nt) lgoD L‑galactonate oxidoreductase; L‑gulonate oxidoreductase
?NC_000913 = 4596231 129 (0.890)coding (242/1023 nt) lgoD L‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Rejected: Frequency below cutoff threshold.

CCGATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4596114‑4596219
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgctgaCCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCCA  <  NC_000913/4596231‑4596133
                                                                                                                                                                                                             
CCGATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATAT                                                                                                           >  1:1276782/1‑100
 CGATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATT                                                                                                          >  2:1838225/1‑100
 CGATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATT                                                                                                          >  1:945104/1‑100
 CGATATTCATGCCTGGGGTGAAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATT                                                                                                          <  1:3753800/100‑1
  GATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTG                                                                                                         <  2:4038335/100‑1
  GATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTG                                                                                                         <  2:3517054/100‑1
  GATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTG                                                                                                         >  1:2072005/1‑100
  GATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTG                                                                                                         <  2:1828293/100‑1
   ATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGC                                                                                                        <  2:2072012/100‑1
   ATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGC                                                                                                        <  1:2065902/100‑1
   ATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGC                                                                                                        <  1:2219754/100‑1
    TATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGTGGAGATTGTTGGGCTGGGTAAAAATATTGCT                                                                                                       <  2:1257253/100‑1
    TATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGTGGAGATTGTTGGGCTGGGTAAAAATATTGCT                                                                                                       <  2:1257125/100‑1
    TATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGTGGAGATTGTTGGGCTGGGTAAAAATATTGCT                                                                                                       <  2:1256988/100‑1
    TATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCT                                                                                                       <  1:4263322/100‑1
     ATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTG                                                                                                      <  2:4401222/100‑1
     ATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTG                                                                                                      >  1:3134958/1‑100
      TTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGA                                                                                                     <  2:3594280/100‑1
                              CATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGACCATTTTTAAGATCAGCAATATTT                                                                             <  1:407349/100‑1
                                 TTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGACCATTTTTAAGATCAGCAATATTTTTA                                                                          <  1:3279472/100‑1
                                  TTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGACCATTTTTAAGATCAGCAATATTTTTAC                                                                         <  1:1174567/100‑1
                                   TTTAGTTATCCTCGTGTTTTAGAACATTAAATTTGTAGGGAGATTTTTGGGCTGGGTAAAAATATTGCTGACCATTTTTAAGATCAGCAATATTTTTACC                                                                        <  1:1925621/100‑1
                                                                                                   TTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGA        >  1:2840598/1‑100
                                                                                                   TTGCTGACCATTTTTAAGATAAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGAAAATAAAAAAAAAATGGTTG         <  1:1857176/99‑1
                                                                                                    TGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGAT       >  1:555170/1‑100
                                                                                                    TGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGAT       >  2:1276786/1‑100
                                                                                                    TGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGAT       >  1:555157/1‑100
                                                                                                    TGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGAT       >  2:3589171/1‑100
                                                                                                    TGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGAT       >  1:4366133/1‑100
                                                                                                     GCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATT      >  1:727881/1‑100
                                                                                                     GCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATT      >  1:4389405/1‑100
                                                                                                      CTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTT     >  2:2521939/1‑100
                                                                                                      CTGACCATGTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTT     <  1:1823988/100‑1
                                                                                                       TGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTC    >  2:1514134/1‑100
                                                                                                       TGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTC    <  1:1190648/100‑1
                                                                                                        GACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCC   >  2:2267919/1‑100
                                                                                                        GACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCC   >  2:1542300/1‑100
                                                                                                        GACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCC   >  1:3703464/1‑100
                                                                                                         ACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCCA  >  2:2331622/1‑100
                                                                                                         ACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCCA  >  2:236542/1‑100
                                                                                                         ACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCCA  <  2:2552324/100‑1
                                                                                                         ACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCCA  >  2:3919931/1‑100
                                                                                                                                                                                                             
CCGATATTCATGCCTGGGGTGGAAATCAACCATTTTTTAGTTATCCACGTGTTTTAGGCCATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4596114‑4596219
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgctgaCCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAATCTCCCCACATATTTCATGGCCTAAAACACGTGGATAACTAAAAAATGGTTGATTTCCA  <  NC_000913/4596231‑4596133

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.