New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4289415203 (1.330)3 (0.020) 3/268 NT 1.5% coding (171/567 nt) nrfB nitrite reductase, formate‑dependent, penta‑heme cytochrome c
?NC_000913 = 4289418 202 (1.370)coding (174/567 nt) nrfB nitrite reductase, formate‑dependent, penta‑heme cytochrome c
Rejected: Frequency below/above cutoff threshold.

GTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4289262‑4289415
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ataccTTCGGTATCTGGTTTATGACAGTCCAGACAGGCGGCATCCGGGTTACGCTGCTGGGTAACTTCGTAGCGATCATCAGATGCCTGCGCCGCAGGGGTAGCGGTAATTCCGTTCAGGCTCCACAACAGGCCCGACGCCAGCACCCCGGCA  <  NC_000913/4289418‑4289271
                                                                                                                                                                                                                                                                                                              
GTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGA                                                                                                                                                          <  2:952207/150‑1
 TTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGAT                                                                                                                                                         >  1:2561558/1‑150
  TATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATTATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATA                                                                                                                                                        <  2:2380458/150‑1
  TATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATA                                                                                                                                                        <  2:1765734/150‑1
  TATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATA                                                                                                                                                        <  2:2380444/150‑1
  TATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATA                                                                                                                                                        <  2:2380560/150‑1
  TATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATA                                                                                                                                                        <  2:467456/150‑1
  TATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATA                                                                                                                                                        <  2:83138/150‑1
   ATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATAC                                                                                                                                                       >  2:1748753/1‑150
    TTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACC                                                                                                                                                      >  1:100309/1‑150
    TTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACC                                                                                                                                                      >  1:1445817/1‑150
    TTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACC                                                                                                                                                      >  1:609451/1‑150
                GTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACCTTCGGTATCTGG                                                                                                                                          <  1:320252/150‑1
                                                                            GCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACCTTCGGTATCTGGTTTATGACA                                                                                                                                 >  1:979888/1‑99
                                                                            GCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACCTTCGGTATCTGGTTTATGACA                                                                                                                                 <  2:979888/99‑1
                                                                                                                                                        CCTTCGGTATCTGGTTTATGACAGTCCAGACAGGCGGCATCCGGGTTACGCTGCTGGGTAACTTCGTAGCGATCATCAGATGCCTGCGCCGCAGGGGTAGCGGTAATTCCGTTCAGGCTCCACAACAGGCCCGACGCCAGCACCCCGGCA  <  2:1209720/150‑1
                                                                                                                                                                                                                                                                                                              
GTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAGCAGCGTAACCCGGATGCCGCCTGTCTGGACTGTCATAAACCAGATACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4289262‑4289415
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ataccTTCGGTATCTGGTTTATGACAGTCCAGACAGGCGGCATCCGGGTTACGCTGCTGGGTAACTTCGTAGCGATCATCAGATGCCTGCGCCGCAGGGGTAGCGGTAATTCCGTTCAGGCTCCACAACAGGCCCGACGCCAGCACCCCGGCA  <  NC_000913/4289418‑4289271

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