New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1207805	33 (0.850)	7 (0.200)	6/254	4.0	18.6%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_tpiA	= 1209602	32 (0.930)					pseudogene (18/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.

GCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_pae_tpiA/1207657‑1207952                                                                                                                                                     |                                                                                                                                                    gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTCCCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACGGGCCAACCGGCCCTCGACAACCTTGGTTTGGGCGAAgg                                                                                                                                                     >  2:91007/1‑149 (MQ=17)    tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                                  >  2:439667/1‑149 (MQ=255)      cGGTGGCAGCCTGGACAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAACCCGTAAGAATCCGTTTGGCGATATTAAATCCGATGCCACTGTGTAAACGGCTCTCGAAAATCTTGGCTTGGGATAAggctctt                                                                                                                                                <  2:527088/149‑6 (MQ=17)                 tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                     >  1:88364/1‑149 (MQ=255)                 tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                     >  1:547460/1‑149 (MQ=255)                        ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGGGAAGGAGCAAAAGTCAATGCAGCAAcag                                                                                                                              >  1:315166/1‑147 (MQ=255)                        ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGGAAAGGACCAAAACTCAATGCAGCAACgg                                                                                                                              >  2:38218/1‑149 (MQ=255)                        ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCACTGCAGCAACgg                                                                                                                              >  2:399809/1‑149 (MQ=255)                         gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCAGCAAGAATCCGTTTAGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCCATGCAGCAACGGc                                                                                                                             >  1:600525/1‑149 (MQ=255)                                 cGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                                                                                                                     >  2:47123/1‑149 (MQ=255)                                        ttacttacCGCACTTCGCGCGCTGCTGTGTTTCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                                                                                              <  2:315166/149‑1 (MQ=255)                                              acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                                        >  1:206522/1‑149 (MQ=255)                                                cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                                                                                      >  2:400169/1‑149 (MQ=255)                                                    ctgcgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCCATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                                                                  <  2:336723/146‑1 (MQ=255)                                                                ctgTTAAGCCGCAAGAATCCGTTTGGCGATAGCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGTACCGGTTTCATAGCTATACCGGTTAt                                                                                      <  2:349849/149‑1 (MQ=255)                                                                   ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                                                                                   >  1:373062/1‑149 (MQ=255)                                                                      aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGACACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                                <  1:489796/149‑1 (MQ=255)                                                                       gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAGCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                               <  1:671462/149‑1 (MQ=255)                                                                       gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                               <  2:442495/149‑1 (MQ=255)                                                                             acgagccGTTTGGCGATATAAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGTTTTGGGAGAAGGAGCAAAACTCAATGCATCAACGGCTACATTAGGACGCCCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                         <  2:545345/144‑1 (MQ=255)                                                                                          cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTAGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTGTTGAGCAATCAGTaa                                                            >  2:46779/1‑149 (MQ=255)                                                                                            atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc                                                          <  1:47123/149‑1 (MQ=255)                                                                                                 aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCa                                                     >  2:367602/1‑149 (MQ=255)                                                                                                      ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                <  1:564566/149‑1 (MQ=255)                                                                                                       gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                               <  1:367602/149‑1 (MQ=255)                                                                                                        aTGGCAGTGTGCAAACGGCTCTCGAAAACCTTCGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggggg                                              <  1:38952/149‑1 (MQ=255)                                                                                                         tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGt                                             <  2:637935/149‑1 (MQ=255)                                                                                                         tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGt                                             <  1:524563/149‑1 (MQ=255)                                                                                                                             tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCt                         <  2:466082/149‑1 (MQ=255)                                                                                                                               gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGAAAAAGCATCTGc                       >  1:255383/1‑149 (MQ=255)                                                                                                                                     cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                 >  1:109236/1‑149 (MQ=255)                                                                                                                                      atGGTTCGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAACGAAGACGCACCGGCTACATAGCTATAAAGGTCATGAAGGGAGGGAGTGAGAAATCAGCAATCATAGGGAGGGGGTGAAACTCCGCAAAGACATCTACCCCTggg                >  2:330793/2‑149 (MQ=255)                                                                                                                                        ggTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg              <  2:554302/149‑1 (MQ=255)                                                                                                                                             gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATg         <  1:400169/149‑1 (MQ=255)                                                                                                                                                    ggAGCAAAGCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATGGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAc  <  1:619826/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_pae_tpiA/1207657‑1207952

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.