New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	1207790 =	19 (0.810)	5 (0.240)	5/254	2.1	20.1%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_pgi	1209619 =	23 (1.100)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913_3_pae_pgi/1209487‑1209759                                                                                                                                     |                                                                                                                                             ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAAACGACACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                                                                                              >  2:311231/1‑149 (MQ=255)       ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGGCACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                                        >  2:48870/1‑149 (MQ=255)       ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGACGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACAAACGTTTATGAAAAGAAAGAAATaa                                                                                                                        >  2:327008/1‑149 (MQ=255)          gTTCGGGGTATCGTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  1:106759/149‑1 (MQ=255)          gTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  1:227623/149‑1 (MQ=255)                       ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                                        >  1:327814/1‑149 (MQ=255)                            cagAAAAAGGGGCAGCATTGCATTTCAGCCAGCGTGTTGGCGGAGTGGCTGAAGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                   <  2:8205/149‑1 (MQ=255)                               aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCACCAGGTAATGCAGAGCCTTCTCGCAAACCAACGTTTATGAAAATGAAGAAATGACAAGCAAATGGCATCATTCCTGCt                                                                                                >  1:366878/1‑149 (MQ=255)                                 aaggaggCACCGTGTCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                              <  2:147640/144‑1 (MQ=255)                                  attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                             >  1:174023‑M2/100‑149 (MQ=255)                                         aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  2:191962‑M2/57‑1 (MQ=255)                                              tggttttCCGTCATCCGGTTGCGGGAGTGGCTGAAGGCATCGAAACCGGGAAGCCAACAGGTAGTGCAGCGCCGTCTCCCAAATCAACGTTTATGAAAATGAAGAAATACCAAGTAGGTGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                                 <  2:25852/145‑1 (MQ=255)                                                     cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCAAAACCAACGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                          <  2:354547/149‑1 (MQ=255)                                                         cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                                      <  1:191947/149‑1 (MQ=255)                                                              tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtcgcc                                                                 <  2:316720/149‑5 (MQ=255)                                                              tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgcc                                                                 >  2:402315/1‑147 (MQ=255)                                                                   ctatgcaaccggtgcgtcctaatgtagacgatgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                            >  1:400262‑M2/67‑149 (MQ=255)                                                                          accggtgcgtcctaatgtagccgttgctgcaacaattttaacgccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                     >  1:99512‑M2/60‑149 (MQ=255)                                                                           gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCAAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                    >  2:47502/1‑149 (MQ=255)                                                                            cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                   <  1:314161/149‑1 (MQ=255)                                                                             tGAAGGCCACGGAACCGGGACAGCAACAGGTAATGTAGAGCCTTCTCCCCAACAAACGTTTATGAAAATGAAGACATAACAAGCAAATGGCATCATTCCTGCTTTTAACAGGGGGGTTTAACATTGTTATTGGCTATGTAggcgagtca                                                  >  2:236971/1‑144 (MQ=255)                                                                                    cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTACGCGTATCAACAAATg                                           >  2:182376‑M2/50‑149 (MQ=255)                                                                                           ccGGGACACCAACAGGTAATGGAGAGCGGTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                    <  2:171265/149‑1 (MQ=255)                                                                                            cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGGTTATGAAAATGAAAAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAAGCGTATCCACAAATGACCACAac                                   >  2:32093/1‑149 (MQ=255)                                                                                             aggACACCAACAGGTAATGCAGAGCCTTCTCACAAACCAACGTTTATGAAAACGAAGAAATAACAAGCATATGGCATCATTCCTGCTTTTACCAGGGAGATTTAACATGCTTATTGGCTATGTCCGCGTATCCACAAATGACCAAAaca                                  >  2:91535/2‑149 (MQ=255)                                                                                                    cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                           >  1:223485/1‑149 (MQ=255)                                                                                                         ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                      <  2:223485/149‑1 (MQ=255)                                                                                                            aaTGCAGAGCCTTTCCCCAAACAAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAACATTCCTGCTTGTTCCAGGGGTATTTAACATGCTTATTGGCTATGTATGCGTATTAACAAATGACCAGAACACAGATCTACAACGTa                   <  1:161954/149‑1 (MQ=255)                                                                                                                       ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc        >  1:196959‑M2/15‑149 (MQ=255)                                                                                                                             caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  2:366878‑M2/141‑1 (MQ=255)                                                                                                                                     |                                                                                                                                             CCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913_3_pae_pgi/1209487‑1209759

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