New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 2258692 =168 (0.880)20 (0.110) 16/188 16.6 10.4% coding (605/942 nt) psuK pseudouridine kinase
?NC_000913_3_bme_pgi 2258736 = 183 (0.990)coding (561/942 nt) psuK pseudouridine kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CATGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/2258794‑2258692
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATT  >  NC_000913_3_bme_pgi/2258736‑2258834
                                                                                                                                                                                                           
CATGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGTGATTGCGCTGTCAGG                                                                                                         >  1:693896/1‑100
CATGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGG                                                                                                         >  2:1091215/1‑100
CATGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGG                                                                                                         <  1:4432838/100‑1
CATGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGG                                                                                                         <  1:4432903/100‑1
  TGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGC                                                                                                       >  1:3429365/1‑100
   GGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGC                                                                                                       >  2:1720880/1‑99
    GAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCG                                                                                                      <  2:2410710/99‑1
                                   CAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGAC                                                                     >  2:4062324/1‑100
                                               CTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTT                                                         >  2:1706399/1‑100
                                                   AGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGA                                                      >  1:647908/1‑99
                                                    GCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACA                                                    >  1:717626/1‑100
                                                    GCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTAGCGCACTTTGACA                                                    >  2:3759246/1‑100
                                                     CCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACA                                                    >  2:1649288/1‑99
                                                     CCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACA                                                    >  2:1662041/1‑99
                                                      CAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACA                                                  >  1:3367729/1‑100
                                                      CAAACCGCCTTGAAGCGGAAAACCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGCGAGTGTG‑GATCTGATTTAGACGGTCGCGCGCTTTGACACA                                                  >  2:4743151/1‑100
                                                       AAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGACTTGAGAGTGTGTGATCTGATTTAGACGGTCGCGCACTTTGACACA                                                  >  2:2089454/1‑100
                                                       AAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGACTTGAGAGTGTGTGATCTGATTTAGACGGTCGCGCACTTTGACACA                                                  >  2:2089409/1‑100
                                                         ACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTT                                               >  1:1847243/1‑100
                                                          CCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGCCGGTCGCGCACTTTGACACATTTC                                              >  2:3626784/1‑100
                                                           CGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTC                                              >  1:589682/1‑99
                                                            GCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCA                                            >  1:1136620/1‑100
                                                             CCTTGAAGCGGAAACCCTGAGTGGGATTGCGCT‑TCAAGGCGGTTTGGCCTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATG                                          >  2:2585685/1‑100
                                                              CTTGAAGCGGAAACCCTGAGTGGGATTGCGCT‑TCAAGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGC                                         >  1:832849/1‑100
                                                                 GAAGCGGAAACCCTGAGTGGGATTGCGCT‑TCAAGGCGGTTTGGCTTGAGAGTGTG‑GATCTGTTCTGGCCGGTCGCGCACTATGCCACATTTCCATGCCG                                       >  2:1643703/1‑99
                                                                       GAAACCCTGAGTGGGATTGCGCT‑TCAAGGCGGTTTGGCTTGAGAGTATG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCG                                 >  1:2029521/1‑99
                                                                                  TGGGATTGCGCTGTCAGGGCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAAT                      <  1:1711689/100‑1
                                                                                                    GCGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCA    <  1:2070208/100‑1
                                                                                                     CGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCAT   <  2:848256/100‑1
                                                                                                     CGGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCAT   <  2:3847985/100‑1
                                                                                                      GGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATT  <  1:1248311/100‑1
                                                                                                      GGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCAT   >  2:515637/1‑99
                                                                                                                                                                                                           
CATGGAAATGTGTCAAAGTGCGCGACCGTCTAAATCAGATCCACACTCTCAAGCCAAACCGCCTTGAAGCGGAAACCCTGAGTGGGATTGCGCTGTCAGGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_bme_pgi/2258794‑2258692
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgGTTTGGCTTGAGAGTGTG‑GATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATT  >  NC_000913_3_bme_pgi/2258736‑2258834

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