New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260480 251597 =5 (0.180)4 (0.160) 4/196 2.0 28.0% coding (497/1023 nt) h16_B0225 Integrase/recombinase protein
?AM260480 251666 = 16 (0.640)coding (566/1023 nt) h16_B0225 Integrase/recombinase protein

CGGCCGGCATGCCGGCCGGCGACACGCGCACCGCGGCATCGCCCACCCAATGGAAGCGCGCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCGTGCCGCTGTTCAGCCTGGCCCATGCGCCGCTACGGGCCTGGCTG  >  AM260480/251540‑251793
                                                                                                                              |                                                                                                                               
ggaccatgcgcagcgccatccggccgggggcgccggcttcgatggcatcgacccgcagcgcccgcgcttcgcctaccttcaggccgccgccgagcagcaccgccagcagggccacatcgcgggcgcTGCGGGTc                                                                                                                          <  1:932569‑M2/8‑1 (MQ=255)
                                     aTCGCCCACCCAATGGAAGCGCGCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGgcgc                                                                                     >  1:1925/1‑134 (MQ=255)
                                                 gacccgcaGCGCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTcc                                                                         <  1:932569‑M1/126‑1 (MQ=255)
                                                           GCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  2:1925‑M1/1‑134 (MQ=255)
                                                           GCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  1:239180‑M1/1‑134 (MQ=255)
                                                           GCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  2:762788‑M1/1‑134 (MQ=255)
                                                           GCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  1:616905‑M1/1‑134 (MQ=255)
                                                              CGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAAc                                                              <  1:110812‑M1/132‑1 (MQ=255)
                                                              CGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAAc                                                              >  2:110812‑M1/1‑132 (MQ=255)
                                                              CGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACgg                                                            <  2:98481‑M1/134‑1 (MQ=255)
                                                                cgcttcgcctaccttcaggccgccgccgagcagcaccgccagcagggccacatcgcgggcgcTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGcc                                                          >  2:537310‑M2/63‑134 (MQ=255)
                                                                CGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGcc                                                          <  1:872416‑M1/134‑1 (MQ=255)
                                                                    gtgGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  2:66786/1‑125 (MQ=255)
                                                                    gtgGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               <  1:66786/125‑1 (MQ=255)
                                                                    gtgGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               <  2:58904/125‑1 (MQ=255)
                                                                    gtgGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  1:58904/1‑125 (MQ=255)
                                                                     cgcctaccttcaggccgccgccgagcagcaccgccagcagggccacatcgcgggcgcTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAg                                                       <  1:537310‑M2/75‑1 (MQ=255)
                                                                             ctgGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCg                                              <  1:147830/133‑1 (MQ=255)
                                                                               ggCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTc                                                                          >  1:896792/1‑103 (MQ=255)
                                                                               ggCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTc                                                                          <  2:896792/103‑1 (MQ=255)
                                                                                        gctCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               <  2:176798/105‑1 (MQ=255)
                                                                                        gctCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  1:176798/1‑105 (MQ=255)
                                                                                         ccgagcagcaccgccagcagggccacatcgcgggcgcTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCGTGCCGCTGTTCAg                                 >  2:467688‑M2/38‑134 (MQ=255)
                                                                                           cggcggcggCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  2:540114/1‑102 (MQ=255)
                                                                                           cggcggcggCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               <  1:540114/102‑1 (MQ=255)
                                                                                                   gcCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCGTGCCGCTGTTCAGCCTGGCCCAt                       <  1:98481/134‑1 (MQ=255)
                                                                                                        aaGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCGTGCCGCTGTTCAGCCTGGCCCATGcgcc                  <  2:932569/134‑1 (MQ=255)
                                                                                                                  aaGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  2:623449/1‑79 (MQ=255)
                                                                                                                  aaGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               <  1:623449/79‑1 (MQ=255)
                                                                                                                        gggcgcTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCGTGCCGCTGTTCAGCCTGGCCCATGCGCCGCTACGGGCctggctg  <  2:872416‑M2/128‑1 (MQ=255)
                                                                                                                         ggcgcTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               <  1:670406‑M2/67‑1 (MQ=255)
                                                                                                                         ggcgcTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACaa                                                               >  2:670406‑M2/6‑72 (MQ=255)
                                                                                                                              |                                                                                                                               
CGGCCGGCATGCCGGCCGGCGACACGCGCACCGCGGCATCGCCCACCCAATGGAAGCGCGCCCGCGATGTGGCCCTGCTGGCGGTGCTGCTCGGCGGCGGCCTGAAGGTAGGCGAAGCGCGGGCGCTGCGGGTCGATGCCATCGAAGCCGGCGCCCCCGGCCGGATGGCGCTGCGCATGGTCCGCGCGGACAACGGCCGAGCCTATACCGTGCCGCTGTTCAGCCTGGCCCATGCGCCGCTACGGGCCTGGCTG  >  AM260480/251540‑251793

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