New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	7 (0.230)	18 (0.680)	16/252	0.6	76.3%	coding (290/630 nt)	stfP	e14 prophage, uncharacterized protein
	?	NC_000913	1209619 =	5 (0.190)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

TAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGG  >  NC_000913/1209483‑1209761                                                                                                                                         |                                                                                                                                               tAAGCCGTGGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                          <  2:32300/143‑1 (MQ=255)              tccatccccactgcatgattactgattgctcaccaccaccattcatacccggtatagctatgaaaccgttgcgtcctaatgtagccgttgctgcattgagtttagctcctgctcccaaaccaaCGTTTATGAAAATGAAGaaa                                                                                                                             <  1:520523‑M2/20‑1 (MQ=255)                   ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgtcgctgcatggagttttgctccttcccccaaaccaaCGTTTATCAAAATGAAGAAATACCa                                                                                                                        >  1:473945‑M2/119‑143 (MQ=255)                    gttATTCCTCAGCATAAAAGGCGGGGCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                       <  1:317617/141‑1 (MQ=255)                            cggacttctgtttgcgcaatagcatcaatcataaccggtatagctatgaaaccggtgcgtcctaatgtatccgttgctgcattgagttgtgctccttctcccaaaccaaCTTTTATGAAAATGAAGAAAGAACAAGCCAATgg                                                                                                               <  1:333970‑M2/34‑1 (MQ=255)                              gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                             >  2:77982‑M2/108‑143 (MQ=255)                                        tgctcaataccaccaatcagaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                   <  2:8146‑M2/46‑1 (MQ=255)                                          ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagctctgctccttctcccaaaccaaCGTTTATGAAAATGAAGGAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                                 >  1:383865‑M2/96‑143 (MQ=255)                                           tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGATTATGAAAATGAAGAAATAACAAGCAAATGGAATCATTACTgctat                                                                                                >  1:284822‑M2/95‑141 (MQ=255)                                                    ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                       >  2:464489‑M2/86‑143 (MQ=255)                                                    ccaatcataaccggtatagctatgaaaccggtgcgtactaatgtagccgttgctgaattgagttttgctccttctcccaaaccaaCGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTCCTTTTACCAgggg                                                                                       >  2:339448‑M2/86‑143 (MQ=255)                                                          ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAAACCTTCTTTTACAAGGGGGATTTa                                                                                 >  1:339264‑M2/80‑143 (MQ=255)                                                            aaccggtatagctatgaaaccggagcgtcctaatgtagccgttgctgcattgagttttgctccttctaacaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAGTCCTGCTGTTACCAGGGGGGTTtaag                                                                               >  1:426844‑M2/78‑142 (MQ=255)                                                                             aaccggtgcgtcctaatgtagccgttgctgcattaaattttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTACTATTACCAGGGGGATTTAACATACTTATTAGCtatat                                                              >  1:204931‑M2/61‑141 (MQ=255)                                                                               gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGcgacgaac                                                            >  1:244676/1‑136 (MQ=255)                                                                                 ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGAGAATTAACAAACTTATTGGAAATGTAagc                                                          >  2:917‑M2/57‑143 (MQ=255)                                                                                      gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTGCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                     >  2:365831‑M2/52‑143 (MQ=255)                                                                                                       tctgcattgagttttgctcctactcccaaaccaaCGTCTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTACCAGGGGGAGTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                    <  1:161505‑M2/109‑1 (MQ=255)                                                                                                         tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGATTTAACATGCTTATTGGCTATGTACGCGTAGCAACAAATGACCAGAAcaca                                  >  1:520401‑M2/33‑143 (MQ=255)                                                                                                         aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTCATTGGCTACGTACGCGTATCAACAAATGACCAGAAcaca                                  >  1:270386/1‑143 (MQ=255)                                                                                                         aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATATAACATGCTTATTGGCTATGTACACGTATCAGCAAATGACCAGAAcaca                                  >  1:429722/1‑143 (MQ=255)                                                                                                             ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTCACGAGACCGAGGGAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACACGACCACAGATc                              >  2:119218‑M2/29‑143 (MQ=255)                                                                                                                 gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTCGCTATGTACGCGTAGCaacaaatgaccagacaacatctctcca                          >  1:226399‑M2/25‑130 (MQ=255)                                                                                                                                cccaaccaaCCGTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAACGACCAGCACACAGATCTACAACGTAATGCGCTGaa           <  1:104584‑M2/134‑1 (MQ=255)                                                                                                                                ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa           >  1:92394‑M2/10‑143 (MQ=255)                                                                                                                                         cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAgg  <  2:244676/143‑1 (MQ=255)                                                                                                                                         |                                                                                                                                               TAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGG  >  NC_000913/1209483‑1209761

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