New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4404334 =22 (0.240)4 (0.050) 3/188 NT 15.5% intergenic (+30/‑52) hflC/yjeT HflB protease modulator specific for phage lambda cII repressor/DUF2065 family protein
?NC_000913 4404361 = 23 (0.260)intergenic (+57/‑25) hflC/yjeT HflB protease modulator specific for phage lambda cII repressor/DUF2065 family protein
Rejected: Coverage evenness skew score above cutoff.

TGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATG  >  NC_000913/4404202‑4404465
                                                                                                                                    |                                                                                                                                   
tGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAATGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAg                                                                                                                       >  2:300638/1‑147 (MQ=255)
tGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAATGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAg                                                                                                                       <  1:300638/147‑1 (MQ=255)
    tATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAg                                                                                                                       >  4:783200/1‑143 (MQ=255)
    tATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAg                                                                                                                       <  3:783200/143‑1 (MQ=255)
          aaCAGCTTCTCTGGCAATGAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGcc                                                                                                                                    >  3:906561/1‑124 (MQ=255)
          aaCAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGcc                                                                                                                                    <  4:906561/124‑1 (MQ=255)
                                      tGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCAcc                                                                                                                                 <  3:71833/99‑1 (MQ=255)
                                      tGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCAcc                                                                                                                                 >  4:71833/1‑99 (MQ=255)
                                                 ccGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTc                                                                                                                         >  1:36832/1‑96 (MQ=255)
                                                 ccGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTc                                                                                                                         <  2:36832/96‑1 (MQ=255)
                                                                    gCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCt                                                  >  4:523123/1‑148 (MQ=255)
                                                                          tGAAGACGCCGACTTCCGTAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCgg                                                                                                              >  4:717822/1‑82 (MQ=255)
                                                                          tGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCgg                                                                                                              <  3:717822/82‑1 (MQ=255)
                                                                            aaGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCAc                                                                                                                                  <  1:660616/60‑1 (MQ=255)
                                                                            aaGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCAc                                                                                                                                  >  2:660616/1‑60 (MQ=255)
                                                                            aaGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGa                                                                                                                    <  3:266680/74‑1 (MQ=255)
                                                                            aaGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGa                                                                                                                    >  4:266680/1‑74 (MQ=255)
                                                                                  ccGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTCTTCTCTTTCtata                                                                                            >  3:1078280/1‑92 (MQ=255)
                                                                                  ccGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGAATCCTCAGGGATGTCGGTGGTTTTCTTTTTCtata                                                                                            <  4:1078280/92‑1 (MQ=255)
                                                                                                    tAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATc                                                                   <  4:396338/99‑1 (MQ=255)
                                                                                                    tAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATc                                                                   >  3:396338/1‑99 (MQ=255)
                                                                                                          aaCGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCt                                                        >  3:554516/1‑104 (MQ=255)
                                                                                                          aaCGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCt                                                        <  4:554516/104‑1 (MQ=255)
                                                                                                                    gTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATg  <  3:523123/148‑1 (MQ=255)
                                                                                                                      aCAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTa                            <  4:344449/120‑1 (MQ=255)
                                                                                                                      aCAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTa                            >  3:344449/1‑120 (MQ=255)
                                                                                                                       cAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTg                                           <  2:953669/104‑1 (MQ=255)
                                                                                                                       cAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTg                                           >  1:953669/1‑104 (MQ=255)
                                                                                                                               tAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAgg                                >  3:840502/1‑107 (MQ=255)
                                                                                                                               tAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAgg                                <  4:840502/107‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
TGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATG  >  NC_000913/4404202‑4404465

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