New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078054 (0.520)5 (0.750) 5/224 0.2 54.1% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 5 (0.750)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAG  >  NC_000913/1209524‑1209723
                                                                              |                                                                                                                         
ggCACC‑TTGGATTTCAGGCAAGCCGTGGGCGGGGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGGGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                        <  1:484038/129‑1 (MQ=255)
 gCACCATTGCATTGCGGCCATGCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTTTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                        <  1:316881/129‑1 (MQ=255)
                       cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattc                                                   >  1:565164‑M2/1‑56 (MQ=255)
                        tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                  >  1:633957‑M2/1‑55 (MQ=255)
                           tGGCGGGGGGGGTGAAGGCCACGGGACCGGGAAACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg                                              >  1:550045‑M2/1‑52 (MQ=255)
                                    ggcgaaaGGCCCCGGAACC‑GGACCCCAACAGGGAAGGCAAAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                    <  1:85552/124‑1 (MQ=255)
                                     gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                                    >  1:574305/1‑129 (MQ=255)
                                       tGAAGGCCCCCGAAACGGGAACCAAAAAGGTAAGTAGGAGCTTTTCCCCAAACC‑ACGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                 <  1:160318/129‑1 (MQ=255)
                                              cACGGAACCTGGACACCCACAGGGAATTCAGAGccttcccccaaaccaaggttttcgagggccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                           >  1:604676‑M2/1‑33 (MQ=255)
                                                                       aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcaaagtccgggacttaactggtctgccc  <  1:445112‑M2/129‑122 (MQ=255)
                                                                              |                                                                                                                         
GGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAG  >  NC_000913/1209524‑1209723

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