New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	31 (0.810)	5 (0.150)	5/252	4.2	14.4%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	32 (0.940)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGC  >  NC_000913/1207652‑1207925                                                                                                                                           |                                                                                                                                        tggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                                                                                     >  1:398065/1‑143 (MQ=255)      gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTggg                                                                                                                                <  2:252478/143‑1 (MQ=255)       cATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGGTATCAAATCGGATGACACTGTGCAAACGGCTCTCCAGAACCCAGGTTCGGga                                                                                                                               >  2:96372/1‑143 (MQ=17)        aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGCAAACCTTGGTTCGGgag                                                                                                                              >  2:352141/1‑143 (MQ=255)        aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAACCGCAAGAATCCGTTTGGCGATCTCAAATCGAATGGCAATGTGCAAACGGCACTCTGAAACCTCGGGTTCGgag                                                                                                                              >  1:340557/1‑143 (MQ=17)          ccGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGCTTGGCGATACCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGaa                                                                                                                            >  1:695763/1‑143 (MQ=255)                        gagggggggCTGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                              <  1:152789/138‑1 (MQ=255)                        gagaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                              <  2:688818/143‑1 (MQ=255)                             ggggggGCACGACCAGTTACGTACCGCACTGCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCACTgcag                                                                                                         <  2:579892/143‑1 (MQ=255)                                  ggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAACACTCTATGCAGCAACg                                                                                                    >  1:203510/1‑143 (MQ=255)                                    cgcGACCAGTGGCCGACCGGACTTCGCGCGCTGCTGTCAAAGCACCAGAATCCGTTTGGCGATATCAAATCGGACGGCACTGTGCAACCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                                                  <  1:220301/141‑1 (MQ=18)                                    cACTACCTGTTACTTACCGGACTGCGCGCGCTGCTGTTAGGCCGCAAGAATCCGTTTGGCGAGATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGGGAGGTAGCAAAACTCAATGCAGCAACGGc                                                                                                  <  2:499429/143‑1 (MQ=21)                                           acttacttacCCCACTCCGCGCGCTGCTGTCACTCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTCCAAACGGCTCTCTAAAACCTTGGTTTGGGAGAAGGAGCAAAACTAAATGCAGCAACGGCTACATTa                                                                                           <  2:28279/141‑1 (MQ=21)                                             ttactgacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                                                                         <  1:282382/143‑1 (MQ=255)                                                    ccGCACTTCGCGCGCTGCTGTTAAGACGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTAAAAACGGATGACGAAAACGATGGTTGGGGAGCAGGAGCAAACATCAATGAAGCAACAGCTAAATTAGGGAGCAcg                                                                                  >  2:107453/1‑142 (MQ=18)                                                             gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                         <  1:556071/143‑1 (MQ=255)                                                              cgcgcTGCTGTTAAGCCGCAAGAATCCGTTAGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTGGGGAGAAGGAGCAAAACTCAAGGCAGCAACGGCTAGATTAGTACGCACCGGTATCGTAg                                                                        >  2:498513/1‑143 (MQ=21)                                                                 gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAAAGGTTCTCTGAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTTATAGCta                                                                     >  2:646336/1‑143 (MQ=25)                                                                        ttgaggcgcgAGTATCCGGTTGTGGATATCAAATCGGAGGGCACAGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGACCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                              <  2:185826/137‑1 (MQ=25)                                                                                 aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                     <  1:607318/143‑1 (MQ=255)                                                                                         gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAGGATTGGTGGTATTg                                             >  1:394638/1‑143 (MQ=255)                                                                                              gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                                        >  1:696134/1‑143 (MQ=255)                                                                                               ggATATCAAGTCGGGTGGCACTGGGCGAACGGCTCTCGAAAACCTTGGTTTGGGCGAAGGAGCTAAACTCAATGCAGCAACGGCTACATTAGGACGCATCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                       <  2:418629/142‑1 (MQ=255)                                                                                                 ctatCAAATCCGATTGCACTGTGCAAACCGCTCACGAAAGCCGTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                     <  2:111320/142‑1 (MQ=255)                                                                                                 atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTCTACGGGTTATGATTGGTGGTGTTGGGCGAtga                                     >  1:253182/1‑141 (MQ=38)                                                                                                      aaaTCGGATGGCACTGTGCAAACGGCTCTCGCAAACCTTGGTTTGGGAGAAGGAGCAAACCTCAATGCAGCAACGGCTACGTTAGGACGCACCGGTTTCAGAGCTATACCGGTTAGGACTGGAGGTATTGCGCCATAAGTGAt                                >  2:46180/1‑143 (MQ=37)                                                                                                                      tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAGCATTCAGTgggggggg                >  2:648115/1‑143 (MQ=255)                                                                                                                          aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGGTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATg            >  2:537584/1‑143 (MQ=255)                                                                                                                           accgCTCTCCAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGc           <  2:535834/140‑1 (MQ=255)                                                                                                                                tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGCGGAATGCCGCaa      >  2:257554/1‑143 (MQ=255)                                                                                                                                    gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCCTTCAGTGGGGGTGGAATGCCGCAAAAGc  >  1:400517/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                        TGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGC  >  NC_000913/1207652‑1207925

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.