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.

CCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGC  >  NC_000913_3_bme_pgi/4150294‑4150569                                                                                                                                            |                                                                                                                                          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                                    <  2:686629/149‑1 (MQ=255)                                                                                              gAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttctt                                    <  1:286210/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)                                                                                                                  aCCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCgg                <  2:224348/149‑1 (MQ=255)                                                                                                                      tCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTgca            <  1:82893/149‑1 (MQ=255)                                                                                                                      tCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTgca            <  1:347631/149‑1 (MQ=255)                                                                                                                          gcAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCaa        >  1:421166/1‑149 (MQ=255)                                                                                                                                gCACGAGGGTTCGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGc  >  2:848735/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                          CCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACG‑CATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGC  >  NC_000913_3_bme_pgi/4150294‑4150569

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.