New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	21 (1.160)	6 (0.380)	6/232	NT	23.8%	coding (290/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	1209619 =	20 (1.260)					pseudogene (37/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAC  >  NC_000913/1209481‑1209753                                                                                                                                           |                                                                                                                                       gATAAGCCTTTGCCAGTTCCTGGTATTCTTCAGCAGAAAAAGCGGCACTATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                    <  2:180955/143‑1 (MQ=255)            gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                         <  2:238098/143‑1 (MQ=255)            gCCAGCTCCGGGTATTGTGCAGCAGAAAGAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTAGCTGAAGGCCACGGAACGGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                         <  2:23784/143‑1 (MQ=255)                         aTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGGACTGACACACCAAAAGAAAAAGGAGAAACAACTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                            >  1:345006/1‑143 (MQ=255)                            cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                         <  2:151332/143‑1 (MQ=255)                             ttcAGCAGATAAATCGTTACCATTTCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                        <  1:112300/143‑1 (MQ=255)                             ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                        >  2:285460/1‑143 (MQ=255)                                     aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGCTTATGAAAAAGAAGAAATAACAAGCAAATAGcaacata                                                                                                >  2:278381/1‑142 (MQ=255)                                      tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAACTAACAAGCAAATGGCATCATTc                                                                                               >  1:203801‑M2/102‑143 (MQ=255)                                            ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGCAATCACACACAACGGGCATCACTCCTGCtt                                                                                         >  1:47356‑M2/96‑143 (MQ=255)                                             tcaataccaccaataataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgcgttttgctccttctgccaaaccaaCGTTTATGAAAATGAGGAACTAAGAGGCAACTGGCAACACTACTGCttt                                                                                        >  2:48644‑M2/95‑143 (MQ=255)                                             gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGAATGTAAAATTAAGGAAAAAACAAGCACATGGCATAATTCCTgatgc                                                                                        >  2:194258/1‑139 (MQ=255)                                             gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                        >  2:266924/1‑143 (MQ=255)                                                ccATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGCAGGCCGCGGGACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCATCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                     <  2:313610/143‑1 (MQ=255)                                                     tcATTTCAGCCAGCGTGTTGGCGGAGTGGATGAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                <  1:320803/142‑1 (MQ=255)                                                      cATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGCTAATGCAGAGCCTTCTACCAAACAAAAAATAATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggga                                                                               >  1:46427/1‑142 (MQ=255)                                                        tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAACAGAGGGa                                                                             >  1:123964/1‑143 (MQ=255)                                                            ataactgttatagctatgaaaccggcgcgtcctaatggagccggtgctgcagtgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                                                                         <  1:57363‑M2/64‑1 (MQ=255)                                                              aaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                       <  1:142724‑M2/66‑1 (MQ=255)                                                                         gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                            >  2:126122/1‑143 (MQ=255)                                                                                    gAAGGCCACGGAACCGCGACACCAGCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTCTCCTTGTGCGCTTTAACATGCTTATTGGCTATGTAcgcg                                                 <  2:176705/143‑1 (MQ=255)                                                                                       ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCATTGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                              <  1:285460/143‑1 (MQ=255)                                                                                                 tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATATAACATGCTTATTGACTATCTACGCGTATCAACAAATGa                                    >  2:38691‑M2/43‑143 (MQ=255)                                                                                                          cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCGCTGCACCCGTATCAACAAATGACCAGAaccc                           >  1:64890/1‑141 (MQ=255)                                                                                                          cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                           >  1:138143/1‑143 (MQ=255)                                                                                                          cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATACCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCAATGTACGCGTATCAACAAATGACCAGAacaa                           >  2:299946/1‑142 (MQ=255)                                                                                                              tggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                       <  1:148335/142‑1 (MQ=255)                                                                                                                                 cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa    >  2:129619‑M2/11‑143 (MQ=255)                                                                                                                                  ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                             <  2:185690‑M2/44‑1 (MQ=255)                                                                                                                                  ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                             >  1:185690‑M2/10‑53 (MQ=255)                                                                                                                                   caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc  >  2:277052‑M2/9‑143 (MQ=255)                                                                                                                                           |                                                                                                                                       GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAC  >  NC_000913/1209481‑1209753

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