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

TACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGG  >  NC_000913/1209381‑1209867                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         ataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                                                                                                                                                                                              >  1:7885‑M2/239‑251 (MQ=255)         aaTACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCATCCACGAAGAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGAAACCAACAGGTAATTCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                                                                                                                                                                      <  2:95581/251‑1 (MQ=255)                   aTAGAGCGACCTGTATCAATCCCACGCCCGTCAGCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCGCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                                                                                                                            <  2:107272/251‑1 (MQ=255)                              aacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                                                                                                                                 >  1:146218‑M2/210‑251 (MQ=255)                                 aTCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATTTCATCATTCCt                                                                                                                                                                                                               <  1:112044/250‑1 (MQ=255)                                            gcccgTCATCCCAGCAACGAATACACTCACCACGTACATAAGGCAATTTTTTTGTCGGATAAGCATTTGCCAGTTCCGGGTATTCTTCAGCAGAAAACGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCCACAGGTAATGCAGAGCCTTCTCCCCAACCAACGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                                                                                                                                   >  2:252016/1‑251 (MQ=255)                                                   aacccAGCCTCGAATAAACTCACCACGTAAAGCATTCAAGTTATTAGTCAGATAAACCTTTGTCGGTGCCGGGGATTCTTCAGCAGAAAAAGCGGCACCATTGCATGTCAGCGAGCCTGTTGGCGGAGTGGCTGAAGGCCGCTGAACCGGTTCACCAACAGGTAATTCAGCGCCTTCTCTCAAAGCAGCGGCTAGCGAAATGGAGAAATAACGATCGAATTGCATCATTCCTGCTTTTTGCAGGGGGAttt                                                                                                                                                                                            <  2:58548/249‑1 (MQ=255)                                                     taaacgcaaccgggaatacaaatgtctttccatccccccccgaggcagatgcttttgcagcattcaacccccactgaatgataactgattgctcaataacaccaatcataaccggtatagctatgacaccggtgcgtcctaatgcagccgttgctgcattgagttttgctccttctcccaaaccaaAGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAGGGGGATTTa                                                                                                                                                                                           >  2:201558‑M2/187‑250 (MQ=255)                                                       cAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGAGGCTGAAGGCCACGGAACCGGGCCACCCACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGTCATCATTCCTGCTTTTAACAGGGGGATTTAAc                                                                                                                                                                                         >  2:361202/1‑250 (MQ=255)                                                        aGCCACGAATAAACTCACCACGTAAATCAGGCAATTTCTTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTACCAt                                                                                                                                                                                       >  1:284078/1‑251 (MQ=255)                                                             cGAATAAACTCACCACGTAAATCAGACAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTACTCTTAAGCAGAACACGCGCCACCATTGCATTTCACACAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCC‑AAACCAACGTTTATGAAACTGACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                                                                                                                                 >  1:203880/1‑251 (MQ=255)                                                                    aCTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAGTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                                                                                                                                           <  1:255917/251‑1 (MQ=255)                                                                                       ccccccataggcagatgcttttgcggcattccacccccactgactgattactgattgctccataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgcgttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAACTCACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                                                                                                                        >  2:25027‑M2/153‑251 (MQ=255)                                                                                         atttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTCATGACAATGAAGAAATAACAAGCAAATCGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg                                                                                                                                                       >  1:123186/1‑250 (MQ=255)                                                                                           ccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                                                                                                                                    <  1:12479‑M2/103‑1 (MQ=255)                                                                                                    ggATACGCATTTGCCAATTCAGGGAATTCTTCCGCAGAAAAAGCGGCACAATTGCATTTCAGCCAGCCTGTTGGCGAAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAAAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                                                                                                                             >  2:346814/1‑249 (MQ=255)                                                                                                    gctgcttttgccgcattccacccccagtgaatgattactgattgctcaatacccccaatcataaccggtatagctatgaaaccggtgagtcctaatgaagccgttgatgcattcagttttgctccttctcccaaaccagCGTTTATGAAAATGACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGTGGATTTAACATGCTTATTGGCTCTGTACGCGTATCAACAAATGACCAGAacac                                                                                                                                             >  2:282229‑M2/140‑249 (MQ=255)                                                                                                      aTAAGCCTTTGCCAGTTCCTGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                                                                                                                         <  1:419161/251‑1 (MQ=255)                                                                                                      aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGTCGGAGTGGCTGAAGGCACCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATATACGCGTATCAACAAATTCCCAGAACACAGAt                                                                                                                                         <  2:58792/251‑1 (MQ=255)                                                                                                                          ggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa                                                                                                                     <  2:21753/251‑1 (MQ=255)                                                                                                                                 aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc                                                                                                               >  1:185108‑M2/111‑250 (MQ=255)                                                                                                                                 aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTGATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACTTAATTCGCTGAACTGTGCa                                                                                                              <  2:362128‑M2/141‑1 (MQ=255)                                                                                                                                   cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACCGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAgg                                                                                                            >  2:63440/1‑251 (MQ=255)                                                                                                                                       agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACGCAGATCTAGAACTTAATGCTCTGAACTGTGCAGGATGc                                                                                                        <  1:346814/251‑1 (MQ=255)                                                                                                                                                         gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAAATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGAc                                                                                      >  1:352359/1‑251 (MQ=255)                                                                                                                                                                                        gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCgg                                                       >  1:98209‑M2/56‑251 (MQ=255)                                                                                                                                                                                                    aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAAc                                           >  2:150867/1‑251 (MQ=255)                                                                                                                                                                                                       gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGc                                        >  2:279492‑M2/41‑251 (MQ=255)                                                                                                                                                                                                                   ggcaaTGCAGAGCCTTCTCCCAACCCAACGGTTATGAAAATGAAGATATAACAAGCAAATTGCATCATTCCTGCTTTTACCATGGGGATTTAACATGCTGATTGGCGATGTACGCGGATCAACAAATGACTAGAATACAGATCTACAACGTAATGGGCTGAACTGTGCAGGATGCTAGCGTATTTTTGAATACAAGATAATCGGCATAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCATGACATTAt                            <  1:52458/248‑1 (MQ=255)                                                                                                                                                                                                                                           ccAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAGGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCGCAGGACAGTATCGGCAGGTGACACTCTGGTTGTCt    <  2:185108/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaCGTTTATAAAAATGAAGAAATAACAAGCAAATGGAATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGGCACTCTGGTTGTCTgg  >  2:256802/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         TACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGG  >  NC_000913/1209381‑1209867

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