New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	5 (0.130)	30 (0.860)	26/252	0.3	87.7%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	= 1209602	4 (0.110)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACG  >  NC_000913/1209475‑1209740                                                                                                                                |                                                                                                                                           ttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                                                                                             >  1:620535‑M2/1‑128 (MQ=255)   gTCGGCTAAGCCCTCCGCAGTTCCCGGTATGCTTCAGCCGAAAGAGCGGCACCGTTCCATTTCAGCCAGCCTGTTCGCGGAGTGGCTGAAGGCCACCGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                           >  2:438503‑M2/1‑126 (MQ=255)     cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                                                                                         <  1:225714‑M2/143‑20 (MQ=255)       gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttt                                                                                                                       <  1:269722‑M2/143‑22 (MQ=255)         tAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                     >  2:225714‑M2/1‑120 (MQ=255)            gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgag                                                                                                                  >  2:176159‑M2/1‑117 (MQ=255)             cgcttcccCGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                 >  2:607232‑M2/4‑116 (MQ=255)             ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                 <  1:586847‑M2/143‑28 (MQ=255)                 tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                                                             <  1:307373/143‑1 (MQ=255)                         ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaca                                                                                                     >  1:96381‑M2/1‑104 (MQ=255)                           gggTATTCGTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                                                                                                   >  1:134774‑M2/1‑102 (MQ=255)                           gggTATCCGTCAGCAGAAAAAGCGGCACCATTGCAATTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                                                   <  1:48726/143‑1 (MQ=255)                                  cctcACCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCATACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                            <  1:297996/141‑1 (MQ=255)                                       gcagAACAAGCGGCACCATTGCACGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                       >  1:182160‑M2/1‑90 (MQ=255)                                       gcagAAAAGGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                       >  2:93104‑M2/1‑90 (MQ=255)                                             aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                 >  1:703961/1‑143 (MQ=255)                                             aaaGCGGCACCATTGCATGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc                                                                                 >  1:405819‑M2/1‑84 (MQ=255)                                                  ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaccgtgccatccgatttcatatcgccac                                                                            <  2:443950‑M2/143‑65 (MQ=255)                                                       cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGATCCGGGACACCAACAGGTAATGCATAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgga                                                                       >  1:290645‑M2/1‑74 (MQ=255)                                                             aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggactcttg                                                                 <  2:491580‑M2/143‑76 (MQ=255)                                                                tCGGCCGGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                              >  2:308716‑M2/1‑65 (MQ=255)                                                                gcgggcagcgtggtgGCGGGGTGGCTGGAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccgtcgcccaagccaaggtgtgcgagagccggttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                              >  2:190863‑M2/6‑65 (MQ=255)                                                                      agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaac                                                        >  1:572389‑M2/1‑59 (MQ=255)                                                                       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                       >  2:231067‑M2/1‑58 (MQ=255)                                                                        ccgattgccgGAGTGTCTCAAGGCCACGGTACCGGTACTCCAACAGGTAATGCAGAGccttctcccaaaccgaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacag                                                      >  2:108507‑M2/5‑57 (MQ=255)                                                                            ttGGCGGAGGGGCTGAAGGCCTCGGGACCGCGACACCAAAAGGTAGTGGATAGcctgctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagc                                                  >  1:446074‑M2/1‑53 (MQ=255)                                                                                       gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccacacggattcttgcgccctaacagcagcgcgcgcagtgc                                       <  1:569542‑M2/143‑102 (MQ=255)                                                                                             ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgacttgatatcgccaaacggattcgggcggcttaacagcagcgcgcggagtgcggtaag                                 <  2:212078‑M2/143‑108 (MQ=255)                                                                                                         gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtg                     >  2:269722‑M2/1‑24 (MQ=255)                                                                                                             cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattcgataccgccaaacggattcttgaggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgc                 >  1:544169‑M2/1‑20 (MQ=255)                                                                                                                     ggtgagtcagagcctcctgccaaaccaaggttttctggaggcgcttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctgaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctc         >  2:622854‑M2/8‑12 (MQ=255)                                                                                                                     ggTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttcacgccaccgcgcgaagttcggtaagtaactggtcgtgccgccccttctc         >  1:491580‑M2/1‑12 (MQ=255)                                                                                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttatcc        <  1:581456‑M2/143‑133 (MQ=255)                                                                                                                        aaTGCAGAGccttctcccaaaccagggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccag      >  1:472824‑M2/1‑9 (MQ=255)                                                                                                                            cAGAGccttctcccaagccaaggttttcgagagccgtttgcacggtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg  >  2:461912‑M2/1‑5 (MQ=255)                                                                                                                                |                                                                                                                                           TTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACG  >  NC_000913/1209475‑1209740

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