New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3995016 =82 (1.240)3 (0.050) 3/150 4.5 4.0% coding (231/951 nt) corA magnesium/cobalt transporter CorA
?NZ_CP009273 = 3995012 71 (1.200)coding (227/951 nt) corA magnesium/cobalt transporter CorA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCAGTTCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAGAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3995113‑3995016
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagaagaagGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  NZ_CP009273/3995012‑3994924
                                                                                                                                                                                           
GCAGTTCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAA                                                                                                   <  2:937625/90‑1
GCAGTTCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAA                                                                                                   <  2:1633092/90‑1
   GTTCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGA                                                                                                >  1:1526561/1‑90
    TTCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAA                                                                                               <  1:9741/90‑1
     TCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAG                                                                                              >  1:20207/1‑90
      CACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAGA                                                                                             <  2:204649/90‑1
                                  ATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAGAAGGAGTGAATATGCAGGCCGTCGTCGTC                                                                 >  2:352215/1‑90
                                    CACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTT                                                               <  1:352215/90‑1
                                           GGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAA                                                        >  1:396262/1‑90
                                                  GCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAA                                                                                               >  2:997084/1‑44
                                                  GCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAA                                                                                               <  1:997084/44‑1
                                                                                          AGAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGC         >  1:782880/1‑90
                                                                                          AGAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGC         <  2:1386240/90‑1
                                                                                           GAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCT        >  1:286878/1‑90
                                                                                           GAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCT        >  2:131616/1‑90
                                                                                           GAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCT        >  1:1470571/1‑90
                                                                                           GAAGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCT        >  2:875824/1‑90
                                                                                             AGAAGGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCT      <  1:143872/90‑1
                                                                                                 GGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  1:222821/90‑1
                                                                                                 GGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  2:1274419/90‑1
                                                                                                 GGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  1:208767/90‑1
                                                                                                 GGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  2:1790067/90‑1
                                                                                                 GGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  1:1887660/90‑1
                                                                                                 GGAGGGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  1:1426095/90‑1
                                                                                                                                                                                           
GCAGTTCACGCTCACGCAGAGTAAACAGACGACCATCACGGATGGTAAATGCCACAGTGGAGTTACCGGCGTGATCTTCCGCATCTTCAAAGAAGAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3995113‑3995016
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagaagaagGAGTGAATATGCAGGCCGTCGTCGTCTTCAAAGAAACGTGCCGATGCTTCGATGTCTTCCAGTTCCGGGCGGGTTGCCAGGCTCTGGCC  <  NZ_CP009273/3995012‑3994924

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