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.

CAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAA  >  NZ_CP009273/4278933‑4279108
                                                                                       |                                                                                        
tcaggccaatggcctgctgggctttctctttggttgagatatccgggatctggatttcatgggtgatgcctttggtcgccagcaggcCAt                                                                                        <  1:208438‑M2/3‑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)
            CCTGCTGGCGACCAAAGGCATCATCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACAt                                                                            >  1:226858‑M1/1‑90 (MQ=255)
            CCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACAt                                                                            >  1:151323‑M1/1‑90 (MQ=255)
             CTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATg                                                                           <  2:122669‑M1/90‑1 (MQ=255)
                    cGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACaaa                                                                    <  2:4334/90‑1 (MQ=255)
                         ctctttggttgagatatccgggatctggatttcatgggtgatgcctttggtcgccagcaggcCATTGGCCTGAACATGGAACAAATCAAg                                                               <  2:40407‑M2/28‑1 (MQ=255)
                              cATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCga                                                          <  1:455408/90‑1 (MQ=255)
                                      aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAgg                                                  >  1:252203/1‑90 (MQ=255)
                                      aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAgg                                                  >  2:637674/1‑90 (MQ=255)
                                         tccgggatctggatttcatgggtgatgcctttggtcgccagcaggcCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACt                                               <  2:86590‑M2/44‑1 (MQ=255)
                                            tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCCGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtca                                            >  2:72959/1‑90 (MQ=255)
                                             ccAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcat                                           >  1:419733/1‑90 (MQ=255)
                                             ccAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcat                                           >  1:387590/1‑90 (MQ=255)
                                               aGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcatca                                         >  1:443445/1‑90 (MQ=255)
                                                gATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                   >  2:417041/1‑47 (MQ=255)
                                                gATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                   <  1:417041/47‑1 (MQ=255)
                                                gATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCaa                                        >  2:237721/1‑90 (MQ=255)
                                                       gATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTg                                 <  1:127730/90‑1 (MQ=255)
                                                                aCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAg                        <  2:537349/90‑1 (MQ=255)
                                                                aCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAg                        <  1:131562/90‑1 (MQ=255)
                                                                aCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAg                        <  1:65310/90‑1 (MQ=255)
                                                                  cAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTg                      >  2:482757/1‑90 (MQ=255)
                                                                   aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                     >  1:538220/1‑90 (MQ=255)
                                                                          aaGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAAc              >  1:250430/1‑90 (MQ=255)
                                                                              ccagcaggcCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGc          <  1:637674‑M2/81‑1 (MQ=255)
                                                                              ccagcaggcCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGc          <  1:406709‑M2/81‑1 (MQ=255)
                                                                              ccagcaggcCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGc          <  2:79040‑M2/81‑1 (MQ=255)
                                                                                agcaggcCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCAc        >  2:174997‑M2/8‑90 (MQ=255)
                                                                                      cCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaaa  >  1:639317‑M2/2‑90 (MQ=255)
                                                                                       |                                                                                        
CAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAA  >  NZ_CP009273/4278933‑4279108

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