New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	4280509 =	30 (0.840)	5 (0.150)	5/268	4.1	15.8%	coding (529/1650 nt)	yjcE	putative cation/proton antiporter
	?	NC_000913_3_bme_pgi	4280529 =	25 (0.730)					coding (549/1650 nt)	yjcE	putative cation/proton antiporter
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTGGCGCTCTCCGGGATTGTTGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGCAGATCGTCCTGCT  >  NC_000913_3_bme_pgi/4280364‑4280654                                                                                                                                                  |                                                                                                                                                  tgtgGCGCTCTCCGGGATTGTTGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCg                                                                                                                                                >  1:388969/1‑149 (MQ=255) tgtgGCGCTCTCCGGGATTGTTGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCg                                                                                                                                                >  1:389007/1‑149 (MQ=255)           tcCGGGATTGTTGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGc                                                                                                                                      >  1:77326/1‑149 (MQ=255)             cGGGATTGTTGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTACAGTTTGCCGTAGCCGTGGCGCTGGGGACGATGATCTTCACCGTCGGCGGTGCGa                                                                                                                                    >  2:173867/1‑149 (MQ=255)             agggATTGAAGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGa                                                                                                                                    >  2:6584/2‑149 (MQ=255)                     ttgGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGaa                                                                                                                            <  2:332177/149‑1 (MQ=255)                         cGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTa                                                                                                                        <  2:176211/149‑1 (MQ=255)                               gcgcATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATATTCACCGTCGGCGGTGCGACGGTCGAATTTAGGaagg                                                                                                                  >  2:312881/1‑147 (MQ=255)                                cgcATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGt                                                                                                                 <  1:237784/149‑1 (MQ=255)                                       cGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCAtt                                                                                                          >  1:79779/1‑149 (MQ=255)                                       cGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGAGGAACTTCACTGTCGTTGGTGCGACGGTCGAATTTATGAAAGTAGCCAtt                                                                                                          >  2:611128/1‑149 (MQ=255)                                                     ggCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGc                                                                                            >  1:208571/1‑149 (MQ=255)                                                                 ggCCAAGCGTTGATGAACAACGCCCCCGGCCTGGACTCCTTGAAGTTTGCCGGAGCAGTGGCGATGGGGACGAAGAGCCTCACCCGCGTCGGTGCAACGCTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCCGTTTTgtggt                                                                                <  2:273608/149‑1 (MQ=255)                                                                        cGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTgtggtgggatga                                                                         >  2:55702/1‑143 (MQ=255)                                                                           tgatgaATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAgctggctg                                                                      >  1:20190/1‑149 (MQ=255)                                                                              tgaatgaCGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTAc                                                                   <  1:285393/149‑1 (MQ=255)                                                                                  tgaCGGTTCCGGGTTGGGGTCGTTGAAGTGTGCGGTAGCAGTGGCGATGGGGATGATGATCTTGACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTTGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGtc                                                               <  2:374361/149‑1 (MQ=255)                                                                                             gTCTGGTGTCGTTGACGTTTGCCGTAGCAGTGGCGATGGGGAAGCTGATATTCATCTTCTGCTTATCGACGGGCGAATTCATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGa                                                    <  2:314573/149‑1 (MQ=255)                                                                                                  gtgtCGTCGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGTAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTATGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCt                                               >  1:589000/1‑149 (MQ=255)                                                                                                                 gCCGTAGCAGTGGCGAGGGGGACGCGGATCTTCGCCGTGGGCGGTGTGACGGTCGAATTTATGAAAGTGGCCATTGGCGGTATTCTCGCCGGGTTTGTGGTGGGCGGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCgg                                <  2:399286/149‑1 (MQ=255)                                                                                                                  ccGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGGGCGACGGTCGAATGTATGAATGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGt                               <  2:233073/149‑1 (MQ=255)                                                                                                                    gTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGtga                             <  1:6584/149‑1 (MQ=255)                                                                                                                                aTGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGAc                 >  1:400317/1‑149 (MQ=255)                                                                                                                                  ggggACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGc               <  1:59001/149‑1 (MQ=255)                                                                                                                                     gACGATGATCTTCACCGTCGGTGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGCAGa            <  1:81081/149‑1 (MQ=255)                                                                                                                                       cGATTATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTTGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGCAGATc          <  1:70733/149‑1 (MQ=255)                                                                                                                                        gatgatCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGCAGATCg         <  2:79779/149‑1 (MQ=255)                                                                                                                                               ttCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGtc                                                               >  2:605794/1‑88 (MQ=255)                                                                                                                                               ttCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGtc                                                               <  1:605794/88‑1 (MQ=255)                                                                                                                                               ttCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGCAGACCGTCctgct  >  2:591601/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TGTGGCGCTCTCCGGGATTGTTGGCGAAGGGCGCATCCCGAAAAAAATCATGGGCATTTTGCAGGGCGAAGCGTTGATGAATGACGCCTCCGGTCTGGTGTCGTTGAAGTTTGCCGTAGCAGTGGCGATGGGGACGATGATCTTCACCGTCGGCGGTGCGACGGTCGAATTTATGAAAGTAGCCATTGGCGGTATTCTCGCCGGTTTTGTGGTGAGCTGGCTGTACGGTCGCTCGCTGCGATTCCTCAGCCGCTGGGGCGGTGATGAACCCGCGACGCAGATCGTCCTGCT  >  NC_000913_3_bme_pgi/4280364‑4280654

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.