New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120779316 (0.170)4 (0.040)
+G		 4/502 9.8 28.6% coding (293/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209613 4 (0.040)pseudogene (43/537 nt) stfE e14 prophage; putative side tail fiber protein fragment
Rejected: Coverage evenness skew score above cutoff.

TGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATT  >  NC_000913/1209387‑1209839
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
tGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                                                                                                                                          <  1:96974/285‑1 (MQ=255)
                                                                       tAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTa                                                                                                   >  1:295217/1‑285 (MQ=255)
                                                                                                        tGCCAGTTCCGGGTATTTTTCAGCAGAAAAAGCGGCACCATTGCATTTCACCCAGCCTGTTTGCGGAGTGGCCGATGGCCACGGAACCGGGACACCAACAGGCAATGCAGTGCCTTCTCCTATACCAACGTTTATGATCTTGAAGAAATAACCAGCAAATGGCATCATTCCTGCTTTTACCAGGGCGATTTAACATGCTTATTGGTCATGTACTCGTATTAACAAATTACCTGAACCCAAATCTACCACGTAATGCGCTGAACGGTGCAGGATGCGCTCTGCttt                                                                  >  2:253122/1‑285 (MQ=255)
                                                                                                                                                                         gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAAAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCCCAAAGTCCGAAAGGCCGGGACTGAAAAAAcctttctctattac  >  1:706804/1‑271 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
TGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATT  >  NC_000913/1209387‑1209839

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