New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1207805	48 (0.580)	8 (0.110)	8/252	5.4	14.9%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_tpiA	= 1209602	49 (0.670)					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.

ATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_pae_tpiA/1207659‑1207952                                                                                                                                                   |                                                                                                                                                    aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                                                                   <  1:864181/149‑1 (MQ=255)  tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                                  <  1:875185/149‑1 (MQ=255)           aGCCTTGAGAAGGGGCGGCACAACCAGTTACTTACCGCACTACGGGCGCTGCTGTTAAGTTGTCAGAATCCGTATGGAGATATCAAATCCGATCGCACTGTGCAACCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaacgcttcg                                                                                                                                         <  2:1262253/149‑8 (MQ=17)               cggggaAGGGGCGGCAGTACCAGTTACTTGCCGCGCTTCCCGCGCTGCTGTTAATCCGCAAGAATTCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                     <  2:659919/145‑1 (MQ=18)                       gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGGGCAAAACTCAATGCAGCAACGGc                                                                                                                             >  1:87858/1‑149 (MQ=255)                           ggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGACAACCGGTGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                                                                                                         >  1:1350702/1‑149 (MQ=25)                            gcAAGACCAGTTACTTACCGCCCTTCGAGCGCTGCCGAAAACAGGGAGGAATTCTGCCGGGGATTGCTGGTCGGATTGCGCGGTGCTATCGGCTCTCGAAAGCTTGGGATTGGGAGAGTGAGCAAAACTCAATGCAGCAACGGCAACAt                                                                                                                        >  2:1369126/1‑149 (MQ=255)                                 aCCAGTTACTTTCCGCACTTCGCGCGCTGCTGTTAGGCCGCAAGAATCCGTTGGGCGATATCAACTCGAAAGCCCACGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                                                                                                   <  1:562600/149‑1 (MQ=21)                                     gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCACAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                                                                               >  2:675183/1‑149 (MQ=255)                                               gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAAGACGCACCGGTTTCAt                                                                                                     >  1:943232/1‑149 (MQ=255)                                                   ttCGGGCGCTTCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                                                                 <  2:812001/149‑1 (MQ=255)                                                    tCGCGCGCTGCGGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                                                                <  1:602355/149‑1 (MQ=255)                                                       cgcgcTGCTGTTGAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                                             <  2:1325533/149‑1 (MQ=255)                                                            tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt                                                                                        >  2:127715/1‑149 (MQ=255)                                                              ctgTTAAGCCTCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                                                                                      <  2:843543/149‑1 (MQ=255)                                                                 ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCGTAGATATACCGGTTATGAt                                                                                   >  1:768778/1‑149 (MQ=255)                                                                    agtcGCAAGAAGCCGTTTGGCGATATCCAACCGGAGGGCACTGTGCAAACGGCTCTCCAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                                <  2:578849/146‑1 (MQ=37)                                                                    aGCCGCAAGAAGCTGTTTGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTTGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                                <  1:296522/149‑1 (MQ=255)                                                                      ccGCAAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                                              <  2:997069/149‑1 (MQ=255)                                                                       cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg                                                                             >  1:573322/1‑149 (MQ=255)                                                                           agagTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCACGAAAACCTTGGTTTGGGAGGAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCCATACCGGTTCTGATGGGAGGCAtt                                                                         >  2:1044935/1‑149 (MQ=37)                                                                           aGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                         <  1:127715/149‑1 (MQ=255)                                                                             gagccGTGTTGCGATATAAAAGGGGATGGCACTGTGCAAACGGCTCGCGAAAACCTTGGTTTGGGAGAAGGAGCAAACCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                                       <  2:1144979/146‑1 (MQ=38)                                                                              aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                                                      <  2:1323576/149‑1 (MQ=255)                                                                               tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGc                                                                     <  2:1086169/149‑1 (MQ=255)                                                                                  gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                                                                                       <  2:1168430/112‑1 (MQ=255)                                                                                  gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                                                                                       >  1:1168430/1‑112 (MQ=255)                                                                                  gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                  <  2:1041038/149‑1 (MQ=255)                                                                                  gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                  <  1:1211475/149‑1 (MQ=255)                                                                                  gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                  <  1:675183/149‑1 (MQ=255)                                                                                             tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAtt                                                       <  2:337125/149‑1 (MQ=255)                                                                                             tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAtt                                                       >  1:456900/1‑149 (MQ=255)                                                                                               aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCa                                                     <  1:1041244/149‑1 (MQ=255)                                                                                                    ggATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                <  1:383817/149‑1 (MQ=255)                                                                                                         gCACTTTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTgg                                           <  2:571931/149‑1 (MQ=255)                                                                                                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGaa                                         >  2:693645/1‑149 (MQ=255)                                                                                                             tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGGTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATg                                       >  1:914407/1‑149 (MQ=255)                                                                                                                     cGGCGCTCGAAAGCCTGGGTTTGGGAGGAGGAGCAAAACGCAATGCAGCGACGGCTACATTAGGACGCACCGGTTTCGTAGCTATAGCGGTTATGATTGGTTGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                               <  2:645349/149‑1 (MQ=255)                                                                                                                           tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCt                         >  2:1417568/1‑149 (MQ=255)                                                                                                                                aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctc                    >  2:139324/1‑149 (MQ=255)                                                                                                                                   cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAAt                                        >  2:572908/1‑126 (MQ=255)                                                                                                                                   cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAAt                                        <  1:572908/126‑1 (MQ=255)                                                                                                                                     tGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg               <  2:1350702/149‑1 (MQ=255)                                                                                                                                          tGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGAt          <  2:1171911/149‑1 (MQ=255)                                                                                                                                              agaAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaa      <  2:87858/149‑1 (MQ=255)                                                                                                                                                aaGGGGGAAAACTAAATGCAGCAACGGCTAAATTAGGAGGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATACCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt    <  2:716095/149‑1 (MQ=255)                                                                                                                                                aaGGAGCAAAACTCAATGCAGCAACGGCTACAATAGGACGCACCGGTTTCATAGCTAAACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAATGGCGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGAGGGATGGAAAt    >  2:1250701/1‑149 (MQ=255)                                                                                                                                                  ggAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAc  >  1:1370433/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_pae_tpiA/1207659‑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.