New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 206661324 (0.400)4 (0.070) 4/284 6.9 14.9% coding (709/1017 nt) insH1 IS5 transposase and trans‑activator
?NC_000913_3_bme_pgi = 2102140 NA (NA)coding (772/1017 nt) insH1 IS5 transposase and trans‑activator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCA  >  NC_000913_3_bme_pgi/2066465‑2066761
                                                                                                                                                    |                                                                                                                                                    
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:569519/150‑1 (MQ=32)
                                                                                                 cGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGc                                                    >  1:348196/1‑150 (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)
                                                                                                                                                gcCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGGGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGCCTTCATg     <  2:511594/150‑1 (MQ=21)
                                                                                                                                                   ccTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCa  >  1:258012/1‑150 (MQ=34)
                                                                                                                                                    |                                                                                                                                                    
TTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGTGGATGCTGTTTCAAGGTTCTTACCTTGCCGGGGCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCCCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACTAGGCTGTGGGTCAGGCCGCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCA  >  NC_000913_3_bme_pgi/2066465‑2066761

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.