New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1207805	6 (0.150)	25 (0.680)	21/254	0.9	85.7%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_tpiA	= 1209602	3 (0.080)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_pae_tpiA/1209466‑1209748                                                                                                                                         |                                                                                                                                                   gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                                                                            >  2:612938‑M2/1‑137 (MQ=255) gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                                                                            <  1:612938‑M2/145‑9 (MQ=255)            gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTTGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                             <  2:313093/149‑1 (MQ=255)                     gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAGAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                    <  1:658065‑M2/149‑33 (MQ=255)                     gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtgttcgagagccgt                                                                                                                    <  1:659400‑M2/149‑33 (MQ=255)                      ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                                                                                   <  2:289314‑M2/149‑34 (MQ=255)                      ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGAGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                                                                                   <  1:210163‑M2/149‑34 (MQ=255)                                ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagcgc                                                                                                         <  2:353543‑M2/149‑44 (MQ=255)                                   cGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccat                                                                                                      <  1:518110‑M2/149‑47 (MQ=255)                                         ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                                >  2:518110‑M2/1‑97 (MQ=255)                                          tcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                               >  1:45352‑M2/1‑96 (MQ=255)                                            ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttcgcaccgtgccatccgatttga                                                                                             <  1:305149‑M2/149‑56 (MQ=255)                                            ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacggtgccatccgatttga                                                                                             <  1:378681‑M2/149‑56 (MQ=255)                                                 cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg                                                                                        <  1:27664‑M2/149‑61 (MQ=255)                                                      aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                   <  2:462410/149‑1 (MQ=255)                                                                 aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                                        >  1:289314‑M2/1‑73 (MQ=255)                                                                  ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccgtccgatttgatatcgccgaacgggtttttgcg                                                                       <  1:63618‑M2/149‑78 (MQ=255)                                                                    gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc                                                                     >  2:27664‑M2/1‑70 (MQ=255)                                                                     cATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccccaaccaaggttttcgagggccgtttgcacagtgccatccgatttgatatcgccacacggattcttgcgggt                                                                    <  1:248823‑M2/149‑81 (MQ=255)                                                                      aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctt                                                                   <  1:688020‑M2/149‑82 (MQ=255)                                                                             ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                                            >  1:148357‑M2/1‑61 (MQ=255)                                                                                  cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                                       >  2:217254‑M2/1‑56 (MQ=255)                                                                                                    aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAAAAATAACAAGCAAATGGCATCACCCGGTCTCTCTCCAGTGGGATTTAACATGCTTGTTGGCTATGTCCGCGTATcacc                                     >  1:20903/1‑147 (MQ=255)                                                                                                           cGGAACCTCGACACCAACCGGTAATTCAAAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggttgaacagcagcgcgcgaagtgcggtaagtaactggtcgt                              >  2:518883‑M2/1‑31 (MQ=255)                                                                                                                ccGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgc                         >  1:292704‑M2/1‑26 (MQ=255)                                                                                                                  gggACACCAACAGGTAATGCAGAGcattatccaaatcccaggttttcgaaacccgtttacacaatgccatcccattttagatcgccaaactgattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                       >  2:63618‑M2/1‑24 (MQ=255)                                                                                                                               gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgacgtgcggtacgtaactggtcgtgccgccccttctcccggctg          <  1:334486‑M2/149‑139 (MQ=255)                                                                                                                               gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg          <  1:216511‑M2/149‑139 (MQ=255)                                                                                                                                    gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcacc     <  1:530503‑M2/149‑144 (MQ=255)                                                                                                                                       tggcctttttccagaccatggttttcgagagccgtttgcactgtgccatccgatttgattgtgccatacggattcttgcggcttaacagcagtgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccgga  >  2:321701‑M2/3‑3 (MQ=255)                                                                                                                                         |                                                                                                                                                   GCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_pae_tpiA/1209466‑1209748

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.