New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3968936 =10 (0.300)9 (0.300) 6/246 3.0 42.8% coding (1020/1179 nt) PP_3498 transposase
?NC_002947 3968992 = 15 (0.500)coding (1076/1179 nt) PP_3498 transposase
Rejected: Coverage evenness skew score above cutoff.

GATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGC  >  NC_002947/3968876‑3969128
                                                                                                                    |                                                                                                                                        
gatgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGTGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGc                                                                                                         >  1:515230/1‑150 (MQ=17)
       cAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATTGGCCCGGTTa                                                                                                  >  1:58961/1‑150 (MQ=17)
           gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGGAAGGAGGCGGCGATAAGGCTGGCTAATCGAAATGTGCGGATTATATGGGCCCTGTTACaga                                                                                              >  1:375955/1‑150 (MQ=12)
               gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                                                                              >  1:238638/1‑146 (MQ=32)
               gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                                                                              <  2:238638/146‑1 (MQ=32)
                 aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATc                                                                                        >  2:318038/1‑150 (MQ=32)
                             gCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGTGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                                                                            <  2:515230/150‑1 (MQ=17)
                                      ggCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcga                                                                   <  1:634129/150‑1 (MQ=32)
                                       gCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                                                                  <  2:568621/150‑1 (MQ=32)
                                                       ccAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcg                                                  >  2:760764/1‑150 (MQ=32)
                                                            tgatCGCCGCCTCTTGCGCTGGCTATCCCGGCTTGTGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTAt                                             <  2:375955/150‑1 (MQ=17)
                                                                  ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGCGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCct                                       >  2:108244/1‑150 (MQ=17)
                                                                      ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                   <  1:760764/150‑1 (MQ=32)
                                                                       tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGCCTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGAc                                  >  2:592427/1‑150 (MQ=17)
                                                                                                gcgcAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAGCCGAACGAAAATCCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgcgc         <  1:226383/150‑1 (MQ=17)
                                                                                                       gaggCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGgcac  <  2:792991/150‑3 (MQ=17)
                                                                                                          gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                      >  1:339269/1‑127 (MQ=33)
                                                                                                          gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                      <  2:339269/127‑1 (MQ=32)
                                                                                                                    |                                                                                                                                        
GATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGC  >  NC_002947/3968876‑3969128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ 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.