New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1207805	11 (0.280)	24 (0.690)	22/248	0.7	79.0%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_tpiA	= 1209602	3 (0.090)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_pae_tpiA/1209454‑1209748                                                                                                                                                     |                                                                                                                                                   cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                                    >  2:257062/1‑149 (MQ=255)            ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGACGGCCACGGACCCGGGACCCCACCAGGTAATGCAAAGccttctcccaa                                                                                                                                         <  2:82217‑M2/149‑12 (MQ=255)              cAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                                                                                       >  2:341067‑M2/1‑136 (MQ=255)                 ttcatttGTCAGATAAGCACTGGGCAGATCCGTGTATGCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTCAGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCaa                                                                                                                                    <  2:127170/146‑1 (MQ=255)                   tatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                                                                                                  >  1:257062‑M2/1‑131 (MQ=255)                   tatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACAGGGACTGTTGCAGGGAATGCAGAGccttctcccaaaccaagg                                                                                                                                  >  1:168730‑M2/1‑131 (MQ=255)                     tttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt                                                                                                                                >  2:632938‑M2/1‑129 (MQ=255)                          cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                           >  2:153464‑M2/1‑124 (MQ=255)                            gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                         >  1:432469‑M2/1‑122 (MQ=255)                                      tgacAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCATGCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACTAACAGGTAATGCAGAGCATTCTCCCAAACCAACGTTTATGAAAATgaagatgt                                                                                                               <  2:437684/146‑4 (MQ=255)                                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacag                                                                                                            <  1:531482‑M2/149‑41 (MQ=255)                                               cGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAAACGGTAATGCAGAGccttctcccaaacccaggttttcgaaagccgtttccacagtcccat                                                                                                      <  2:284564‑M2/149‑47 (MQ=255)                                                  gTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg                                                                                                   <  1:441736‑M2/149‑50 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcccatttga                                                                                             <  2:226785‑M2/149‑56 (MQ=255)                                                             cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg                                                                                        <  1:498149‑M2/149‑61 (MQ=255)                                                                        gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatatgatatcgccaaccgggtg                                                                             <  1:222797‑M2/149‑72 (MQ=255)                                                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttg                                                                         >  2:531482‑M2/1‑74 (MQ=255)                                                                             aTTGGATTTCAGCCAGCGTGGTGGCGGAGTGGCTGAAGTCCACGGAATCGGGACACCAACAGTTAATGCAGAGccttctcctaaaccaaggttttcgagatccgtttgcacagtgccatccgatttgatatcgccaaacggtttcttgc                                                                        >  2:430229‑M2/1‑73 (MQ=255)                                                                                         ccagccTGTGTGCGTAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                                            >  2:101515‑M2/1‑61 (MQ=255)                                                                                          cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcag                                                           <  1:433802‑M2/149‑90 (MQ=255)                                                                                                  tGGCGGAGTTGCTTAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                   <  1:42563/149‑1 (MQ=255)                                                                                                         gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                            <  1:509032‑M2/149‑105 (MQ=255)                                                                                                                      aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg                                                                                                          <  1:609712‑M2/74‑43 (MQ=255)                                                                                                                      aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg                                                                                                          >  2:609712‑M2/1‑32 (MQ=255)                                                                                                                            ccGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgc                         >  2:509032‑M2/1‑26 (MQ=255)                                                                                                                              gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGTTTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                       >  1:1244/1‑149 (MQ=255)                                                                                                                                    ccAACAGGTAATGCAGAGccttctcccaacccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctc                 >  1:282947‑M2/1‑18 (MQ=255)                                                                                                                                      aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca               <  1:684852‑M2/149‑134 (MQ=255)                                                                                                                                                gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcacc     <  2:460293‑M2/149‑144 (MQ=255)                                                                                                                                                   gagccttgtcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggccgtgccgccccttctccaggctggcaccgga  >  2:220616‑M2/1‑3 (MQ=255)                                                                                                                                                   gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccgga  >  2:384581‑M2/1‑3 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_pae_tpiA/1209454‑1209748

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.