New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	3 (0.150)	10 (0.560)	8/252	0.9	77.9%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	3 (0.170)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGC  >  NC_000913/1207649‑1207931                                                                                                                                              |                                                                                                                                              tggtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGGGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                                                                                              <  1:381503‑M1/143‑1 (MQ=255) tggtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                                                                                              <  1:226928‑M1/143‑1 (MQ=255)  gcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttattcattttcataaacgTTg                                                                                                                                             <  2:261214‑M1/3‑1 (MQ=255)         gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAACTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCACGGTTTGGG                                                                                                                                      >  2:359753‑M1/1‑143 (MQ=255)               cattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcctgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                                                                                <  1:215229‑M1/16‑1 (MQ=255)                attacgttgtagatatgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaataactggatctattgcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTAGTTTGGGAGAAGGa                                                                                                                               <  2:309215‑M1/17‑1 (MQ=255)                   ccAGCCTGGAGAAGGGGCGGCACGACCAGTTATTTACCTCAGGTCGCGGGCTGCTGTTAAGCCGCAAGAATCCGTTTGTCGATATCAAATCGGATGGCACTTTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCa                                                                                                                            >  2:334048/1‑143 (MQ=17)                    cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                                                                                                           <  2:379460‑M1/143‑6 (MQ=255)                     ggcCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCACTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGCTATCTGTGTTTTGGTAGAAGGctctgc                                                                                                                          <  2:228294‑M1/142‑7 (MQ=255)                     aggctgtAGAGGGGGCGGCACGACCAGTTACTTACCGGACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                                                                                          <  2:345562‑M1/140‑7 (MQ=255)                         tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGCAGAAGGgcctgcatta                                                                                                                      >  2:17620‑M1/1‑133 (MQ=255)                           gagagtgggCTGTACGACAAGTTACTTACCGCACTTCGCGCGCGGCTGTTAAGCCGCAAGAATCCGTTTGGCGATTTCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                    <  2:177595/137‑1 (MQ=18)                                 gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATACCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGctatgcagtaccagttgg                                                                                                              >  2:96110‑M1/1‑125 (MQ=255)                                  gtcggcacGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAATCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCACGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                                                                                             <  1:7572‑M1/141‑20 (MQ=255)                                               gttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATtat                                                                                                <  1:385635‑M1/48‑2 (MQ=255)                                                           acatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttattcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCATCTACGGCTCAATTAGGACGCAccggct                                                                                    <  2:259031‑M1/60‑3 (MQ=255)                                                                      tgctgTTAAGCCGCAAGAATCCGATTGGCGATGGCGAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaaca                                                                         >  2:335348‑M1/1‑88 (MQ=255)                                                                            tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgaattacctgatgatgtccagattccgtggccttcagacactacgcaaacagggtgg                                                                   >  1:21855‑M1/1‑82 (MQ=255)                                                                               aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGATTGGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGCTATg                                                                <  2:247070‑M1/80‑1 (MQ=255)                                                                                                     tatCACATCGGATGGCGCTCTGCATACGGCTCTCGAAAGCCTTGGTTTGGGAGAAGGctctgcattacccgttggtgtcccggttccgtggcctttagccactccgccaacaggctggctgaaatgcaatggtgccgcttttt                                          <  2:258163‑M1/143‑87 (MQ=255)                                                                                                        ggtgtgatgtgcttgttatttcttcatttgcataagcgTGGGTTTGGTAGAAGGAGCAAAACGCAATGCAGCAACGGCTATATTAGGACGGATCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGTGCAATCAGTaa                                       >  2:245088‑M1/39‑143 (MQ=255)                                                                                                                     tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                          <  2:42733‑M1/118‑1 (MQ=255)                                                                                                                           cAAACGGCTCTCGAAAACCTTGTTTTGGGAGAAGGctctgcattacctgtttttgtcccgtttccttggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaact                    <  1:17620‑M1/143‑109 (MQ=255)                                                                                                                              tcgtctctcgaaAACCTTGGTTCGGGAGAAGGctctgcattccctgttggtgtcccggttccgtggccttccgccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggc                 <  2:209185‑M1/139‑112 (MQ=255)                                                                                                                               tcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc                <  2:83720‑M1/128‑1 (MQ=255)                                                                                                                                 gCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaa              >  1:85245‑M1/1‑29 (MQ=255)                                                                                                                                      cataaaCCTTGGTTTGGGTGAAGGAGCAATACTCATTGTAGCAACGGCTATATTATGACGCACCGGTTTCATAGCTATACAGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAg         <  1:105805/140‑1 (MQ=255)                                                                                                                                           acgTTGGTGTGGCAGAAGGAGCAAAACTTGCTGCAGCAACGGCGACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAagcagct    >  2:136559‑M1/4‑140 (MQ=255)                                                                                                                                             gTTGGTTTTCAAGAAGAAGCAAAACTCAATGCAGTATCTGCTACATTAGGACGCACCGGTTTCATAGCTATAACGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc  >  1:123485‑M1/2‑143 (MQ=255)                                                                                                                                             gTTGGTTCGGGAGACGGAGCTAAACTCAATGCAGCACCGGCTACATTAGGACGCGCCGGTTTCATAGCTATACCGGTTATTATTTGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc  >  2:20293‑M1/2‑143 (MQ=255)                                                                                                                                              |                                                                                                                                              TGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGC  >  NC_000913/1207649‑1207931

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