New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	3344992 =	40 (0.840)	3 (0.070)	3/262	6.8	7.2%	coding (276/1434 nt)	rpoN	RNA polymerase, sigma 54 (sigma N) factor
	?	NC_000913_3_bme_pgi	3345024 =	40 (0.910)					coding (308/1434 nt)	rpoN	RNA polymerase, sigma 54 (sigma N) factor
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGCTGCTTGAGCAAATCGACACTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTA  >  NC_000913_3_bme_pgi/3344850‑3345140                                                                                                                                               |                                                                                                                                                     cGCTGCTTGAGCAAATCGACACTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCt                                                                                                                                                <  2:702052/149‑1 (MQ=255)      ctTGAGCAAATCGACACTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAGGAGCTGCCGCTCGGTGACAGTTGGGAAACCATTTACACCGCTGGTac                                                                                                                                           >  1:424601/1‑149 (MQ=255)       ttGAGCAAATCGACACTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTaca                                                                                                                                          <  1:602023/149‑1 (MQ=255)                tCGACACTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCACCCCGc                                                                                                                                 <  2:432973/149‑1 (MQ=255)                     acTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAg                                                                                                                            >  2:818972/1‑149 (MQ=255)                       gaaggaagaaTTCGACACCCGCGAAAGGCAAGACAGTGAAACGCTGGACACCGCCGAAGAGCACGAAAAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCg                                                                                                                          <  2:360165/145‑1 (MQ=255)                         aggaagaaATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGt                                                                                                                        <  2:502380/147‑1 (MQ=255)                         acgaagaaATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGt                                                                                                                        <  1:349410/147‑1 (MQ=255)                         aTGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGt                                                                                                                        <  1:150059/149‑1 (MQ=255)                         aTGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGt                                                                                                                        <  2:689868/149‑1 (MQ=255)                                  accacacCCGCCAAACGCAAGACAGTCAAACCCTGGACACCGCCCACGCGCTCGAACAAAAAGAGATGCCGGAACAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCACTAGTACACCACCCGGCACCAGCGGTGACTACAtt                                                                                                               <  2:6759/146‑1 (MQ=255)                                       acCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTgacga                                                                                                          >  1:95673/1‑149 (MQ=255)                                         ccGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTAAACCGCTGGTCAACCATCCAGCCCCAGCGGTGATTAAATTAacgacg                                                                                                        >  2:420608/1‑149 (MQ=255)                                               aaCGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGc                                                                                                  <  2:95673/149‑1 (MQ=255)                                                      gACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTa                                                                                           >  2:131304/1‑149 (MQ=255)                                                         aGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAg                                                                                                      <  1:824757/135‑1 (MQ=255)                                                         aGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAg                                                                                                      >  2:824757/1‑135 (MQ=255)                                                           tGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAgg                                                                                      <  1:748908/149‑1 (MQ=255)                                                            gAAACGGTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAggg                                                                                     <  2:109078/149‑1 (MQ=255)                                                               aCGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGATGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGAATAGATTGACGCCGAGCTGCCGGTCTACCAGGGCGa                                                                                  >  1:463608/1‑149 (MQ=255)                                                                   tGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAacg                                                                              >  2:122712/1‑149 (MQ=255)                                                                   tGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGCCTACAAGGGCGAAacg                                                                              >  2:121669/1‑149 (MQ=255)                                                                          cgccgACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGa                                                                       >  2:232844/1‑149 (MQ=255)                                                                                        aaCAAAAGGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGCACCACATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTAc                                                         <  2:382227/149‑1 (MQ=255)                                                                                         aCAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTAcc                                                        >  2:31728/1‑149 (MQ=255)                                                                                           aaaaaGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACAt                                                                                                                >  2:387213/1‑91 (MQ=255)                                                                                           aaaaaGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACAt                                                                                                                <  1:387213/91‑1 (MQ=255)                                                                                           aaaaaGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCCAAACGACGCCGCCTCTGCAAGATTACCTg                                                      >  2:404028/1‑149 (MQ=255)                                                                                             aaaGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGCCGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGAt                                                    >  1:645193/1‑149 (MQ=255)                                                                                                          aaGAGCTGCCGGTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAg                                       <  2:345717/149‑1 (MQ=255)                                                                                                           agagCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAGACGACGCAGACCTTGCAGGATTACTCGATGTGGCAGATTGAGc                                      >  2:562040/1‑149 (MQ=255)                                                                                                                    gCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGCATTACCTGATGTGGCAGGTTGAGCTGACACCGt                             >  2:789397/1‑149 (MQ=255)                                                                                                                        gATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGGCACCGTTTTc                         >  2:324831/1‑149 (MQ=255)                                                                                                                             cAGTTGGGACACCATTTCCACCGCTGGGAACCCATCCGGGACCAGCGGTGACCACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATTTGGCAGGTTGAGCTGACACCGTTTTCCGaca                    <  2:53932/149‑1 (MQ=255)                                                                                                                               gTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGCGCTGACACCGTTTTCCGACACt                  >  2:272036/1‑149 (MQ=255)                                                                                                                                   ggACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGAcc              >  1:165771/1‑149 (MQ=255)                                                                                                                                   ggACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGAcc              <  1:31728/149‑1 (MQ=255)                                                                                                                                      cacCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGGCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACcgcg           <  2:129127/149‑1 (MQ=255)                                                                                                                                       acCATTGACCCCGCTGGGACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCGTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCgcgc          <  1:79204/149‑1 (MQ=255)                                                                                                                                         cATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACATTGACCGCGCTa        >  2:159900/1‑149 (MQ=255)                                                                                                                                               cacCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTa  <  1:232844/149‑1 (MQ=255)                                                                                                                                               cacCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGAAGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCCGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGGTa  >  1:41566/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     CGCTGCTTGAGCAAATCGACACTCATGAAGAAATCGACACCCGCGAAACGCAAGACAGTGAAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGAAACGACGCAGACCTTGCAGGATTACCTGATGTGGCAGGTTGAGCTGACACCGTTTTCCGACACTGACCGCGCTATTGCTA  >  NC_000913_3_bme_pgi/3344850‑3345140

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