New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	1207790 =	27 (0.860)	5 (0.180)	4/252	3.6	18.5%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_tpiA	1209619 =	20 (0.720)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAG  >  NC_000913_3_pae_tpiA/1207644‑1207924                                                                                                                                                   |                                                                                                                                       cAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCTGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTg                                                                                                                                      <  1:362350/149‑1 (MQ=255)   gCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                                                                                    >  1:402216/1‑149 (MQ=255)    cGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCAAAAACCTTGGtt                                                                                                                                   >  2:84584/1‑149 (MQ=255)    cGTGGTGGAGGAATCCGGTGGCAGCCTGGAGGAGGGGCGGCGCGACTAGTTACTTACCGGACTTCGCGCGCTTCTGTGAAGCCGCAAGAATCCGTTGGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGtt                                                                                                                                   <  2:186952/149‑1 (MQ=18)        gtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCGCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGGTATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTggg                                                                                                                               >  1:127320/1‑149 (MQ=255)             ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAg                                                                                                                          <  2:148410/149‑1 (MQ=255)             ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAg                                                                                                                          <  1:400585/149‑1 (MQ=255)               cATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACCGACCGCACTTCGCGCGCTGCTGTTAAGCGGCAAGAATCCGTTGGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                        <  2:511442/149‑1 (MQ=255)                               ggAGACGGGGCGGCACGACCAGTTACCTACCGCACTTCTCGTGCTTCTGTTAAGCCGCAAGAATCCGTTTGGCGCGATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                                        <  1:97111/149‑1 (MQ=18)                                     ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTGTCGAAAACATTGGTTTGGGAGAAGGAGCAAAAGTCAATGCAGCGACgg                                                                                                  >  1:11213/1‑149 (MQ=255)                                        gcggcACGAGCAGTGACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                                                               <  2:383470/149‑1 (MQ=255)                                           gcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAAAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                                            >  2:266567/1‑149 (MQ=255)                                                  cAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAAAGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAAGAcc                                                                                     >  2:102698/1‑148 (MQ=255)                                                         ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                                                              >  2:538008/1‑149 (MQ=255)                                                               cACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                        >  1:135693/1‑149 (MQ=255)                                                                         gctgttGTTAAGCCGCAAGAAACCGTGTGGCCAAATCAAATCGGATGGCACGGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg                                                              <  1:375451/149‑1 (MQ=25)                                                                          ccgcggTTAAGCCGCAAGAATCCGGTTGGCGATATCAAATTGGATGGCACTGTGCAAGTGGCTCTCGAAATCCTTGGTTTGGGAGAACGAGCAATTCTCAATGCTGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATGCCGGt                                                             <  2:50126/144‑1 (MQ=255)                                                                                     ccGCAAGAACCCGTTTGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACACAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                  <  1:174379/149‑1 (MQ=255)                                                                                      cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCACAAGGCGCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGGTAACTTAGGAAGCACCGCTTTACTAGCTATACCGGTCATGATtggtgg                                                 >  2:241238/1‑149 (MQ=21)                                                                                       gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGGGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGAGCGCACCGGTTTCATAGCTATACCGGTTATGATTagtggt                                                >  1:48596/1‑149 (MQ=255)                                                                                          aGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                             >  1:62488/1‑149 (MQ=255)                                                                                             aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                          <  1:538008/149‑1 (MQ=255)                                                                                                        gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGCAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGAACGCACCGGTTTCATAGCTATACAGGTAATGCTTGGTGGTATTGAGCAATCAGTAAt                               >  2:509398/1‑149 (MQ=38)                                                                                                                             gtgAAAACGCCTCTCGAAAACCTTGGTGTGGGAGAAGGAGAAAAACTCAATGAAGCAACGGCTACATTAGGACGCACCCATTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGc          <  1:94561/149‑1 (MQ=255)                                                                                                                                     ggCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCacaaa  >  2:38316/1‑148 (MQ=255)                                                                                                                                          tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                                                          >  2:10378/1‑72 (MQ=255)                                                                                                                                          tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                                                          <  1:10378/72‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                       CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAG  >  NC_000913_3_pae_tpiA/1207644‑1207924

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