New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 396792033 (0.630)6 (0.120) 6/238 NT 20.4% coding (4/1179 nt) PP_3498 transposase
?NC_002947 = 3967989 15 (0.300)coding (73/1179 nt) PP_3498 transposase

GACCGCACTGAGTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGAGCTTTGCACAACCATCTGCGTAGACCTGGCCAAACAGGTCTTCCAATTAGCAGGCGAGGACGCTACTGGTCGGGTGATCTACGAAGATCGC  >  NC_002947/3967842‑3968132
                                                                                                                                                   |                                                                                                                                               
gACCGCACTGAGTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGAcc                                                                                                                                               >  1:837179‑M2/1‑148 (MQ=255)
     cACTGAGTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACGCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGaa                                                                                                                                          >  1:1234093/1‑150 (MQ=255)
      aCTGAGTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGaaa                                                                                                                                         >  1:1413024/1‑150 (MQ=255)
        tGAGTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGaaaaa                                                                                                                                       >  2:933035/1‑150 (MQ=255)
           gTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCa                                                                                                                                    <  1:982951/150‑1 (MQ=255)
               gCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGt                                                                                                                                >  1:792465/1‑150 (MQ=255)
                        ggCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAacac                                                                                                                       <  2:1057289/150‑1 (MQ=255)
                          ccTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGACGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATa                                                                                                                                  <  1:1154127/137‑1 (MQ=255)
                          ccTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATa                                                                                                                                  >  2:1154127/1‑137 (MQ=255)
                               tGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTccc                                                                                                                <  1:933035/150‑1 (MQ=255)
                                      tAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTc                                                                                                         <  2:831496/150‑1 (MQ=255)
                                             aTGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGAt                                                                                                  <  1:239113/150‑1 (MQ=255)
                                                  tgcccACCCCGATTATGTGCAATCTATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGa                                                                                             <  2:792465/147‑1 (MQ=11)
                                                    tCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCATCAACTGATGTGGAGc                                                                                           >  2:813160/1‑150 (MQ=255)
                                                         cccGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGAGCTTTGc                                                                                      >  2:1322158/1‑150 (MQ=38)
                                                                 ggTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGAGCTTTGCACaa                                                                                  >  2:1011418/1‑146 (MQ=31)
                                                                 ggTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGAGCTTTGCACaa                                                                                  <  1:1011418/146‑1 (MQ=31)
                                                                  gtgtAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGAGCTTTGCACAACCATc                                                                             >  1:686551/1‑150 (MQ=11)
                                                                          gATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGAccatcgattacacctgatcggggtggacacacatcggtttagcggccaagaggccctgcccggcgaactcagtgcg                                                                     <  2:1331532‑M2/150‑77 (MQ=255)
                                                                             ggACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGAccatcgattacacctgatcggggtggacacacataggtttagcggccaagaggccctacccggcgaactcagtgcggtc                                                                  <  1:783299‑M2/150‑80 (MQ=255)
                                                                               aCTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGAccatcgattacacctgatcggggtggacacacatcggtttagcggccaagaggccctgcccggcgaactcagtgcggtccg                                                                <  2:292499‑M2/150‑82 (MQ=255)
                                                                                                           cACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGAccatcgattacacctgatcggggtggacacacatcggtttagcggccaagaggccctgcccggcgaactcagtgcggtccgcgtgggatggtcttgagcaagtcttcgg                                    >  1:1331532‑M2/1‑41 (MQ=255)
                                                                                                                        gaagaaGCCAACGGCAGCGAGGGCGAGTccatcgattacacctgatcggggtggacacacatcggtttagcggccaagaggccctgcccggcgaactcagtgcggtccgcgtgggatggtcttgagcaagtcttcggggctgatatgccc                       <  2:214227‑M2/150‑123 (MQ=255)
                                                                                                                                             ggCGAGAccatcgattacacctgatcggggtggacacacataggtttagcggccaagaggccctgcccggcgaactcagtgcggtccgcgtgggatggtcttgagcaagtcttcggggctgatatgccccaccacctgggccaccgcgct  >  2:783299‑M2/1‑7 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GACCGCACTGAGTTCGCCGGGCAGGGCCTCTTGGCCGCTAAACCGATGTGTGTCCACCCCGATCAGGTGTAATCGATGGACTCGCCCCCGCCGAGGCATCACAGTGCCACGGTGGTGTTCGAAGAAGCCAACGGCAGCGAGGGCGAGACCATGAAAAAGCATAGTTCGAAACACGATTCCCTGTCTTCCACTGATGTGGAGCTTTGCACAACCATCTGCGTAGACCTGGCCAAACAGGTCTTCCAATTAGCAGGCGAGGACGCTACTGGTCGGGTGATCTACGAAGATCGC  >  NC_002947/3967842‑3968132

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