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.

GCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATGGAAGAAGATGATCTCTGCGATGACCAATTTGCCCGATAATATT  >  NC_000913/4404270‑4404508
                                                                                           |                                                                                                                                                   
gCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCt                                                                                             >  4:523123/1‑148 (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)
                                                                 cACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGaa                                                                             <  2:1117331/99‑1 (MQ=255)
                                                                 cACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGaa                                                                             >  1:1117331/1‑99 (MQ=255)
                                                                  aCCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGt                                                                          >  2:987736/1‑101 (MQ=255)
                                                                  aCCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGt                                                                          <  1:987736/101‑1 (MQ=255)
                                                                  aCCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATgg                                            <  4:987349/131‑1 (MQ=255)
                                                                  aCCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATgg                                            >  3:987349/1‑131 (MQ=255)
                                                                    cGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCg                                                                                                                    <  2:393472/57‑1 (MQ=255)
                                                                    cGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCg                                                                                                                    >  1:393472/1‑57 (MQ=255)
                                                                            cAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATGGAAGAAGATGATCTCTGCGATGACCAAt                 >  3:841016/1‑148 (MQ=255)
                                                                                  tGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATGGAAGAAGATGATCTCTGCGATGACCAATTTGccc           >  1:813913/1‑148 (MQ=255)
                                                                                           ttttctttttctATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATGGAAGAAGATGATCTCTGCGATGACCAATTTGCCCGATAATAtt  <  4:841016/148‑1 (MQ=255)
                                                                                           |                                                                                                                                                   
GCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTGGAAGGTTTAGGGCCGATGCTTTACCCGAAGGCATGGAAGAAGATGATCTCTGCGATGACCAATTTGCCCGATAATATT  >  NC_000913/4404270‑4404508

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