New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	3 (0.070)	21 (0.550)	20/250	1.0	86.3%	coding (305/630 nt)	stfP	e14 prophage, uncharacterized protein
	?	NC_000913	= 1209602	4 (0.100)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAC  >  NC_000913/1209461‑1209739                                                                                                                                              |                                                                                                                                          aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCCGAGc                                                                                                                                          <  1:645548‑M2/143‑2 (MQ=255)     ggCAGTTGATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAAGAGGGAATGCAGAGccttc                                                                                                                                      <  1:710945‑M2/143‑6 (MQ=255)       cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                    >  2:645548‑M2/1‑136 (MQ=255)                         aGCCTTTGCCAGTTCCGGGTATCCTGCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                  >  2:307084‑M2/1‑118 (MQ=255)                           ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                <  2:667870‑M2/143‑28 (MQ=255)                           ccTTTGCCAGTTCCGGGTATGCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                >  2:391431‑M2/1‑116 (MQ=255)                             tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                                                                                              <  2:692739‑M2/143‑30 (MQ=255)                             tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACAACAACAGGTAATGCAGAGCCTTCTCaccaaccacagctttcgatagg                                                                                                              >  1:242340/1‑129 (MQ=255)                                    gtgccGGGGACTCTTCAGCAGAAAAAGCGGCACCAGTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                                                                       <  1:674640/140‑1 (MQ=255)                                                     gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                      >  2:461568‑M2/1‑90 (MQ=255)                                                           aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc                                                                                <  2:220220‑M2/143‑60 (MQ=255)                                                                       ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagggtttcgagagccgtttgcacagtgccatccgatttgatctcgccaaacggatt                                                                    <  1:57704‑M2/143‑72 (MQ=255)                                                                        tGCATTTCAGCCAGCCTGTTGGCAGAGTGCCTCAACGCCTCTGTACCTGGACACCAACATGTGATGCAGAGCCTTCTCCCCAACCAACGTGTCTGAAAACGACGAGACAACTAGCAAATGGCATCATTCATGCTTTCACCAgg                                                                   >  1:269454/1‑143 (MQ=255)                                                                           aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGCGACACCAACAGGTAATGCAGAGccttctcccgaaccaagggtttcgagagccgtttgcacagtgccatccgagttgacatcgccacaccgagtgtgg                                                                <  1:258299‑M2/143‑76 (MQ=255)                                                                            tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                               <  2:467852‑M2/143‑77 (MQ=255)                                                                                     gccTGTTGGCGGTGTGTCTGCACCCCCCAGAACCGGGACACCAACAGGTAATGCAGAGccttctcccacaccaaggttttcgagagccgtttgcaccgtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                      >  1:486195‑M2/1‑58 (MQ=255)                                                                                             gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagcgacgataatcatcgcccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                              <  2:173344‑M2/143‑94 (MQ=255)                                                                                                gAGTGGCTGAAGGCCAGGGCTCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                                           >  1:566747‑M2/1‑47 (MQ=255)                                                                                                              cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattggatatcgccaaacggctccttccggcctaacgccagcgcgcgcagtgcggtaagtca                             <  2:315657‑M2/143‑111 (MQ=255)                                                                                                               aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagcgccatccgatttgatatcgccaaacggattcttgcggcttaaccgccgagcgcggagtgcggtcagtacc                            <  1:44151‑M2/143‑112 (MQ=255)                                                                                                                  gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTCTCATGCTTATCGACTATGTACCCGTATCAACaaa                         >  2:50487/1‑143 (MQ=255)                                                                                                                   aaacgggcaagcAACAGGTAATGCAGAGCCTTCTCACAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                        <  1:270847/132‑1 (MQ=255)                                                                                                                      cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgcagtgcggtaagtaactggtcgt                     >  1:667870‑M2/1‑25 (MQ=255)                                                                                                                        ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                   <  1:676285‑M2/143‑121 (MQ=255)                                                                                                                             tcaccttgtattgcaGAGccttctcccaaaccaagtttttcgagagccgtttgcacactgccatccgattcgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc              >  1:80405‑M2/8‑18 (MQ=255)                                                                                                                              cACCAGGTAATGCAGAGccttctcccaagccaccgttttcgagggccgtgtgcacggtgccatccgatttgatagcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc             >  1:657288‑M2/1‑17 (MQ=255)                                                                                                                                aCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccctt           <  2:436807‑M2/143‑129 (MQ=255)                                                                                                                                aCAGGTAATGCAGAGccttctcccaaaccaaggtgttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccctt           >  2:30372‑M2/1‑15 (MQ=255)                                                                                                                                         gCAGAGccttctcccagaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggct  >  2:618191‑M2/1‑6 (MQ=255)                                                                                                                                              |                                                                                                                                          ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAC  >  NC_000913/1209461‑1209739

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.