New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3726887 =125 (1.390)3 (0.040) 3/224 NT 2.7% intergenic (+5/+37) wecH/yiaA O‑acetyltransferase for enterobacterial common antigen (ECA)/YiaAB family inner membrane protein, tandem domains
?NC_000913 3726917 = 100 (1.230)intergenic (+35/+7) wecH/yiaA O‑acetyltransferase for enterobacterial common antigen (ECA)/YiaAB family inner membrane protein, tandem domains
Rejected: Frequency below/above cutoff threshold.

TGGACTGGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATGCCGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3727041‑3726887
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaccgggcttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTG  >  NC_000913/3726917‑3727059
                                                                                                                                                                                                                                                                                                          
tcGACTGGCTTTCTTCTTAAGCGTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATTCCTTAATTAACCCACCAAAAGAAACACAGGTTACCCAGGTAGAATACAGCTATTAACTCACGTAAGCCCGGGCAGTCCAATGTGAC                                                                                                                                                        <  6:324469/146‑1
TGGACTGGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCTGAAAAATATTCGTGATGCAGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGA                                                                                                                                                         <  8:342467/147‑1
  GACTGGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATGCCGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACC                                                                                                                                                       <  3:149945/147‑1
    CTGGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATGCCGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGG                                                                                                                                                    >  3:66274/1‑148
     TGGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATGCCGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGC                                                                                                                                                   >  4:44536/1‑148
      GGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATGCCGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGCT                                                                                                                                                  >  8:302909/1‑148
                                                               CGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTT                                                                                          >  2:60156/1‑147
                                                                        CCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCA                                                                                >  4:135923/1‑148
                                                                                                 CAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTA                                                       >  8:326323/1‑148
                                                                                                                                                 GACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCT                                                        <  1:348024/99‑1
                                                                                                                                                 GACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCT                                                        >  2:348024/1‑99
                                                                                                                                                  ACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTT        <  7:331294/146‑1
                                                                                                                                                      GGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTG  <  8:56811/148‑1
                                                                                                                                                      GGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCT   >  4:320126/1‑147
                                                                                                                                                                                                                                                                                                          
TGGACTGGCTTTCTTCTTAAGCTTGTTTGGTGCAGTAGCGGTGCAGAAAAATATTCGTGATGCCGGAATAAACCCACCAAAAGAAACACAGGTTACCCAGGAAGAATACAGCGAATAACTCACGTAAGCCCGGTCAGTCCAATGTGACCGGGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3727041‑3726887
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaccgggcttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTG  >  NC_000913/3726917‑3727059

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