New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 1207805	19 (0.810)	7 (0.330)	7/254	1.6	26.0%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_pae_pgi	= 1209602	23 (1.100)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_pae_pgi/1207661‑1207952                                                                                                                                                 |                                                                                                                                                    ccGGTCCCAGCCTGGAGAACGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                                                                                                 <  2:244377‑M1/149‑5 (MQ=255)    gTGCCAGCCTGGAGAA‑GGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                                                                                                             >  2:286075/1‑149 (MQ=255)     tGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                                                                                                             <  1:369676/149‑1 (MQ=255)                            aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                                      <  1:61440/149‑1 (MQ=255)                              gACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                                                                                                    >  1:353218/1‑149 (MQ=255)                              catttgttgatacgcgtacagagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                                                                                                    >  2:182376‑M1/101‑149 (MQ=255)                                 cAGTTACTTACCGGACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                                                                                                                 <  2:190355/149‑1 (MQ=255)                                    ttacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCCATGCAGCCACGGCTACATTAGGACGccc                                                                                                              >  2:10670/1‑147 (MQ=255)                                        tacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGCGTTAAAATTGTTGCAGCAACGGCTACATTAGGACGCACCGGt                                                                                                          >  1:99512‑M1/91‑149 (MQ=255)                                          acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                                        <  2:60137/149‑1 (MQ=255)                                          acCGCACTTCGCGCGATGCTGTTAAGCCGCAAGAACCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                                        <  1:94974/149‑1 (MQ=255)                                               acatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCATCGTCTACATTAGGACGCACCGGTTGCATAg                                                                                                   >  1:400262‑M1/84‑149 (MQ=255)                                                          tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt                                                                                        >  1:89416/1‑149 (MQ=255)                                                            ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                                                                                      <  2:57417/149‑1 (MQ=255)                                                                tAAGCCGCAAGAATCCGTTTGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                                                  <  1:158337/149‑1 (MQ=255)                                                                   gccgccATAATCTGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGCCTCTCGAAAACCTTGGGTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc                                                                               <  2:150353‑M1/149‑71 (MQ=255)                                                                   gccgccAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                               <  2:109917/149‑1 (MQ=255)                                                                      gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCCTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaat                                                                            >  1:178826‑M1/1‑76 (MQ=255)                                                                         ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                         <  2:191962‑M1/92‑1 (MQ=255)                                                                             tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgg                                                                          <  1:369429‑M1/144‑76 (MQ=255)                                                                             tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgg                                                                          >  2:369429‑M1/1‑69 (MQ=255)                                                                               cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                                                                   >  2:57166/1‑149 (MQ=255)                                                                                aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                  >  1:174023‑M1/51‑149 (MQ=255)                                                                                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                                       <  2:211983‑M1/149‑95 (MQ=255)                                                                                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAtt                                                       <  1:222631/149‑1 (MQ=255)                                                                                                  ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtaacttgcagctactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacc                                                <  1:183623‑M1/149‑102 (MQ=255)                                                                                                   gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                               <  2:353218/149‑1 (MQ=255)                                                                                                    aTGGCAATGTGCTAACGGCTCTCAAAAAACCTGGTTTGTGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggggg                                              <  2:92283/149‑1 (MQ=255)                                                                                                             tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTCGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc                                     >  2:325784/1‑149 (MQ=255)                                                                                                                              aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctc                    <  2:248522/149‑1 (MQ=255)                                                                                                                              aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctc                    >  2:76326/1‑149 (MQ=255)                                                                                                                                    ggTTTGGGAGAAGGctctgcagtacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaatt              <  2:259983‑M1/149‑136 (MQ=255)                                                                                                                                      tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaacaggcaaaggcttatccgacaaataaattgc            >  1:241238‑M1/1‑12 (MQ=255)                                                                                                                                           gggaagGctctgcattacctggtggggtcccggtgccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgat       <  2:161954‑M1/147‑143 (MQ=255)                                                                                                                                              aaGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgattta    >  1:381190‑M1/1‑4 (MQ=255)                                                                                                                                              AAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt    <  1:238910‑M1/149‑1 (MQ=255)                                                                                                                                                ggctctgcattacctgttggtgtccctgttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacg  <  1:166603‑M1/149‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATAC  >  NC_000913_3_pae_pgi/1207661‑1207952

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.