New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	11 (0.380)	10 (0.370)	9/452	1.0	56.7%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	5 (0.190)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

TTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGG  >  NC_000913/1209376‑1209867                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         cagtgataggtgcgttgacattgtcataatttgcaacaacggcaacacaggcattatgaaacgcaaccgggaagacaactgtatttccatcccccccagaggcagatgctgttgcggcattccaccctcagtgaatgattactgattgctcaatgccgccaatcatagccggtagagctatgaaactggtgcggcctaatgttgccgttggtgcgttgatttttgcgtcttctcccaaaccaaCGTTTATg                                                                                                                                                                                                                                                   <  2:16353‑M2/8‑1 (MQ=255)      ataggcgggataacattgtaataatttgcaacaccggcaacactgacattattaaacgcaaacgggaatacacctgtatttccatcccccccagaggcagatgcttttgcggcattccccccccactgaatgattactgattgctcaataccaccaatactaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCATTTATGAAAAt                                                                                                                                                                                                                                              >  2:46197‑M2/239‑251 (MQ=255)             cgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                                                                                                                                       <  1:1314‑M2/20‑1 (MQ=255)              gctgacattgtcataatttgcaacaacggcaacacaggcattagtaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                                                                                                                      <  1:265987‑M2/21‑1 (MQ=255)                               gaccgttATCCATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCTCCACGTAAATCAGTCAATTTCTTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                                                                                                                                                                     <  2:173988/245‑1 (MQ=255)                                       acggcaacacaggcattagtaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccatccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAGGAAAAGGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                                                                                                                             <  2:181876‑M2/46‑1 (MQ=255)                                            aacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaacccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                                                                                                                         >  2:198801‑M2/201‑250 (MQ=255)                                                                              actgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                                                                                                                                      >  1:197400‑M2/167‑251 (MQ=255)                                                                                      aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                                                                                                                                              <  1:140668/251‑1 (MQ=255)                                                                                       tCAGGCAATTTATTTGTCGGATAAGCCTTTGCAAGTTCCGGAAAATCTTCCGCAGAAAAAGCGGCACCATTACATTTCAGCCCGCCTGTAGGCGGAATGGCTGAAGGCCACGGAACAGGAACACCACCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTGAACATGCTTATTGGCTATGTACGCGTATca                                                                                                                                                              >  1:212257/1‑250 (MQ=255)                                                                                                                          ccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAACTGAAGAAATAACACGCAACTCGCATCATACCTGCTTTTACCAGGGGGATTTAACAGGCTTATTGTCTAAGTACGCGTATCAACAAATGACCAGACCACAGATCTACACCGTAATgcg                                                                                                                          >  2:112677‑M2/123‑251 (MQ=255)                                                                                                                                  tctgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCCTGCTTTGACCAGGGGGATTTAACATGCTGATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                                                                                                                              <  2:13000‑M2/125‑1 (MQ=255)                                                                                                                                  actgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                                                                                                                              >  1:13000‑M2/115‑239 (MQ=255)                                                                                                                                  aTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAGCTACAACGTAATGCGCTGAACt                                                                                                                   <  1:255403/250‑1 (MQ=255)                                                                                                                                                                                                                              gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAg                      <  2:273922/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACCCGGGATGTCTgg  >  2:241990/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         TTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGG  >  NC_000913/1209376‑1209867

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