New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	18 (0.310)	19 (0.350)	17/450	1.2	57.8%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	11 (0.210)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAA  >  NC_000913/1209370‑1209869                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           ccactgcatggataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccgttgcgtcctaaggtagccgttgctgcattgagttttgctccttcgcccaaaccaaCg                                                                                                                                                                                                                                                           <  2:472904‑M2/2‑1 (MQ=255)    ccGTTGCATACCCCTGAATACTTAAAAGAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACTCCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                                                                                                                                        <  1:257336/251‑1 (MQ=255)     ggcagtgagaggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatttagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                                                                                                                                       <  2:248486‑M2/6‑1 (MQ=255)                            aaaTAGAGCGACCTGTATCAATCCCCCGACAGTCATACCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTGGTCGGACAAGCCTTTGCCAGTTCCGGGTATGCTTCAGAAGAAAAAGCGGCACCATGTCATTGCAGCCATCCTGTTGGCGGAGTGGCTGAAGGCCACCGTACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGaataaata                                                                                                                                                                                                                                       <  1:576986/244‑1 (MQ=255)                               cataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                                                                                                                                             <  1:541846‑M2/32‑1 (MQ=255)                                         tGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                                                                                                                                   >  1:32675/1‑251 (MQ=255)                                                     ccgcccgTCATCCCATCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTTTCGGATTAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                                                                                                                                       <  1:65294/251‑1 (MQ=255)                                                              aTCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCATCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                                                                                                                              <  2:472764/251‑1 (MQ=255)                                                                                       gtatttcaatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcctaaccggtatagatatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                                                                                                                                     >  1:213875‑M2/164‑251 (MQ=255)                                                                                                        atttGTCGGGTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCGCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGTCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAGTGGCTATGTACGCGTATCATCAAATGACCAg                                                                                                                                                    <  2:248134/251‑1 (MQ=255)                                                                                                          aggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAGTAACaa                                                                                                                                                                                                                                   <  1:164505‑M2/26‑1 (MQ=255)                                                                                                          aggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                                                                                                                                   >  2:164505‑M2/145‑170 (MQ=255)                                                                                                           tGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                                                 >  2:413438/1‑251 (MQ=255)                                                                                                                                    gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg                                                                                                                         >  2:506446/1‑250 (MQ=255)                                                                                                                                    ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGACAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg                                                                                                                         >  2:239949‑M2/119‑250 (MQ=255)                                                                                                                                     cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa                                                                                                                        >  1:365078‑M2/118‑250 (MQ=255)                                                                                                                                        actgaatgattactgattgctcaataccaccaatcataaccggtatagctatgcaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGCACACAGATCTACAACGTAATGCGCTGAACtg                                                                                                                    >  1:451999‑M2/115‑251 (MQ=255)                                                                                                                                        actgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtg                                                                                                                    >  1:245139‑M2/115‑251 (MQ=255)                                                                                                                                                         tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTg                                                                                                    >  1:8255‑M2/98‑250 (MQ=255)                                                                                                                                                                     cATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAACTGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGAc                                                                                        >  2:494586/1‑250 (MQ=255)                                                                                                                                                                            gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAg                                                                                <  2:32675/251‑1 (MQ=255)                                                                                                                                                                                  gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCggccc                                                                          >  2:283321‑M2/73‑249 (MQ=255)                                                                                                                                                                                      agctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACaa                                                                        >  1:384642‑M2/69‑249 (MQ=255)                                                                                                                                                                                         tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTcc                                                                   >  1:217031‑M2/66‑251 (MQ=255)                                                                                                                                                                                                            aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCt                                                                                                     >  2:60796‑M2/47‑198 (MQ=255)                                                                                                                                                                                                            aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCt                                                                                                     <  1:60796‑M2/152‑1 (MQ=255)                                                                                                                                                                                                               gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGAtttt                                                                                               <  2:340017‑M2/158‑1 (MQ=255)                                                                                                                                                                                                               gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGAtttt                                                                                               >  1:340017‑M2/44‑201 (MQ=255)                                                                                                                                                                                                               gtagccgttgcagcattgagttttgctccttctcccacaacaaCGTTTATGAACATGAAGAAATAACAAGCAACTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTAATTGGCTATGTACGCGAATCAACAAAAGACCAGAACAAAGATCTACAACGTAATGCGCTGAACTGTGAAGGATGCGCGCTGATTGTAGAAGACAAGATAAGCGGCACAAAGTCAGAAAGGCAGGGACTGaaaaaa                                              >  2:255766‑M2/44‑250 (MQ=255)                                                                                                                                                                                                                               gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATc                             >  1:29497/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGGGACACTCTGGtt          >  2:456852/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aaCGTTTATGAAAATGAGGAACTAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTAATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGAGCTGAACTGTGCAGGATGCGGGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGAACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGCAACTCTGGTTGTCTgg    >  2:268075/1‑251 (MQ=255)                                                                                                                                                                                                                                                          cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAGCAACAAATGCCTAGAACACAGATCGACAACGTAATGCTCTGAACTGTGCAGGATGCGAGCTGATTTTTAAAGATAAGATAAGCGGCACAAAGTCCGAAAGCCCTGGACTGAAAAAACGGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGaa  <  2:351513/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           CCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAA  >  NC_000913/1209370‑1209869

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