New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 206655723 (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.

TGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCT  >  NC_000913_3_bme_pgi/2066412‑2066614
                                                                                                                                                 |                                                         
tGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACATTGCCGTGTCGCTCgg                                                       >  2:942816/1‑150 (MQ=12)
 ggCTTTCACGAAGCCGAACTTTCGCTTGATCATGCGAAATGCGTGCTCCACCTTGGCCCGGAAGCAGGCTTTCATGTATTCGATGATGATGCCCGTTTTGTTCTTGCGTGGATGCTGGTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGc                                                      <  2:766350/150‑1 (MQ=11)
   cTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGa                                                    >  2:142810/1‑150 (MQ=32)
   cTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGa                                                    >  1:775163/1‑150 (MQ=32)
         cGAATCCGAACTGTCGGCTGATGATGCGAAATGGGCGCTCCTCCTTGGCCCGGATGCTGGCTTTCATGTATTCGTTCTTGCTGGCCGTCTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGcc                                              <  2:405974/150‑1 (MQ=11)
         cGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGAAATCGATGTTGATGGCCGTATTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGTGATCAGcc                                              <  2:868609/150‑1 (MQ=255)
                        gCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCt                               >  1:106172/1‑150 (MQ=32)
                           tgatgaTGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGTTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCgg                            >  2:679494/1‑150 (MQ=17)
                              tgatgCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCa                         >  1:782045/1‑150 (MQ=32)
                                   cGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGctcc                    <  1:812690/150‑1 (MQ=32)
                                   cGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGctcc                    >  2:600574/1‑150 (MQ=32)
                                    gAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGctcct                   >  1:503733/1‑150 (MQ=32)
                                    gAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGctcct                   >  1:477642/1‑150 (MQ=32)
                                      aaTGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTcg                 <  1:382220/150‑1 (MQ=32)
                                      aaTGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTcg                 <  1:378936/150‑1 (MQ=32)
                                          ggTGCTCGACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCt             <  1:714984/150‑1 (MQ=17)
                                          ggTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCt             <  1:749377/150‑1 (MQ=32)
                                          ggTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTGTTGTTCTTGCGGGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCt             <  1:830988/150‑1 (MQ=17)
                                            tGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCtgt           <  1:451129/150‑1 (MQ=32)
                                             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)
                                                                                                                                                 |                                                         
TGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCTTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCT  >  NC_000913_3_bme_pgi/2066412‑2066614

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.