New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	0 (0.000)	24 (0.750)	24/568	0.3	92.3%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	4 (0.120)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCT  >  NC_000913/1209325‑1209865                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ctacacacagtaatacaaacgaattggttggatatccccccgttgccactgcattcataggtgagctgacattgtcataagttgcaacaacggcaacacaggcattattaaacgcaaccaggaatacaacggtatttccatcccccccagaggcgggtgcttttgcgacattccacccccactgaatgattacagattgctcaataccgccagtcataaccggtatagctgtgaaaccggtgcgtcctaatgtagccgtggctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                                                                                                                                  <  1:304669‑M2/7‑1 (MQ=255)            aTCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCACGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCCATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGCTGGCGGAGTGGCTGACGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCATTCTCCCAAACACACGTTATAGCAAATgaaga                                                                                                                                                                                                                                       >  2:59767/1‑301 (MQ=255)                 aaacgaagtggttgtatatccccccgttgccactgcattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcactaccaccaatcataaccggtatagctatgaaacagttccgtcctatggtagccgttcctcctttgagttttgctcctccacccaacccaaCGTTTATGAAAATGAACAAATaa                                                                                                                                                                                                                                  >  1:92250‑M2/279‑301 (MQ=255)                  aacgaagtggttgtatatccccccgttgccactgcattgataggtgcgctgacattgtaataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcatttagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAtaaa                                                                                                                                                                                                                                 >  1:115621‑M2/278‑300 (MQ=255)                         aCGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                                                                                                                          >  2:87674/1‑301 (MQ=255)                                atatccccccgttgccactgcattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgatttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                                                                                                                                    >  2:73240‑M2/264‑300 (MQ=255)                                   tccccccgttgccactgcattgataggtgctctgtcattgtcataatttgcaactacggcaccacattcattattaaacgcaaccgggaatacaattgtattaccactccccccagaggcagatgctttttcggcattccacccccactgaatgattactgattgttcaataccaccaatcataaccggtatagctattaaaccgctgagctctaatgtagacgttgctgcattgagttttgctccttctcccaaaccacAGTTTATGAAAATGAAGCAATAACAAGCAAATGGCattatt                                                                                                                                                                                                                >  2:125497‑M2/261‑301 (MQ=255)                                        ccgttgccactgcattgataggggcgctgacattgtcagaatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                                                                                                                           <  2:262376‑M2/46‑1 (MQ=255)                                             gccactgcattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                                                                                                                                      >  2:19205‑M2/251‑301 (MQ=255)                                                                          gtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcacaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                                                                                                                                         >  2:221171‑M2/222‑301 (MQ=255)                                                                                  ttgaaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                                                                                                                                 <  2:115624‑M2/88‑1 (MQ=255)                                                                                     caacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaatcccaccaatcataaccggtatagctatgaaaccggtgcgtcctactgtagccgttgctgcatggcgttttgctccttctcccaaaccaaCGTTTATGAAATTGAAGAAATAACAAGCAAATGGCATTACTTCTCCTTTTCACAGGGGGATTACACATGCTTATTGCTTATGTAAGCTTAt                                                                                                                                                              >  1:226445‑M2/211‑301 (MQ=255)                                                                                      ccaaccacacccttccgatttattaaacgccaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccactcccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                                                                                                                             <  1:208593‑M2/92‑1 (MQ=255)                                                                                                           attaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACTAGGGGGTTTTACCATGCTTAttggctatgtccgcgtatcagatcggaagagaacacgt                                                                                                                                          >  2:208599‑M2/189‑281 (MQ=255)                                                                                                                     accgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCATCAAATGACCAGAACACAGATCTAAAACGt                                                                                                                              >  2:276946‑M2/179‑301 (MQ=255)                                                                                                                                     taattccatcccccccagaggcagatgctttttcgtcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgtg                                                                                                              <  1:19204‑M2/139‑1 (MQ=255)                                                                                                                                              cccccccagagtcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCg                                                                                                     <  2:92253‑M2/148‑1 (MQ=255)                                                                                                                                                                             ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCCGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGca                                                                       <  1:87671/300‑1 (MQ=255)                                                                                                                                                                                cccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAATACAGATCTCCAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAACGACACaaa                                                                   >  2:28997‑M2/120‑301 (MQ=255)                                                                                                                                                                                      tgattcACCAGAAAAAGCGGCAACGTTGCCTTTCAGCCAGCCGGTTGGCGGAGTGGCTGAAGGCCACAGAACCTGGACACCAACAGGTAATTCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAAGAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTGACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGa                                                             <  1:59765/298‑1 (MQ=255)                                                                                                                                                                                                            aaaccaccaatcataaacgtcattgctatgaaaccggggcgtcctaattgagccgttgctgcattgagttttgctccttcgcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTg                                       <  1:276938‑M2/210‑1 (MQ=255)                                                                                                                                                                                                             aaccacaaatcataaccgggatagctatgaaaccgttgcggcataatgtagccggtgctgcaatgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGc                                      <  1:221163‑M2/211‑1 (MQ=255)                                                                                                                                                                                                                   caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGa                                >  2:285936‑M2/85‑301 (MQ=255)                                                                                                                                                                                                                              ggttaagctatgaaaccggtgcttcctaatgtagccggtgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGc                      <  1:285928‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                                    gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGTTAAGCGGCCCAAAGTCCGAATGGCCGGGCCTGAAAAAACTGctcaggacattactggaaggcaa                >  1:153952‑M2/68‑289 (MQ=255)                                                                                                                                                                                                                                     ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGaca              >  2:7252‑M2/67‑301 (MQ=255)                                                                                                                                                                                                                                             ccggtgcgacctaatgtagccgttgctgaattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGtt      <  1:73238‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                                  gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGATTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCt  <  1:28996‑M2/247‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       CATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCT  >  NC_000913/1209325‑1209865

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