New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	1207790 =	6 (0.200)	25 (0.960)	22/250	0.3	85.7%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_pgi	1209619 =	3 (0.120)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_pae_pgi/1209479‑1209764                                                                                                                                             |                                                                                                                                                  agatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                                                                                                           <  1:198542‑M2/9‑1 (MQ=255)  gatgctttcgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                                                                                                          <  1:64420‑M2/10‑1 (MQ=255)     tAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                                                       <  1:211267/149‑1 (MQ=255)                  cccacccccactgaatgattactaatggcgcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                          <  1:224995‑M2/26‑1 (MQ=255)                                    ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                        <  2:274830‑M2/44‑1 (MQ=255)                                        tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtaaccgttgctgcattgagatttgctccttctcccaaaccaaCGTTTATGAAAATAGAAAAATAACAACCAAAAGACATCATTCCTGCtt                                                                                                    >  2:128669‑M2/102‑149 (MQ=255)                                               tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             >  1:41511‑M2/95‑149 (MQ=255)                                               gCAACATTGCATTTCAGCCAGCCTGTTGGGGGGGTGGCTGAAGGCCACGGAACGGAGACACCAACCTGTAATGCCGAGTCTTCTCCCAGACCAACGATTATGAAACTGTAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAg                                                                                             >  2:85546/1‑149 (MQ=255)                                                caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                            >  2:444553‑M2/94‑149 (MQ=255)                                                 aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                           <  1:179714‑M2/57‑1 (MQ=255)                                                 aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                           >  1:127058‑M2/93‑149 (MQ=255)                                                           atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                 <  2:306388‑M2/67‑1 (MQ=255)                                                                         agctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                                   <  1:317589‑M2/81‑1 (MQ=255)                                                                           ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                 >  1:265467‑M2/67‑149 (MQ=255)                                                                             ggggTGGCTGAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                               <  2:440098/146‑1 (MQ=255)                                                                                      gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAGGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACCTGCTTATGGGCTATGTACGCGTATcaa                                                      <  2:114728‑M2/94‑1 (MQ=255)                                                                                              taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                              <  2:262704‑M2/102‑1 (MQ=255)                                                                                               aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                             <  2:41511‑M2/103‑1 (MQ=255)                                                                                                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGAATCATACCTATCAGTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                         <  2:110674‑M2/107‑1 (MQ=255)                                                                                                        gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                    >  1:365380‑M2/38‑149 (MQ=255)                                                                                                        gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACCCGTATCAACAAATGACCAGAACACAg                                    >  1:472293‑M2/38‑149 (MQ=255)                                                                                                          tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                  <  2:127058‑M2/114‑1 (MQ=255)                                                                                                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAATTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAACt                                >  1:136356‑M2/34‑149 (MQ=255)                                                                                                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATCGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACa                             >  2:410632‑M2/31‑149 (MQ=255)                                                                                                                            tccttctcccaagccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGAATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCTTATGAAAAAATGACCAGAACACAAGTCTACAACGTAttgcgcgg                >  2:74528‑M2/18‑147 (MQ=255)                                                                                                                            tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACTTAATGCGCTg                >  2:442582‑M2/18‑149 (MQ=255)                                                                                                                            tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                 >  2:139652‑M2/18‑148 (MQ=255)                                                                                                                            tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                 <  1:139652‑M2/131‑1 (MQ=255)                                                                                                                               ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGTAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc             >  1:449883‑M2/15‑149 (MQ=255)                                                                                                                                  tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt          <  2:75774‑M2/138‑1 (MQ=255)                                                                                                                                          caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTCAAACGTAATGCGCTGAactgtgcgagagg  >  1:291246‑M2/4‑143 (MQ=255)                                                                                                                                             |                                                                                                                                                  CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_pae_pgi/1209479‑1209764

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