New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3512437 =0 (0.000)23 (1.190) 18/170 0.1 100% intergenic (‑14/+387) aroK/hofQ shikimate kinase AroK/DNA uptake porin HofQ
?NZ_CP009273 = 3646590 NA (NA)pseudogene (3/379 nt) BW25113_RS18165 hypothetical protein

ATCTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3512511‑3512437
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGT  <  NZ_CP009273/3646590‑3646528
                                                                                                                                          
ATCTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAA                                                  <  1:272192/90‑1
  CTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGT                                                >  2:98697/1‑90
     ATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCC                                             <  1:329314/90‑1
      TGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCT                                            <  1:390968/90‑1
         GGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGAT                                         >  1:253221/1‑90
         GGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGAT                                         <  2:447779/90‑1
          GTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATA                                        >  1:1278/1‑90
          GTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATA                                        >  1:52246/1‑90
          GTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAACCGGGTTTATCATTAACGAATAGGCTTAGGGAAAGCGCGAACAAGTCCCTGATA                                        >  2:196299/1‑90
            TTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATG                                      <  1:72740/90‑1
            TTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATG                                      <  2:458649/90‑1
                 TTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATC                                 <  1:94459/90‑1
                        CCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTG                          >  2:351440/1‑90
                                 CTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGA                 <  1:351440/90‑1
                                           GAAGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAAC       <  1:50085/90‑1
                                             AGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAG     <  2:387776/90‑1
                                             AGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAG     <  2:69803/90‑1
                                             AGCGGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAG     >  1:11388/1‑90
                                                GGGTTTATCATTAACGAATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGT  <  1:308977/90‑1
                                                                 ATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGC               >  1:314169/1‑60
                                                                 ATAGTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGC               <  2:314169/60‑1
                                                                    GTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGT  >  1:529306/1‑70
                                                                    GTCTTAGGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGT  <  2:529306/70‑1
                                                                                                                                          
ATCTTATGAGGTTTCAGTTCATGTCCTGCGGCGCTCTCTGAGCGAAGCGGGTTTATCATTAACGAATAGTCTTAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3512511‑3512437
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGT  <  NZ_CP009273/3646590‑3646528

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