New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	1203246 =	5 (0.150)	24 (0.780)	23/292	0.3	89.7%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	CP009273	1205075 =	1 (0.030)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGC  >  CP009273/1204837‑1205342                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                            tACCCCTGAATACTTAAAATAGAGCGCCCTGTCTCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAATTCAGTCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                                                                                                                                                         <  1:309984/275‑1 (MQ=255)      tgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaatacccccaatcatcaccggtatagctatgaaaccggtgcgtcctcatgtagccgttgctgctatgagttttgctccttctcccaaaccacCGCTCAAGAAAACGAAGAAATAACAACCAcatgccataaatc                                                                                                                                                                                                                                    >  2:604351‑M2/234‑270 (MQ=255)                                                            aaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                                                                                                                                                                                                                        >  1:525917‑M2/180‑201 (MQ=255)                                                                gggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccaaccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                                                                                                                                                    <  2:525909‑M2/26‑1 (MQ=255)                                                                                                                                       tactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                                                                                                                                                             >  2:219597‑M2/105‑185 (MQ=255)                                                                                                                                           gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                                                                                                                                                         <  1:219599‑M2/85‑1 (MQ=255)                                                                                                                                                    ataccaccaatcgtaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                                                                                                                                                          >  2:129433‑M2/92‑111 (MQ=255)                                                                                                                                                        caccaatcgtaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                                                                                                                                                      <  1:129434‑M2/24‑1 (MQ=255)                                                                                                                                                             atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                                                                                                                                                                                        >  2:119623‑M2/83‑120 (MQ=255)                                                                                                                                                               cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                                                                                                                                                 >  1:489989‑M2/81‑205 (MQ=255)                                                                                                                                                                 taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                                                                                                                                                    <  1:119624‑M2/42‑1 (MQ=255)                                                                                                                                                                  aaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAg                                                                                                                                                                     >  2:81323‑M2/78‑182 (MQ=255)                                                                                                                                                                   accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcg                                                                                                                                             <  2:489982‑M2/129‑1 (MQ=255)                                                                                                                                                                    ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                                                                                                                                                       >  2:536535‑M2/76‑130 (MQ=255)                                                                                                                                                                      ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                                                                                                                                                 <  1:81324‑M2/109‑1 (MQ=255)                                                                                                                                                                        tatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                                                                                                                                                   <  1:536543‑M2/59‑1 (MQ=255)                                                                                                                                                                                    tccggtgtgtcctaatgtagccgttgctgcattgagttttgctccttctccaaaaccaaCGTTTATAAAAATGAAGAAATACNAAGCAAATNGAATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATGGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAACGTCCGAAGGGCCGGGACTGAAAAAACTGCTCAgg                                                      <  1:460005‑M2/216‑1 (MQ=255)                                                                                                                                                                                                            tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                                                                                                                                                         >  2:463225‑M2/36‑104 (MQ=255)                                                                                                                                                                                                            tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                                                                  >  1:288123‑M2/36‑143 (MQ=255)                                                                                                                                                                                                                gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                                                                                                                                                                     <  1:463232‑M2/73‑1 (MQ=255)                                                                                                                                                                                                                gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                                                                                                                                              <  2:288119‑M2/112‑1 (MQ=255)                                                                                                                                                                                                                      agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATTTAcgcg                                                                                                                                                                                      >  1:260282‑M2/26‑113 (MQ=255)                                                                                                                                                                                                                       gttttgctcattctcccaaaccaaCGTTTAAGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTACCAGGGGGATTTAACANGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAACCTGCTCAGGACATTAGCGGCAGGTGACACTCTGGTTGTCTGGaa                   <  1:502856‑M2/251‑1 (MQ=255)                                                                                                                                                                                                                          ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                                                                                                                                                  <  2:260278‑M2/92‑1 (MQ=255)                                                                                                                                                                                                                           tgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                                                                                                                                            >  2:324702‑M2/21‑118 (MQ=255)                                                                                                                                                                                                                               ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                                                                                                                                                        <  1:324706‑M2/102‑1 (MQ=255)                                                                                                                                                                                                                                              aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                                                                  >  2:327713‑M2/2‑109 (MQ=255)                                                                                                                                                                                                                                              aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGCTTTTTGAAGACAAGATAAGCGGCACAAAGTCAGAAAGGCCGGGACTGAAAAAACTTACAAGGACATTATCTTCAGGTCACACGCGGGttgtctgacagctgtttcgcctctttc  >  2:336090‑M2/2‑256 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                            TACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGC  >  CP009273/1204837‑1205342

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