New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	1207790 =	2 (0.140)	16 (1.270)	15/228	0.1	94.8%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_bme_pgi	1209619 =	0 (0.000)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGC  >  NC_000913_3_bme_pgi/1209470‑1209746                                                                                                                                                      |                                                                                                                                cccagacgcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                 <  1:90223‑M2/1‑1 (MQ=255)         cagatgcttttgctgcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctattaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccacCGTTTATGa                                                                                                                         <  2:91127‑M2/9‑1 (MQ=255)           gatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                       >  1:79289‑M2/140‑150 (MQ=255)            atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcatagccggtatagctatgaaaccggttcttcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTGTATGaaaa                                                                                                                      <  1:130851‑M2/12‑1 (MQ=255)            atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                                                      <  1:63992‑M2/12‑1 (MQ=255)                                             ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                                              >  4:29598‑M2/106‑109 (MQ=255)                                             ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                                              <  3:29598‑M2/4‑1 (MQ=255)                                                        tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                          >  3:32909‑M2/95‑150 (MQ=255)                                                         caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                         <  2:79289‑M2/57‑1 (MQ=255)                                                                accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                              >  4:79895‑M2/87‑138 (MQ=255)                                                                accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                              <  3:79895‑M2/52‑1 (MQ=255)                                                                      cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctcctactcccaacccaaTGTTTATGATAATGAAGAAATAACAAGTGAATGGTATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                            <  1:60778‑M2/70‑1 (MQ=255)                                                                         aaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                         <  4:32909‑M2/73‑1 (MQ=255)                                                                            cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                      >  3:46462‑M2/75‑150 (MQ=255)                                                                                                                     ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTa             >  2:17005‑M2/34‑150 (MQ=255)                                                                                                                                tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgc  <  1:17005‑M2/128‑1 (MQ=255)                                                                                                                                                     aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                   <  2:66791‑M2/79‑1 (MQ=255)                                                                                                                                                     aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                   >  1:66791‑M2/2‑80 (MQ=255)                                                                                                                                                      |                                                                                                                                TTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGC  >  NC_000913_3_bme_pgi/1209470‑1209746

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