New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 2066557	23 (0.380)	3 (0.050)	3/284	7.5	12.0%	coding (765/1017 nt)	insH1	IS5 transposase and trans‑activator
	?	NC_000913_3_bme_pgi	= 2066606	22 (0.380)					coding (716/1017 nt)	insH1	IS5 transposase and trans‑activator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGCCTTC  >  NC_000913_3_bme_pgi/2066457‑2066755                                                                                                                                                      |                                                                                                                                                      gCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCtgtg                                                                                                                                                       >  1:154416/1‑150 (MQ=32)    ccACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGgcg                                                                                                                                                    >  2:726782/1‑150 (MQ=32)         ttGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCt                                                                                                                                               <  2:892421/150‑1 (MQ=11)         ttGGCCCGGATGCGGGCTTTCATGTATTCGATGTTGATGGCCGGTTTGGTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCt                                                                                                                                               <  2:922177/150‑1 (MQ=2)                                                                                                       cGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTa                                                 <  1:449110/150‑1 (MQ=31)                                                                                                          cGATCAGCCAGTCCACCTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCAGCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGc                                              <  1:566937/150‑1 (MQ=17)                                                                                                          cGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGc                                              >  1:348196/1‑150 (MQ=32)                                                                                                          cGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGc                                              <  1:569519/150‑1 (MQ=32)                                                                                                            aTCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCtg                                            >  1:477748/1‑150 (MQ=32)                                                                                                               aGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGCGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCAGGCTCGATGGCCGCGGTGGTGACTAGGCTGTgg                                         >  1:26446/1‑150 (MQ=12)                                                                                                                 ccAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTCCCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGt                                       >  1:556363/1‑150 (MQ=17)                                                                                                                   aggccACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCa                                     <  1:679494/147‑1 (MQ=17)                                                                                                                    gTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAg                                    >  1:488143/1‑150 (MQ=32)                                                                                                                      ccACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCCGCTGAGTGAGGTCATGCTCGTTGGCCGCGGTGGGGACTAGGCTGTGGGTCAGGc                                  >  1:562656/1‑150 (MQ=12)                                                                                                                         caTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGc                               >  2:836093/1‑150 (MQ=33)                                                                                                                              aCCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCACGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTg                          >  1:716252/1‑150 (MQ=21)                                                                                                                               ccTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTgg                         <  1:284122/150‑1 (MQ=34)                                                                                                                                    gCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGATCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCATGCCGCTCATGGCATCg                    >  1:737576/1‑150 (MQ=255)                                                                                                                                           cctcGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCaa             >  1:134373/1‑150 (MQ=34)                                                                                                                                            ctcGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAACTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAAt            >  2:677982/1‑150 (MQ=21)                                                                                                                                              cgcgCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGCGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGCTTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACGAAtgt          >  2:173745/1‑150 (MQ=255)                                                                                                                                                cgcTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTgg        >  2:867501/1‑150 (MQ=34)                                                                                                                                                   tgtgGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGcc     <  1:187263/150‑1 (MQ=34)                                                                                                                                                      ggCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCCTGCAGCAGATTCCCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGCCACCAATGTGGGCCTTc  >  1:673653/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGCCTTC  >  NC_000913_3_bme_pgi/2066457‑2066755

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