New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 403263193 (0.990)3 (0.030) 3/244 NT 3.2% coding (140/615 nt) yigZ UPF0029 family protein
?NC_000913 = 4032664 91 (1.010)coding (173/615 nt) yigZ UPF0029 family protein
Rejected: Frequency below/above cutoff threshold.

GAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4032480‑4032631
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacccGCGACCCACGCCACGCAATGGTGCCGGGCATCGGGGTGTTCTGCCCGCACCGATTCAACAAACGCTTTCGCCGCCTCAACGCCATCGGTATGCGCCAACATCGTTATGAAGCGGCTCTTTTTGATCTCTTCAACGACCGTGACCGG  <  NC_000913/4032664‑4032519
                                                                                                                                                                                                                                                                                                          
GAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACA                                                                                                                                                       >  4:336053/1‑149
 AAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGAGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACAC                                                                                                                                                      >  1:332592/1‑149
  AACTGGCGTGATGGAAAGCTGGGTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACC                                                                                                                                                     >  2:152431/1‑149
                                      ACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAAGCGGTGCGGGCAGATCACCCGCGACCCACGCCACGCAATGGTGCCGGGCATCGGG                                                                                                                 <  1:368972/149‑1
                                        CGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCGCGACCCACGCCACGCAATGGTGCCGGGCATCGGGGT                                                                                                               <  6:273180/149‑1
                                             ACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCGCGACCCACGCCACGCAATGGTGCCGGGCATCGGGGTGTTCT                                                                                                          <  6:80616/149‑1
                                                                                                                                                    ACCCGCGACCCACGCCACGCAATGGTGCCGGGCATCGGGGTGTTCTGCCCGCACCGATTCAACAAACGCTTTCGCCGCCTCAACGCCATCGGTATGCGCCAACATCGTTATGAAGCGGCTCTTTTTGATCTCTTCAACGACCGTGACCG   <  4:311009/149‑1
                                                                                                                                                     CCCGAGACCCACGCCACGCAATGGTGCCGGGCATCGGGGTGTTATGCCCGCACCGATTCAACAAACGCTTTCTCCTCCTCAACGCCATCGGTATGCGCCAACATCGTTATGAATCGGCTCTTTTTGATCTCTTCAACGACCGTGCCCGG  >  6:270997/1‑149
                                                                                                                                                                                                                                                                                                          
GAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTTGGCGCATACCGATGGCGTTGAGGCGGCGAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4032480‑4032631
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacccGCGACCCACGCCACGCAATGGTGCCGGGCATCGGGGTGTTCTGCCCGCACCGATTCAACAAACGCTTTCGCCGCCTCAACGCCATCGGTATGCGCCAACATCGTTATGAAGCGGCTCTTTTTGATCTCTTCAACGACCGTGACCGG  <  NC_000913/4032664‑4032519

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.