New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	4 (0.100)	19 (0.650)	12/98	0.7	88.3%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	2 (0.070)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209553‑1209681                                                                   |                                                               ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                               <  1:40235‑M2/2‑1 (MQ=255)  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                              >  1:1304056‑M2/66‑68 (MQ=255)  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                              >  1:829014‑M2/66‑68 (MQ=255)  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                              >  1:636617‑M2/66‑68 (MQ=255)  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagtttagctccttctcccacaccaaCGt                                                              >  1:1470770‑M2/66‑68 (MQ=255)  catgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaacgc                                                              <  2:1152083‑M2/3‑2 (MQ=255)                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                         >  2:48418‑M2/45‑68 (MQ=255)                       tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                         >  2:1093492‑M2/45‑68 (MQ=255)                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                        >  1:1214235‑M2/44‑68 (MQ=255)                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAAAACa                                        >  2:950221‑M2/44‑68 (MQ=255)                        gtagccgttgctgcattgagtattgctccttctcccaaaccgaCGTTTATGAAAATGAAGAAATAACa                                        >  2:935483‑M2/44‑68 (MQ=255)                          agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                      >  2:288472‑M2/42‑68 (MQ=255)                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                    >  2:696217‑M2/40‑68 (MQ=255)                             gACACCAACAGGTAATGCAGAGCCATCTCACAAACAAACGTTTATGAAAATGAAGAAATAACAAGCaa                                   >  2:1021981/1‑68 (MQ=255)                              gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                  >  1:766298‑M2/38‑68 (MQ=255)                                 gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                               <  1:243424‑M2/34‑1 (MQ=255)                                        tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        <  2:766294‑M2/41‑1 (MQ=255)                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:167204‑M2/46‑1 (MQ=255)                                             gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   >  2:736987/1‑68 (MQ=255)                                                gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                >  1:127689‑M2/20‑68 (MQ=255)                                                  tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa              <  2:1214231‑M2/51‑1 (MQ=255)                                                          ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg      <  1:1384964‑M2/59‑1 (MQ=255)                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:109802‑M2/60‑1 (MQ=255)                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:707996‑M2/60‑1 (MQ=255)                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:48419‑M2/60‑1 (MQ=255)                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAGCATTCCTGCTTTTACCAGGGGGa     <  2:388035‑M2/60‑1 (MQ=255)                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  <  1:736991‑M2/63‑1 (MQ=255)                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  >  2:870458‑M2/6‑68 (MQ=255)                                                                   |                                                               GCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209553‑1209681

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