New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	50 (0.510)	13 (0.150)	13/250	5.6	22.2%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	= 1209602	47 (0.540)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

TCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAA  >  NC_000913/1207660‑1207949                                                                                                                                                  |                                                                                                                                                 tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                               <  3:377956/149‑1 (MQ=255)  ccGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctct                                                                                                                                              >  7:92699/1‑145 (MQ=255)   cGGTGCCAGCCTGGAGAAGTGGCGGAACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCTCAAGAATCCGTTTGGCGATATCAACTCGGCTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGgcgctg                                                                                                                                             >  2:245927/1‑147 (MQ=17)         cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGTATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGTGAGAAGGAGCAAAACTc                                                                                                                                        >  8:100601/1‑148 (MQ=255)            ccTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                                                                                                    <  5:353339/149‑1 (MQ=255)             cTGGAGAAGCGGCGGCACGACAATTTTCTTACCTGACTTCGCGCGCTTCTGTTAATCCGCAAGAATCCGTTTGTCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                                   <  4:296529/149‑1 (MQ=18)                gagtAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCAGCAAGAATCCGTTTGTCGATATCAAATCGGATGGCACTGTGCAATCGTCTCTCTAAAACCTTGGTTTGTGAGTAGGAGCAAAACTCTATgcag                                                                                                                                 <  8:57722/148‑1 (MQ=255)                  gaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                                   >  2:10630/1‑144 (MQ=255)                  gaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                                   <  1:10630/144‑1 (MQ=255)                     ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                           <  5:34235/149‑1 (MQ=255)                               gACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCACAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                                                                                                                  >  5:282773/1‑148 (MQ=255)                                 acAGTTAGTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAATAATCCGTTTGGCGATATCAAATCGGATTGCACTGTGCAAACGGCTCTCGAATACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                                                                          <  5:133389/137‑1 (MQ=255)                                           acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                                     >  5:26313/1‑149 (MQ=255)                                           acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                                                                                       <  8:575/147‑1 (MQ=255)                                            ccGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCCCCGGttt                                                                                                     >  6:356600/1‑148 (MQ=255)                                              gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                                  <  5:356600/149‑1 (MQ=255)                                                 cTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                                                               <  2:351519/149‑1 (MQ=255)                                                 cTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                                                               <  6:118010/149‑1 (MQ=255)                                                     gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                                                                                           >  1:211666/1‑149 (MQ=255)                                                               gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                                                                 >  5:160560/1‑149 (MQ=255)                                                                ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                                                                                <  5:124189/149‑1 (MQ=255)                                                                  aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                                               <  7:142359/148‑1 (MQ=255)                                                                   aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtg                                                                              <  3:244285/148‑1 (MQ=255)                                                                   aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGAAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                        >  8:29326/1‑138 (MQ=255)                                                                   aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGAAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                        <  7:29326/138‑1 (MQ=255)                                                                    gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                            <  4:186523/149‑1 (MQ=255)                                                                     ccacAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                                           <  6:282773/146‑1 (MQ=255)                                                                     ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                      <  1:237498/90‑1 (MQ=255)                                                                     ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                      >  2:237498/1‑90 (MQ=255)                                                                       gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg                                                                          >  4:10598/1‑148 (MQ=255)                                                                        cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGAACCGGTTTCTTAGCTATACCGGTTATGATTGGTGGTa                                                                        >  5:73094/1‑149 (MQ=255)                                                                         aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                                       >  5:167227/1‑149 (MQ=255)                                                                              tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                                                                        >  2:374918/1‑143 (MQ=255)                                                                              tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                                                                        <  1:374918/143‑1 (MQ=255)                                                                                cGTTTGGCGATATCAAATCGGATGGCCCTGTGCAAACGGATCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACCTTAGGACGCCCAGGTTTCATATCTCTACCGGTTATGATTGGTGTGAGTGAGCaa                                                                >  5:277207/1‑149 (MQ=37)                                                                                     ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                                            <  3:10598/148‑1 (MQ=255)                                                                                          tatCAAATCGGAGGTCACTGTTCAAACGGATCTCGAAAACCTTGCTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGTCTACATTAGGACGCACCGGTTTCATCGCTATACCGGTTATGATTGGTGGTATTGATCAATCAGTAATc                                                       <  6:277207/148‑1 (MQ=255)                                                                                                   ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTGTACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                                              <  7:100601/148‑1 (MQ=255)                                                                                                      tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGt                                          >  3:92647/1‑149 (MQ=255)                                                                                                         cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAAGACGCACCGGTTTCATAGCTATACCGGtt                                                                                     <  2:340809/103‑1 (MQ=255)                                                                                                         cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAAGACGCACCGGTTTCATAGCTATACCGGtt                                                                                     >  1:340809/1‑103 (MQ=255)                                                                                                                cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCg                                 >  3:124174/1‑148 (MQ=255)                                                                                                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                            <  6:26313/149‑1 (MQ=255)                                                                                                                        tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                             >  6:68394/1‑128 (MQ=255)                                                                                                                        tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                             <  5:68394/128‑1 (MQ=255)                                                                                                                                   ttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                            <  3:97944/118‑1 (MQ=255)                                                                                                                                   ttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                            >  4:97944/1‑118 (MQ=255)                                                                                                                                           ggAGAAGGATCAAAACTCAATGCATCAACTGCTACATTAGGACGCACCGGTTTCATACCTATACCGGTTATGATTGGGGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATgg     <  6:167227/149‑1 (MQ=255)                                                                                                                                           ggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCTGTTATGATTGGTGTTATTGAGCAATCAGTAATCATTCAGTGGGGCTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATgg     <  2:211666/149‑1 (MQ=255)                                                                                                                                           ggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAACATACAGTGGGGGTGGAAGGCCGCAAAAGAATCTGCCTCTGGGGGGGATg      >  5:234793/1‑148 (MQ=255)                                                                                                                                              gaAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaaa  >  5:99223/1‑149 (MQ=255)                                                                                                                                               aaGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaaa  >  8:59477/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAA  >  NC_000913/1207660‑1207949

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