New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	3 (0.110)	20 (0.770)	20/566	0.3	85.4%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	= 1209602	4 (0.150)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AAAATAGAAATTAATCGTTGCATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGC  >  NC_000913/1209305‑1209895                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                      aaaaTAGAAATTAATCGTTGCATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTTGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                                                                                                                                                                                      <  2:41307‑M2/299‑2 (MQ=255)                                aTCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGCTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                                                                                                                                                                                     >  2:36069‑M2/1‑267 (MQ=255)                                             aCGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGTGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggatttcgagagccgtttgcacagtgccat                                                                                                                                                                                                                                                         <  2:7044‑M2/300‑47 (MQ=255)                                                            cATGAGCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTTGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg                                                                                                                                                                                                                                           >  1:41268‑M2/1‑239 (MQ=255)                                                             aTGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagggtttcgagagccattggcacagtgccatccgatttgagatcgcca                                                                                                                                                                                                                                        <  1:73634‑M2/301‑64 (MQ=255)                                                                 ttcagacgtgtgctcttccgatctcttaaaatagagcgacctgtatcaaTCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac                                                                                                                                                                                                                                     >  2:108439‑M2/25‑234 (MQ=255)                                                                                      aTACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc                                                                                                                                                                                                                <  2:270789‑M2/300‑88 (MQ=255)                                                                                         cTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacagatcggaagagcgtcgtgtagggg                                                                                                                                                                                                            <  1:108331‑M2/301‑92 (MQ=255)                                                                                              aaTAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                                                                                                                                                                                                       >  2:73703‑M2/1‑205 (MQ=255)                                                                                               aTAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                                                                                                                                                                                      >  1:11921‑M2/1‑204 (MQ=255)                                                                                                       aCCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggta                                                                                                                                                                                              <  2:229406‑M2/301‑106 (MQ=255)                                                                                                                                      ccACGAATAAACTCCA‑CACGGAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca                                                                                                                                                                  >  1:270457‑M2/1‑165 (MQ=255)                                                                                                                                         cGAATAAAC‑CACACACGTAAATCAGTCAATTTATTTCCCCGAGAACCCTGTGCCATTTCCGTGTTTTCTGCAGCAGAAAAAGCGGGACCACAGCAGTTCAGCCAGCCTGTTGGCGGAGTGGCTGCAGGCCACGGAACCGCGTCACCAACAGGGAATGCAGAGcgttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccgtacggattcttgcggcttaaccgcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg                                                                                                                                                            >  1:7044‑M2/1‑162 (MQ=255)                                                                                                                                                    ac‑cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgtcttaacagcagcgcgcgaagtgtggtaagtaacttgtcgtgccgccccttctccagcctgccaccggatgc                                                                                                                                                  <  1:102445‑M2/300‑150 (MQ=255)                                                                                                                                                                                        ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGa                                                                                                              >  2:34963/1‑301 (MQ=255)                                                                                                                                                                                              ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgacttgatatcgccaaacggattcttgcggcttaacagccgcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcacaggatgcctccaccacgctgcaccgttccccttgccacatgtccaagt                                                                                                         <  1:111574‑M2/300‑191 (MQ=255)                                                                                                                                                                                                          aTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccag                                                                                            <  2:54882‑M2/301‑204 (MQ=255)                                                                                                                                                                                                             cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGATAACAGTCCGaaa                                                                                         >  2:255069/1‑301 (MQ=255)                                                                                                                                                                                                                                              gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtacgtaactggtcgtgtcgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgccaaagtagtcatcgtccagaccggtggcagccttgccggtttgggggttaccacg                                                        <  1:48895‑M2/301‑240 (MQ=255)                                                                                                                                                                                                                                                  gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggtctgggggttaccacgggta                                                    <  2:5447‑M2/301‑244 (MQ=255)                                                                                                                                                                                                                                                             gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggtctgggggttaccacgggtaaaaccgttct                                          >  1:5447‑M2/1‑47 (MQ=255)                                                                                                                                                                                                                                                               gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACAct                                          >  2:149719/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                 aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGAt                     <  1:149565/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                ggTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggtctgggggttaccaagggtaaaacctttcttacccgcgccgaacttacccttctgcgcggctttcg      <  1:199257‑M2/301‑290 (MQ=255)                                                                                                                                                                                                                                                                                                    aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcccagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggtctggggtttaccacgggtaaaaccgttcttacccgcgccgaacttatccttctgcgcggttttcgtgtc  >  2:48944‑M2/1‑8 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                      AAAATAGAAATTAATCGTTGCATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCAC‑CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGC  >  NC_000913/1209305‑1209895

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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.