New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3079260 =35 (1.350)3 (0.140) 3/112 2.1 10.4% coding (71/444 nt) cmtB putative mannitol‑specific enzyme IIA component of PTS
?NC_000913 3080360 = 22 (1.000)coding (1276/1992 nt) tktA transketolase 1, thiamine triphosphate‑binding
Rejected: Frequency below/above cutoff threshold.

TAATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3079336‑3079260
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcaggaagcCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGT  >  NC_000913/3080360‑3080425
                                                                                                                                               
TAATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAGTCTCAGTTATACATTCAGCAAAAGATTG                                                                             >  2:465674/1‑68
 AATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGG                                                                            <  2:352968/68‑1
    TTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAG                                                                         >  1:729900/1‑68
    TTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAG                                                                         >  2:430360/1‑68
        GCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCTGGAAG                                                                     >  1:218784/1‑68
         CGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC                                                                    <  2:60518/68‑1
         CGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC                                                                    >  2:94663/1‑68
                                     CAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAAT                                        >  1:81549/1‑68
                                               TATACATTCAGCAAAAGATTGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATA                              <  2:81549/68‑1
                                                         GCAAAAGATTGGCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGC                    <  1:226822/68‑1
                                                                    GCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACC             >  1:53966/1‑64
                                                                    GCAGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACC             <  2:53966/64‑1
                                                                      AGGAAGCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTG       >  1:429342/1‑68
                                                                         AAGCCACCGTGCAGGGAGTTACCGTTAGCAATCGCGGTCATACCGAACTCTCGAACTCCGTTGTGGAT    <  1:258585/68‑1
                                                                           GCCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGT  <  2:962479/68‑1
                                                                                                                                               
TAATTTATGCGGCTTAGTGATTATTTTCCAGAATCATCAATCTCAGTTATACATTCAGCAAAAGATTGGCAGGAAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3079336‑3079260
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcaggaagcCACCGTGCAGGGAGATACCGTTAGCAATCGCGGTCATACCGAACTCGCGAACACCGTAGTGGATGT  >  NC_000913/3080360‑3080425

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