New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	21 (0.480)	15 (0.390)	14/244	NT	41.4%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	= 1209602	24 (0.620)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913/1209465‑1209748                                                                                                                                          |                                                                                                                                                   ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCCCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                                                                                                         <  2:482126‑M2/149‑12 (MQ=255)   cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAGAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                                                                                       <  1:616926‑M2/149‑14 (MQ=255)   cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                                                                                       <  2:457050‑M2/149‑14 (MQ=255)   cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                                                                                       <  1:411778‑M2/149‑14 (MQ=255)   cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                                                                                       <  1:558695‑M2/149‑14 (MQ=255)               cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCCGCCTTTTGGCGGGGTGGCTGAAGGCCACGGAACCGGGACCCCAACAGGTAATCCGGAGccctctcccaaaccaaggtttcgga                                                                                                                           <  2:127292‑M2/149‑26 (MQ=255)                 gATACGCCTTTTCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAACCTTCTCCCAAACCGACGGTTATaaaa                                                                                                                         >  2:718317/1‑149 (MQ=255)                      gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                    >  1:330255/1‑149 (MQ=255)                         tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaagaa                                                                                                                 >  2:715589/1‑149 (MQ=255)                         tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGATGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaagaa                                                                                                                 >  2:627290/1‑149 (MQ=255)                                 ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                                                         <  2:436652/149‑1 (MQ=255)                                 ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                                                         <  2:90717/149‑1 (MQ=255)                                   ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcca                                                                                                       <  2:563155‑M2/149‑46 (MQ=255)                                             ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                             >  2:598105/1‑149 (MQ=255)                                                agcagAAAAAGCGGTACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                          >  2:605541/1‑149 (MQ=255)                                                  cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                        >  1:651714/1‑149 (MQ=255)                                                      aaaaGCGTCACCATTGCACTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                    <  2:411778/149‑1 (MQ=255)                                                      aaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                    <  1:457050/149‑1 (MQ=255)                                                      aaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                    <  2:558695/149‑1 (MQ=255)                                                              cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattc                                                                            <  2:4293‑M2/149‑73 (MQ=255)                                                               aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                                           <  2:222677‑M2/149‑74 (MQ=255)                                                               aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                                           <  2:293190‑M2/149‑74 (MQ=255)                                                                  agttcatgtCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                                        >  1:66377‑M2/5‑73 (MQ=255)                                                                        tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctta                                                                  <  2:494747‑M2/149‑83 (MQ=255)                                                                        tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                  <  2:651714/149‑1 (MQ=255)                                                                           cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                               >  2:707219/1‑149 (MQ=255)                                                                             gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttcacagc                                                             <  1:106141‑M2/149‑88 (MQ=255)                                                                                       tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                   <  1:496669/149‑1 (MQ=255)                                                                                         gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcacgaggt                                                 <  1:322118‑M2/149‑100 (MQ=255)                                                                                          cGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgggagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaccgcagcgcgcgaagtg                                                <  1:114751‑M2/149‑101 (MQ=255)                                                                                            gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCCTCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTGTGTAcg                                              >  1:292977/1‑149 (MQ=255)                                                                                             aGTGGCTGAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                             <  2:726050/149‑1 (MQ=255)                                                                                                 gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaag                                         <  1:590797‑M2/149‑108 (MQ=255)                                                                                                    gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaa                                      >  1:653710‑M2/1‑39 (MQ=255)                                                                                                     aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                     <  2:330255/149‑1 (MQ=255)                                                                                                     aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGAACATGCAACTCTTCCTGCTTTTACCAGGGGGATTTAACATGTTTATTGGCTGTTTACGCTTATcaac                                     >  1:375777/1‑149 (MQ=255)                                                                                                     aaGGCCACGGAACCGGAATAGTAGCAATTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGTCATCATTCCTTCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                     <  2:547735/149‑1 (MQ=255)                                                                                                     aaGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                     <  2:361863/149‑1 (MQ=255)                                                                                                       ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactg                                   >  2:616926‑M2/1‑36 (MQ=255)                                                                                                        gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcgggaagtaactgg                                  <  2:644167‑M2/149‑115 (MQ=255)                                                                                                                   gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                       <  1:184249‑M2/149‑126 (MQ=255)                                                                                                                      acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgacgtgcggtaagtaactggtcgtgccgcccctt                    <  1:233975‑M2/149‑129 (MQ=255)                                                                                                                        acCAACAGGTAATGCGAAGCCTTCTCCCAAACCAACGTATATGAAAATGACGAAATAACAAGCAAATGGCATCAGTCCCCCTTTTACAAGTGGGATTAACCTTGATTAGTGGCTATGTACGCGGATCAACAAATGACCAGATCACAGAt                  <  2:375777/149‑1 (MQ=255)                                                                                                                          cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                >  1:738213/1‑149 (MQ=255)                                                                                                                                   aTGCAGAGCCTTCTCCCAAACCAACGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCAGCTGTTACAAGGGGGATTTAACACGCTACTTGGCAATGTACGCGAAACAACAAATAACCAGAACACAGACCTACAACGTaa       >  2:132605/1‑149 (MQ=255)                                                                                                                                      ttgcgccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccg    >  2:387096‑M2/5‑5 (MQ=255)                                                                                                                                        gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcgc  <  1:715589/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                   GGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913/1209465‑1209748

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.