New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	1207790 =	5 (0.160)	19 (0.680)	19/256	0.5	75.3%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_bme_pgi	1209619 =	8 (0.290)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_bme_pgi/1209479‑1209764                                                                                                                                             |                                                                                                                                                   cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                                                                                            <  1:61124/149‑1 (MQ=255) agatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgccttgaggtttgctccttctcccaaaccaaCGTTTATGa                                                                                                                                            >  2:297897‑M2/141‑149 (MQ=255)   atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                          <  2:438181‑M2/11‑1 (MQ=255)        tttgcggcattccccccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatctagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                                     <  1:155594‑M2/16‑1 (MQ=255)                   ccacccccactgaatcgttactgagtgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                                          <  1:476769‑M2/27‑1 (MQ=255)                           actgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGc                                                                                                                 >  1:327215‑M2/115‑149 (MQ=255)                               ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGcatca                                                                                                             >  2:222585/1‑149 (MQ=255)                                  agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAGGGCCACGGAACCGGGAAACCAACAGGTAATGCCGAGCCTTCTCCCAAACCAACGTTTCTGAAAATGAAGAAATAACAAG‑CAAACGGCGTCATTc                                                                                                          >  1:159175/1‑149 (MQ=255)                                     agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTg                                                                                                       <  1:354696/149‑1 (MQ=255)                                           ttgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTa                                                                                                 >  1:517706‑M2/99‑149 (MQ=255)                                             actcaagagcaacactcacacccggtacagtgatgaaaccggtgcctcctaatgtagccgttgctgcattgagttttgcgcctgctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAC‑CAAATGGCATCATTCCTGCTTTTAcc                                                                                               <  2:444169‑M2/53‑1 (MQ=255)                                                caataccaccaatcataaccggtatagctatgaaaccgttgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAgg                                                                                            <  1:437651‑M2/56‑1 (MQ=255)                                                 aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCttt                                                                                                   >  2:399564‑M2/93‑141 (MQ=255)                                                 aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCttt                                                                                                   <  1:399564‑M2/49‑1 (MQ=255)                                                       accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAT‑CAGATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                     >  1:238649‑M2/87‑149 (MQ=255)                                                           atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                 >  2:382248‑M2/83‑149 (MQ=255)                                                              ataaccggtatagctatgcaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                                              <  1:141923‑M2/70‑1 (MQ=255)                                                                      tgtgGGCGGAGGAGCGGAAGGACACGGAACCGCGAAACCAACAGGTAAGGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGCAATAACAAC‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                      <  2:164186/149‑1 (MQ=255)                                                                              gAGTGGCTTAAGGCCACGGAACCGGGACACCAACAGGTACTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGACGAAATAACAAG‑CAAATTGCATCATTCCTGCTTTTACCAGTGCGATTTAACATTCGTATTGGCTATGTAct                                                              >  2:436048/1‑148 (MQ=255)                                                                                    cggggcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                        <  1:149003‑M2/92‑1 (MQ=255)                                                                                          gccctcatgtggccgttgcgtcattgcgttttgctccgtctcccaacctaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATGCTCGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                  <  1:242898‑M2/98‑1 (MQ=255)                                                                                                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                >  2:280273‑M2/34‑149 (MQ=255)                                                                                                                 ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAAGCGGATCTACAAc                           >  2:175493‑M2/29‑149 (MQ=255)                                                                                                                 ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                           <  1:222585/149‑1 (MQ=255)                                                                                                                     gttgtgctcctgcttccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCTTCATTCCTGCTTTTACCATGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                       <  2:221024‑M2/125‑1 (MQ=255)                                                                                                                          gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACCCGTATCAACAAATGACCAGAACACAGATCTACAACGTAACgcgc                  >  2:482442/1‑149 (MQ=255)                                                                                                                           ctccttcccacaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCAACATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                 <  2:126375‑M2/131‑1 (MQ=255)                                                                                                                               ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc             >  2:359592‑M2/15‑149 (MQ=255)                                                                                                                                  tcccaaaccaaCGTTTATGAAAATGAAGAAAT‑ACACGACAAAGGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt          >  1:540338‑M2/12‑149 (MQ=255)                                                                                                                                      aaaccaaCGTTTATGAAAATGAAGAAATAACACG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg      <  1:232769‑M2/142‑1 (MQ=255)                                                                                                                                          caaCGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCGGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATg  >  1:411718‑M2/4‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                   CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAG‑CAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_bme_pgi/1209479‑1209764

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