New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4289734141 (1.170)3 (0.030) 3/270 NT 2.2% coding (490/567 nt) nrfB nitrite reductase, formate‑dependent, penta‑heme cytochrome c
?NC_000913 = 4289754 134 (1.150)coding (510/567 nt) nrfB nitrite reductase, formate‑dependent, penta‑heme cytochrome c
Rejected: Frequency below/above cutoff threshold.

GTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4289582‑4289734
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGATGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACATCGTGCGGCCAGAACGCTTTTTGCAACT  <  NC_000913/4289754‑4289606
                                                                                                                                                                                                                                                                                                              
GTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTT                                                                                                                                                          >  2:571412/1‑150
GTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTT                                                                                                                                                          >  2:749630/1‑150
 TATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTG                                                                                                                                                         <  2:536123/150‑1
  ATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGC                                                                                                                                                        >  1:1188089/1‑150
   TGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGCG                                                                                                                                                       <  2:571022/150‑1
                           TGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGCGTCGATCGCTGTGGCAATCGACGCA                                                                                                                               <  1:985752/150‑1
                                    CGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGCGCTGATCGCTGTGGCAATCGACGCAAATCTTAAT                                                                                                                      <  1:1770308/150‑1
                                      TTCTGGCCGCACGATGTCCACGAAACCACCGGGGCGTGCGCCAGCTGCCATTCCCTGCATCCACACCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGCGCTGATCTCTGTGGCAATCGACGCAAATCTTAATCC                                                                                                                    <  1:107553/150‑1
                                                                                                                                                      GCGCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGATGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACATCGTGCGGCCAGAACGCTTTTTGCAA    >  1:1870685/1‑150
                                                                                                                                                      GCGCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGAGGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGCCATCGTGCGGCCAGAACGCTTTTTGCAA    >  2:1265633/1‑150
                                                                                                                                                       CGCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGATGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACATCGTGCGGCCAGAACGCTTTTTGCAAC   <  2:1837177/150‑1
                                                                                                                                                       CGCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGATGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACATCGTGCGGCCAGAACGCTTTTTGCAAC   <  2:747681/150‑1
                                                                                                                                                       CGCTGAGCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTACCGTCGGCCTCGGATCTTGCGGTGGATGCAGGGAAGGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACAGCGTGCGGCCAGAACGCTTTTTGCAAC   <  1:1315020/150‑1
                                                                                                                                                        GCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGATGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACATCGTGCGGCCAGAACGCTTTTTGCAACT  <  2:760091/150‑1
                                                                                                                                                                                                                                                                                                              
GTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTTCTGGCCGCACGATGTCCACGTAACCAAAGTGGCGTGCGCCAGCTGCCATTCCCTGCATCCACAGCAAGATACGATGCAGACGTTAAGTGACAAAGGACGGATTAAGATTTGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4289582‑4289734
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgCTGATCGCTGTGGCAATCGACGCAAATCTTAATCCGTCCTTTGTCACTTAACGTCTGCATCGTATCTTGCTGTGGATGCAGGGAATGGCAGCTGGCGCACGCCACTTTGGTTACGTGGACATCGTGCGGCCAGAACGCTTTTTGCAACT  <  NC_000913/4289754‑4289606

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