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

ATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAA  >  NC_000913/1209469‑1209743                                                                                                                                      |                                                                                                                                              atttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccc                                                                                                                                      <  1:306260‑M2/143‑10 (MQ=255)             gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGttt                                                                                                                          >  1:205858/1‑143 (MQ=255)                             tcccggGTATTCCTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                          >  1:109199‑M2/3‑106 (MQ=255)                               ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                        >  1:209093/1‑143 (MQ=255)                                                     aGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttggtatcgc                                                                                  <  1:280906‑M2/143‑62 (MQ=255)                                                           aCCATTGCATTTCAGCCAGCCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg                                                                            >  1:214254‑M2/1‑76 (MQ=255)                                                            cccatGCACTTCGTCCGGCCTGTTGGTGGATTTGCTGGAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                           <  2:301708/139‑1 (MQ=255)                                                             cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGACGCGATAACAAGCAAATGGCATCATTCCTGCTCTTaac                                                                          >  1:302318/1‑141 (MQ=255)                                                                         cgtagcaTGTTGGCGGAGTGGCTGACGGCCACAGTACCGCGTCTGCACCAGGAAATGCAAAGCATTTTTCCAAACCAACGTTTATGAAAATGAAGAAATAACTAGTATATTGTATCATTCCTGCTTTTACCAGGGGGATTTaa                                                              <  2:81134/140‑1 (MQ=255)                                                                                tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                                       >  1:51446‑M2/1‑55 (MQ=255)                                                                                                         ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCCTTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTACGCGTATccacaa                              >  1:297237/1‑143 (MQ=255)                                                                                                            aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcggatttgatgtggccaaccggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtc                           <  2:173730‑M2/143‑117 (MQ=255)                                                                                                              cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctgaaccgcagcgcgcgaagggcggtaagtaactggtcgt                         <  1:275363‑M2/143‑119 (MQ=255)                                                                                                                    acCAACAGGTAATGCAGAGCCTTCTCTCAAACCTTCGTTTATGAAAATGAAGAAATAACACGCAAATGGCATCATTCCTGCATCTACCCACTCGATTTAACCTGCTCATCGGCTATGCACCCGCAACAACAAATAACCAGAac                   >  2:156794/1‑143 (MQ=255)                                                                                                                    acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTCACGTGCTTACAGACTATGTACGCGTATCAACAAATGACCAGacc                   >  2:65857/1‑141 (MQ=255)                                                                                                                      cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc                 <  1:202698‑M2/143‑127 (MQ=255)                                                                                                                          aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt             >  2:231071/1‑143 (MQ=255)                                                                                                                                   agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTACCAACAAATGACCAGAACACAGATCACCAACGt    >  2:140041/1‑143 (MQ=255)                                                                                                                                     agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa  <  2:244662/143‑1 (MQ=255)                                                                                                                                      |                                                                                                                                              ATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAA  >  NC_000913/1209469‑1209743

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.