New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1207805	62 (1.380)	5 (0.130)	5/254	5.2	10.2%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_bme_pgi	= 1209602	33 (0.830)					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.

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAT  >  NC_000913_3_bme_pgi/1209454‑1209732                                                                                                                                                     |                                                                                                                                   cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTCAAGGCCCTGGCACCGTGACACCAACAGGTAATGCcgag                                                                                                                                    >  2:363391/1‑149 (MQ=255)    gTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGCAGCTGAAGGCCACGGAACCGGGACACCAACAGGTGATGCAGAGCCt                                                                                                                                 >  1:397246/1‑149 (MQ=255)       aaTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtct                                                                                                                              >  1:21152/1‑149 (MQ=255)            ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTTAGCTGAAAAAGCGGCACCATTGCATGTCAGCCAGCCTGTTGACGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGCGCCTTCTCCCaa                                                                                                                         >  2:754888/1‑149 (MQ=255)             gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGAACCCAACAGGTAATGAAGAACCCTCTCCCaac                                                                                                                        >  2:246471/1‑148 (MQ=255)             gCAATTTACTTGTCCGGAATGCTTTTGTCAGTTCCGTGTATTCTTCCGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGCCTGAACTCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCgaa                                                                                                                        <  2:203156/149‑1 (MQ=255)                  ttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                                                                   <  2:21152/149‑1 (MQ=255)                      ttGTCGGATAAGCCTTTGACAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACCGTTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                               <  2:640686/149‑1 (MQ=255)                                aGCCTTCGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                                                     <  1:184570/149‑1 (MQ=255)                                  ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGACGGAGTGGCTGAAGGCCACGGAACCGGGCCACCAACAGGTAATGCAGAGCCTTCTCCCCAACCAACGTTTACGAAAATgaac                                                                                                   >  2:180045/1‑148 (MQ=255)                                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAGTGAACAAATAAc                                                                                            >  1:323036/1‑149 (MQ=255)                                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                                                            <  1:501814/149‑1 (MQ=255)                                           gTTCCGGAGATTCTTCAGCAGAAAAAGCCGCACCATTGCATTTCAGCCAGCCTGTTGGCTCGCTGGCGGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                          <  1:183651/149‑1 (MQ=255)                                                gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGaaaa                                                                                     >  1:46321/1‑149 (MQ=255)                                                     ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGGAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                <  2:531126/149‑1 (MQ=255)                                                     ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                >  1:129998/1‑149 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                             >  1:561625/1‑149 (MQ=255)                                                        ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCATGCCTTTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCACACAAACGTTTATGAAAATGAAGCAATACCAAGCAACTGGcatca                                                                             >  2:369285/1‑149 (MQ=255)                                                         tcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                            >  1:318681/1‑149 (MQ=255)                                                          cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGCCTTACGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                           <  1:182522/149‑1 (MQ=255)                                                           agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAAAAGGAAATGCAGAACCTTCTTCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                          <  2:367336/149‑1 (MQ=255)                                                             cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                        >  1:534759/1‑149 (MQ=255)                                                                            cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                         >  2:172980/1‑149 (MQ=255)                                                                             aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGTAACCCACGGCACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTCTGAAAATGAAGAAATAAAAAGCAAATGGCATCCTTCCTCCTTTTTCCAggggac                                                        >  2:563891/1‑147 (MQ=255)                                                                              ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                       <  2:776638/149‑1 (MQ=255)                                                                                  aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAGGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                            <  2:358062/140‑1 (MQ=255)                                                                                  aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAGGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                            >  1:358062/1‑140 (MQ=255)                                                                                  aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGCCCTCATTCCTGCGTTTACCAGGCGGATTTaa                                                   >  2:84889/1‑149 (MQ=255)                                                                                        gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTCTGAAAACGGAGAAATAACAACCAAATAGCATCATTCCTGCATTTACCAGGGGTATCTAACATGCt                                             >  2:282187/1‑149 (MQ=255)                                                                                          cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACACGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGGTTTAACATGCTTa                                           >  1:386050/1‑149 (MQ=255)                                                                                                            gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                         >  1:25249/1‑149 (MQ=255)                                                                                                              tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                             <  1:755302/143‑1 (MQ=255)                                                                                                              tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                             >  2:755302/1‑143 (MQ=255)                                                                                                                 aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca                    <  1:596099/149‑1 (MQ=255)                                                                                                                    ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                 <  1:779982/149‑1 (MQ=255)                                                                                                                      aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa               <  1:447291/149‑1 (MQ=255)                                                                                                                      aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa               <  1:172980/149‑1 (MQ=255)                                                                                                                           aCCGGGACAACAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa          <  2:306511/149‑1 (MQ=255)                                                                                                                                   acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt  <  1:255097/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                   CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAT  >  NC_000913_3_bme_pgi/1209454‑1209732

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.