New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 447605335 (0.600)10 (0.180) 9/276 5.7 23.1% coding (243/321 nt) PP_3966 transposase
?NC_002947 = 4476057 NA (NA)coding (239/321 nt) PP_3966 transposase
Rejected: Coverage evenness skew score above cutoff.

AAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTGGCATTGATGCCGTGGTTGATG  >  NC_002947/4475986‑4476206
                                                                       |                                                                                                                                                     
atgcgcGAATAAATTGCGCACAACCCTCGGGATCCCAGGCTTGCCATTTCACCGTGATCACCTGCCCGGTTACCGGCAGCTGGATGATCACTGACGCTTCGACTGGCCGTTTAGGGGCGGCTGTTCTCCGGACGAGCGCTGGTAGTGAAg                                                                         <  2:310324/145‑1 (MQ=1)
            aTTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGCCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCa                                                             <  1:695302/150‑1 (MQ=11)
             ttGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGACCACATCGCCAGTGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGAGGAAAGCTGGTAGTGCAGCCGCTGCGGGCAg                                                            >  1:83286/1‑150 (MQ=11)
             ttGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCAATACTCGCTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAg                                                            <  1:90204/150‑1 (MQ=9)
                 gcACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGAtct                                                        <  2:193291/150‑1 (MQ=32)
                   acaACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCg                                                      >  2:330635/1‑150 (MQ=32)
                   acaACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCg                                                      >  1:6323/1‑150 (MQ=32)
                    caACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGa                                                     <  1:1219718/150‑1 (MQ=32)
                    caACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGa                                                     <  2:1086788/150‑1 (MQ=32)
                      aCCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATa                                                   <  2:572689/150‑1 (MQ=32)
                            gggATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCgg                                             >  1:541052/1‑150 (MQ=32)
                              gATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCNACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCa                                           >  1:904108/1‑150 (MQ=17)
                                tCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGc                                         <  2:855315/150‑1 (MQ=32)
                                 ccGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGcc                                        >  2:603897/1‑150 (MQ=32)
                                   gAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGGTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCAt                                      >  1:50437/1‑150 (MQ=17)
                                   gAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGGGGCAGCCAt                                      >  1:618747/1‑150 (MQ=17)
                                   gAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCAt                                      >  1:74799/1‑150 (MQ=32)
                                      gCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCTCGTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTg                                   <  1:428424/150‑1 (MQ=12)
                                       cTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGGGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGc                                  <  2:211290/150‑1 (MQ=17)
                                          gCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGa                               >  1:834730/1‑150 (MQ=32)
                                          gCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACCCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGa                               >  2:273449/1‑150 (MQ=17)
                                           ccATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGAt                              <  1:1191149/150‑1 (MQ=32)
                                           ccATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTAGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTTCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGAt                              <  1:1086480/150‑1 (MQ=17)
                                            cATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATg                             >  1:479456/1‑150 (MQ=32)
                                               ttCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACg                          <  2:311153/150‑1 (MQ=32)
                                                 cACTGTGCTCACTTGCCCGGCCACGAGGAACTCGATGATCACTGACGCTTAGACTGGCCGTTTAGGTTCGGCTTTTAGGGGGACGAAAGCGGGTAGCGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGtt                        <  2:70205/150‑1 (MQ=11)
                                                  aCCGTAACCACTTGCCGCGCCACGGGCAACGCGTTGATACCGGACGCTTGGCCTGGCCGTGTAGGTGCGGCTTTTAGCGGGACGAAAGCCGGTAGTCAAGCCGCTGCGGGCAGATCTCGATAGACCGGCAGCCATTTGCGGATGACGTTg                       <  2:242404/150‑1 (MQ=11)
                                                   ccGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTgg                      <  2:904108/150‑1 (MQ=32)
                                                      tGATCACTTGCCCGGCAACGGGCAACTCGAGGAGCACTGACGCTTAGAGTGGCCGTTTAGGTGGGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTGGCGGATGACGTTGCCac                   <  2:197187/150‑2 (MQ=11)
                                                       gATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTGGCAtt                  <  1:15618/150‑1 (MQ=32)
                                                               gCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCACATCTCGAAAGAGCGGCAGCCTTTTGCGGAGTGCGTTTGCTTTGATGCCGt          >  2:13033/1‑150 (MQ=11)
                                                                 ccGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTGGCATTGATGCCGTgg        <  1:1141587/150‑1 (MQ=32)
                                                                  cGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGACAGCCATTCGCGGATGACGTTGGCATTGATGCCGTGGt       >  2:341567/1‑150 (MQ=17)
                                                                   ggCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTGGCATTGATGCCGTGGtt      <  2:1215827/150‑1 (MQ=32)
                                                                       aCGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGGTTAGGTGCGGCTTTTAGCGTGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTGGCATTGATGCCGTGGTTGATg  <  1:489307/150‑1 (MQ=17)
                                                                       |                                                                                                                                                     
AAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGCCACGGGCAACTCGATGATCACTGACGCTTCGACTGGCCGTTTAGGTGCGGCTTTTAGCGGGACGAAAGCTGGTAGTGAAGCCGCTGCGGGCAGATCTCGATAGAGCGGCAGCCATTTGCGGATGACGTTGGCATTGATGCCGTGGTTGATG  >  NC_002947/4475986‑4476206

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