New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4278944 =30 (1.290)3 (0.140) 3/152 1.8 10.8% coding (1253/1437 nt) nrfA ammonia‑forming nitrite reductase cytochrome c552 subunit
?NZ_CP009273 4279020 = 23 (1.110)coding (1329/1437 nt) nrfA ammonia‑forming nitrite reductase cytochrome c552 subunit
Rejected: Frequency below/above cutoff threshold.

CGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACA  >  NZ_CP009273/4278855‑4279033
                                                                                         |                                                                                         
ttttacgtgcctgttcttcccactgcgggatcaccgttttgatgaagtcctgcttctcggccttgatttgttccatgttcaggccaatgg                                                                                           >  1:639317‑M1/90‑90 (MQ=255)
      gtgcctgttcttcccactgcgggatcaccgttttgatgaagtcctgcttctcggccttgatttgttccatgttcaggccaatgGCctgct                                                                                     >  2:174997‑M1/84‑90 (MQ=255)
        gcctgttcttcccactgcgggatcaccgttttgatgaagtcctgcttctcggccttgatttgttccatgttcaggccaatgGCCTGCTgg                                                                                   <  2:79040‑M1/9‑1 (MQ=255)
        gcctgttcttcccactgcgggatcaccgttttgatgaagtcctgcttctcggccttgatttgttccatgttcaggccaatgGCCTGCTgg                                                                                   <  1:637674‑M1/9‑1 (MQ=255)
        gcctgttcttcccactgcgggatcaccgttttgatgaagtcctgcttctcggccttgatttgttccatgttcaggccaatgGCCTGCTgg                                                                                   <  1:406709‑M1/9‑1 (MQ=255)
          atGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCGCGCACCAAACGGGCGCGCCTGgcggcg                                                                                 >  1:511896/1‑90 (MQ=255)
          atGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCg                                                                                 >  1:188767/1‑90 (MQ=255)
            gcacgcacCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGAc                                                                               >  2:422246/1‑90 (MQ=255)
               cgcacCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCaa                                                                            >  2:612877/1‑90 (MQ=255)
                gcacCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCaaa                                                                           >  2:470048/1‑90 (MQ=255)
                 cacCGGAAGAAGGTTTACGGATGCTCGGTACGTCGATGGATAAAGCGGCGGATGCACACACCAAAACGGCGGGCCGGCTGGCGTCCaaaa                                                                          >  1:549855/1‑89 (MQ=255)
                 cacCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAg                                                                          >  2:430886/1‑90 (MQ=255)
                 cacCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACAGGCGCGCCTGCTGGCGACCCAAg                                                                          >  1:291663/1‑90 (MQ=255)
                   ccGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGc                                                                        >  2:381962/1‑90 (MQ=255)
                    cGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAACCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGTCa                                                                       <  1:470048/90‑1 (MQ=255)
                      gaagaaGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATc                                                                     <  2:3346/90‑1 (MQ=255)
                         gaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCAcc                                                                  <  1:261893/90‑1 (MQ=255)
                           agGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCa                                                                >  2:178469/1‑90 (MQ=255)
                               ttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGaa                                                            <  1:430886/90‑1 (MQ=255)
                               ttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGaa                                                            <  1:633964/90‑1 (MQ=255)
                                            gTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCgg                                               <  2:410794/90‑1 (MQ=255)
                                             tACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGa                                              <  2:317853/90‑1 (MQ=255)
                                             agtcctgcttctcggccttgatttgttccatgttcaggccaatgGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGa                                              <  2:86590‑M1/46‑1 (MQ=255)
                                                 gCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATAtc                                          <  2:89511/90‑1 (MQ=255)
                                                 gCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATAtc                                          <  2:328393/90‑1 (MQ=255)
                                                   gATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATAtctc                                        <  1:383966/90‑1 (MQ=255)
                                                             cttgatttgttccatgttcaggccaatgGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAagag                              <  2:40407‑M1/62‑1 (MQ=255)
                                                                gcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCATCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGaaa                           >  2:226858/1‑90 (MQ=255)
                                                                gcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGaaa                           >  2:151323/1‑90 (MQ=255)
                                                                      gcacccacCAAACTGGCGCGCCTTCTGGCGACCAAAGGCTTCACCCATGAAAACCAGATCCCGGATATCTCAACCAAAAAAAAAACCcag                     >  1:209748/1‑90 (MQ=255)
                                                                       cacgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCcagc                    <  2:78943/90‑1 (MQ=255)
                                                                         cgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCcagcag                  <  2:627530/90‑1 (MQ=255)
                                                                              cAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAt             <  2:188767/90‑1 (MQ=255)
                                                                               aaaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAtt            >  1:136985/1‑90 (MQ=255)
                                                                                aaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTg           <  1:349053/90‑1 (MQ=255)
                                                                                   tGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc        <  2:113492/90‑1 (MQ=255)
                                                                                    ggCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCt       <  2:628975/90‑1 (MQ=255)
                                                                                      atgGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGa     <  1:208438‑M1/87‑1 (MQ=255)
                                                                                         GCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACa  <  2:262627‑M1/90‑1 (MQ=255)
                                                                                         |                                                                                         
CGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACA  >  NZ_CP009273/4278855‑4279033

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