New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 1207805	6 (0.200)	31 (1.190)	28/250	0.1	88.2%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_pgi	= 1209602	3 (0.120)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAC  >  NC_000913_3_pae_pgi/1209461‑1209728                                                                                                                                              |                                                                                                                               aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                         <  1:75774‑M2/149‑8 (MQ=255)  tCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                                                        >  1:60196‑M2/1‑141 (MQ=255)  tCAGGCAATTTATTTGTCGGAAAAGCCTTTGCCCGTTCCGAGTATTCTTCATCAGAAAAAGCGGCACCAATGCATTTCAGCCAGCCTGTTGGCGTAATGGCTGAAGGCCGCGGAACCGGGACACCAACAGTTAATGCAGAGCCTTCTcc                                                                                                                        >  2:425559/1‑149 (MQ=255)         atttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                                                                                 >  2:39385‑M2/1‑134 (MQ=255)           ttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                               >  2:415036‑M2/1‑132 (MQ=255)                 gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                         >  2:256806‑M2/1‑126 (MQ=255)                 gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                         >  1:301875‑M2/1‑126 (MQ=255)                       tAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                   <  1:211267/149‑1 (MQ=255)                        aaGCGTTTTCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCCACCAAAAATGGAGAGccttctcccaaaccaaggttttcgagagcc                                                                                                  >  1:384255‑M2/1‑119 (MQ=255)                          gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                <  2:233475‑M2/149‑33 (MQ=255)                                 ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGTATCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaca                                                                                         <  1:385016‑M2/149‑40 (MQ=255)                                        cGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccat                                                                                  <  2:210457‑M2/149‑47 (MQ=255)                                         gggTATTCTTCAGCAGAAAAAGCGGCCCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                 >  1:37692‑M2/1‑102 (MQ=255)                                            tATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                                                                              >  1:210457‑M2/1‑99 (MQ=255)                                                  tcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGGGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgat                                                                        >  1:148317‑M2/1‑93 (MQ=255)                                                           aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACGGGGACACCAACAGGTAATGCAGAGccttcgcccaaaccaaggttttcgagagccgtttgcaaagtgccatccgatttgatatcgccaaa                                                               <  2:476857‑M2/149‑66 (MQ=255)                                                           aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaa                                                               <  2:192252‑M2/149‑66 (MQ=255)                                                                 gCAACATTGCATTTCAGCCAGCCTGTTGGGGGGGTGGCTGAAGGCCACGGAACGGAGACACCAACCTGTAATGCCGAGTCTTCTCCCAGACCAACGATTATGAAACTGTAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAg                                                         >  2:85546/1‑149 (MQ=255)                                                                         gCATTTCACCCAGCCTGTTGGCGGAGTGGGCGAAGGCCACGGACCCGGGACACCAACAGGTAATGCAGCGccttctcctaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc                                                 >  1:209900‑M2/1‑70 (MQ=255)                                                                           aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctt                                               >  1:106079‑M2/1‑68 (MQ=255)                                                                           aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctt                                               <  1:284966‑M2/149‑82 (MQ=255)                                                                              tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGACACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaac                                            >  1:196491‑M2/1‑65 (MQ=255)                                                                                  ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                        <  2:71967‑M2/149‑89 (MQ=255)                                                                                   cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcag                                       >  2:385016‑M2/1‑60 (MQ=255)                                                                                    agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatccccacatggagtcttgcggcttacaagcagc                                      <  2:214255‑M2/149‑91 (MQ=255)                                                                                       cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                   <  1:275234‑M2/149‑94 (MQ=255)                                                                                               ggggTGGCTGAAGGCGACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                           <  2:440098/146‑1 (MQ=255)                                                                                               ggAGTTGTTGAAGGTCACGAAACCGGGACACCAAATGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgc                           >  1:425559‑M2/1‑48 (MQ=255)                                                                                               gaagtGGCGAAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggatttttgcggcttaacagcagcgcgcgaagtgc                           >  2:304407‑M2/3‑48 (MQ=255)                                                                                                gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                          >  1:233475‑M2/1‑47 (MQ=255)                                                                                                 aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcgg                         >  2:132831‑M2/1‑46 (MQ=255)                                                                                                       tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagta                   <  1:293107‑M2/149‑110 (MQ=255)                                                                                                           ggtcacgGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactg               >  1:115153‑M2/4‑36 (MQ=255)                                                                                                                cGGAACCGGGACACTAACAGGCAATTCACAGccttatccgcaacgaaggttttcgagaggcgtttgcagagtgccatccgattttatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagttcggtaagtaactggtcgt          >  2:65197‑M2/1‑31 (MQ=255)                                                                                                                 ggAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtg         >  1:452608‑M2/1‑30 (MQ=255)                                                                                                                      cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcc    >  1:192252‑M2/1‑25 (MQ=255)                                                                                                                       gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggactcttgcggcttaccagcagcccgcgaactgcgataagtaacagatcgtgccgccc   <  1:76636‑M2/149‑126 (MQ=255)                                                                                                                        ggACACCAACAGGTAATGCAGAGccttctcccaaacccaggtattcgagagccggttgcccagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcccgcgaagtgaggtcagtaactggtcgtgccgcccc  <  2:409193‑M2/149‑127 (MQ=255)                                                                                                                                              |                                                                                                                               ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAC  >  NC_000913_3_pae_pgi/1209461‑1209728

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