New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	4150405 =	49 (0.870)	7 (0.130)	6/260	6.5	13.4%	intergenic (‑140/+42)	eptC/ppc	LPS heptose I phosphoethanolamine transferase/phosphoenolpyruvate carboxylase
	?	NC_000913_3_bme_pgi	4150433 =	45 (0.850)					intergenic (‑168/+14)	eptC/ppc	LPS heptose I phosphoethanolamine transferase/phosphoenolpyruvate carboxylase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGC  >  NC_000913_3_bme_pgi/4150259‑4150553                                                                                                                                                   |                                                                                                                                                      aaTGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGaaaa                                                                                                                                                     <  1:588014/149‑1 (MQ=255)      aTAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTc                                                                                                                                                >  1:138268/1‑149 (MQ=255)      aTAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTc                                                                                                                                                <  2:688441/149‑1 (MQ=255)                 cTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCGTCGCGCAAAAGCa                                                                                                                                     >  1:855908/1‑149 (MQ=255)                 cTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCa                                                                                                                                     <  1:589092/149‑1 (MQ=255)                 cTGTCTTGAAAAAAGATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCa                                                                                                                                     <  2:448409/149‑1 (MQ=255)                   gTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACg                                                                                                                                   >  2:286210/1‑149 (MQ=255)                    tCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGa                                                                                                                                  >  1:82585/1‑149 (MQ=255)                     cTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAg                                                                                                                                 <  1:133078/149‑1 (MQ=255)                      ttGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAACACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAgg                                                                                                                                >  2:758788/1‑149 (MQ=255)                                    ccGAAAGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAgg                                                                                                                  <  1:918669/149‑1 (MQ=255)                                            aaCGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCCCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAg                                                                                                          >  2:611275/1‑149 (MQ=255)                                             aCGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGc                                                                                                         <  2:275556/149‑1 (MQ=255)                                               gACAATTCCTTAAGGATATCTGAAAGTACATTTTAAAGGTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCg                                                                                                       >  2:966072/1‑149 (MQ=255)                                               gACAATTCCTTAAGGATATCTGAAAACATGTTCTTTGTTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAAGGTTTGCAGAAGAGGAAAATTAGcgg                                                                                                       >  2:966116/1‑147 (MQ=255)                                                       cTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACg                                                                                               >  1:686629/1‑149 (MQ=255)                                                        ttAAGGATATCTGAAGGTATATTCAGCATTTGAATAAAATGAAGACAGAAATATATTTAAAAGGAGGGGGTAAGAACAGAAGTATTTCATAAAACCCTCGCGCAAACGCGAAGGAGTATCACGAAGAGGAAAACTAGCCGGTATTAcca                                                                                              >  2:52990/1‑147 (MQ=255)                                                        ttAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑c                                                                                             <  1:562322/149‑1 (MQ=255)                                                              atatCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCt                                                                                       >  2:836206/1‑149 (MQ=255)                                                                     aaGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCaa                                                                                >  2:994670/1‑149 (MQ=255)                                                                          atatTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAGGAGTAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCg                                                                           <  1:291763/149‑1 (MQ=255)                                                                           tatTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATAc                                                                                         <  1:272639/134‑1 (MQ=255)                                                                           tatTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATAc                                                                                         >  2:272639/1‑134 (MQ=255)                                                                             ttCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCa                                                                        >  2:396032/1‑149 (MQ=255)                                                                               cAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAATATTAGCCGGTATTACG‑CCTACCTGCCGCCATCCCGGCAAt                                                                      >  1:380684/1‑149 (MQ=255)                                                                                aGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATa                                                                     <  1:564611/149‑1 (MQ=255)                                                                                          tAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAACACCCTCGCGCAAAAGCACGGGGGTTTGCAGAAGCGGAAGATTAGCCGGTATTACG‑CATACCTCCAGCACTCCAGGCAATACTGACCAtga                                                           >  2:158940/1‑147 (MQ=255)                                                                                          tAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTa                                                           >  1:309082/1‑149 (MQ=255)                                                                                           aaaaTGCAGACAGAAATATATTGAAACCGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTaa                                                          <  1:759747/149‑1 (MQ=255)                                                                                            aaaGGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAAc                                                         <  2:408210/149‑1 (MQ=255)                                                                                              acgCAGACAGACACATAGTGAAAAGGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGc                                                       <  2:321090/147‑1 (MQ=255)                                                                                                gCAGACAGAAATATATGGACAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCtt                                                     <  1:370056/149‑1 (MQ=255)                                                                                                   gacagaAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGtt                                                  <  2:82585/149‑1 (MQ=255)                                                                                                               ttGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGAt                                      >  1:986024/1‑149 (MQ=255)                                                                                                                  aaaaCGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCg                                   >  1:983054/1‑149 (MQ=255)                                                                                                                   aaaCGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCCACACGT‑CCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCgg                                  >  1:734703/1‑149 (MQ=255)                                                                                                                        aGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCAGTAACGCTTGTTCGACGCGAGGATCCGGTTCCt                             >  1:950049/1‑149 (MQ=255)                                                                                                                             gtTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCct                        <  1:611275/149‑1 (MQ=255)                                                                                                                             gtTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCct                        >  1:594082/1‑149 (MQ=255)                                                                                                                                aGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCCCGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGACcttct                     <  2:138268/149‑1 (MQ=255)                                                                                                                                 gAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttctt                    <  1:286210/149‑1 (MQ=255)                                                                                                                                 gAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttctt                    <  2:686629/149‑1 (MQ=255)                                                                                                                                  aaCAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCgg                                  <  1:108442/134‑1 (MQ=255)                                                                                                                                  aaCAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCgg                                  >  2:108442/1‑134 (MQ=255)                                                                                                                                    cAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCttttt                 <  2:983054/149‑1 (MQ=255)                                                                                                                                       aaGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTg              >  1:964420/1‑149 (MQ=255)                                                                                                                                        aGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGc             <  2:919456/149‑1 (MQ=255)                                                                                                                                            tttCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTgg         >  1:512604/1‑149 (MQ=255)                                                                                                                                                  aaaaCCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAg   <  1:758788/149‑1 (MQ=255)                                                                                                                                                   aaaCCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGc  <  1:776985/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      AATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGC  >  NC_000913_3_bme_pgi/4150259‑4150553

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