New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078053 (0.200)7 (0.540) 7/232 NT 75.1% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 2 (0.150)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGT  >  NC_000913/1209463‑1209741
                                                                                                                                            |                                                                                                                                           
cAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGTATTTCAGCCAGCCTGTTGGCGGAGTGGCTGCAGGCCGCGGAACCGGGCCTCCAACAGGTTATGCAGAGcct                                                                                                                                          <  1:98102‑M2/143‑4 (MQ=255)
               tCGGATAAGCCCTTGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                                           <  1:134601/143‑1 (MQ=255)
                             gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATACAGAGgcttctcccaaaccaaggttttcgagaaccgg                                                                                                             <  2:144644‑M2/143‑33 (MQ=255)
                              ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCCACAGCCATTTGCACCAACGGGTACTGCAGAGccttctcccaaaccaaggttttcgcgagccgtt                                                                                                            <  1:220730‑M2/143‑34 (MQ=255)
                                         tATTCGTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCGGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcc                                                                                                 >  1:36057‑M2/1‑99 (MQ=255)
                                            tgctcagcAGAAAAGGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                              >  1:168866‑M2/4‑96 (MQ=255)
                                                  gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgcttt                                                                                        <  2:153580‑M2/143‑54 (MQ=255)
                                                                     tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCATAGccttctcccaaaccaaggttttcgagagcctttttcacagtgccatccgatttgatatcgccaaacggattc                                                                     >  2:265716‑M2/1‑71 (MQ=255)
                                                                               ccagcgTGGTGGCCGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                           <  2:254623/143‑1 (MQ=255)
                                                                                                tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagt                                          <  2:153905‑M2/143‑100 (MQ=255)
                                                                                                          gacccGGAACCTGTACACCAACAGGTAATGCAGAGCCTTCTCGCCAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                 <  2:127619/139‑1 (MQ=255)
                                                                                                                                         agagccttgtcccaaacttaggttttcgagagccgtttgcacagtgccattggatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg  >  2:42872‑M2/1‑4 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGT  >  NC_000913/1209463‑1209741

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