New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 729706NA (NA)12 (0.110) 10/154 3.5 NA pseudogene (103/1506 nt) BW25113_RS03635 RHS repeat protein
?NZ_CP009273 = 729729 NA (NA)pseudogene (126/1506 nt) BW25113_RS03635 RHS repeat protein
Rejected: Coverage evenness skew score above cutoff.

CCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGA  >  NZ_CP009273/729658‑729818
                                                                       |                                                                                         
ccGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACg                                                                         <  2:1898884/90‑1 (MQ=31)
 cGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACgg                                                                        >  2:1944394/1‑90 (MQ=32)
 cGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACgg                                                                        >  1:123501/1‑90 (MQ=32)
 cGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACgg                                                                        >  2:2799434/1‑90 (MQ=32)
     cAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaa                                                                    <  1:2779538/90‑1 (MQ=35)
      aGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaaa                                                                   >  2:1271086/1‑90 (MQ=33)
      aGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaaa                                                                   >  1:1913237/1‑90 (MQ=33)
         aCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTTAGGGAGaaaaac                                                                >  2:2444836/1‑88 (MQ=255)
             tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                            >  1:2407065/1‑90 (MQ=33)
             tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                            <  1:164586/90‑1 (MQ=35)
             tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                            <  2:1695663/90‑1 (MQ=35)
               aGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTc                                                          <  1:1626294/90‑1 (MQ=35)
                 cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                                        >  1:3291646/1‑90 (MQ=33)
                 cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                                        >  2:531917/1‑90 (MQ=33)
                 cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                                        >  1:2957603/1‑90 (MQ=33)
                  aTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATc                                                       <  2:3112950/90‑1 (MQ=35)
                   tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGTAGGCTGACGGAGAAAACCGACCTCATcc                                                      <  2:65815/90‑1 (MQ=25)
                   tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc                                                      <  2:2829458/90‑1 (MQ=35)
                    gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                                                     >  2:3110434/1‑90 (MQ=33)
                    gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                                                     >  2:3049523/1‑90 (MQ=33)
                    gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc                                                      >  2:688032/1‑89 (MQ=33)
                    gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc                                                      <  1:688032/89‑1 (MQ=35)
                       gCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                                  <  1:2649147/90‑1 (MQ=35)
                        ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                                                 <  2:1364704/90‑1 (MQ=35)
                        ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                                                 >  1:2684829/1‑90 (MQ=33)
                        ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                                                 <  2:391178/90‑1 (MQ=35)
                         cGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAg                                                <  1:2937041/90‑1 (MQ=35)
                          ggATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAACCGACCTCATCCCGGAAggg                                              >  1:743520/1‑90 (MQ=16)
                          ggATAACCGTATCGACCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAgg                                               >  1:511965/1‑90 (MQ=18)
                              aaCCGTATCGCCCGTGACGGGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGtt                                           <  1:1942814/90‑1 (MQ=25)
                                  gTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc                                       <  1:2300812/90‑1 (MQ=33)
                                  gTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc                                       <  1:3316686/90‑1 (MQ=33)
                                    aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGc                                     >  2:940561/1‑90 (MQ=33)
                                     tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGGTGACGGAGAAAACCCACCGCTAACCGGTAGGGTTTATTCGCa                                    >  2:1578830/1‑90 (MQ=11)
                                      cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   <  2:3193559/90‑1 (MQ=33)
                                      cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   <  1:1815014/90‑1 (MQ=33)
                                      cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   <  2:2913645/90‑1 (MQ=33)
                                      cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   <  1:1205962/90‑1 (MQ=33)
                                      cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   <  2:3124899/90‑1 (MQ=33)
                                       gCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACg                                  <  2:2968755/90‑1 (MQ=33)
                                       gCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACg                                  <  1:1201941/90‑1 (MQ=33)
                                         ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGGGGAAACCGACCTCATCCCGGAAGGGGTTATCCACACggg                                >  2:193398/1‑89 (MQ=11)
                                         ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGga                                >  1:2322676/1‑90 (MQ=33)
                                          cGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgat                               >  1:3002712/1‑90 (MQ=33)
                                                gcgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGGGAAAACCGACCTCATCCCGGAAGGGGTTATCCCCCCGGATgaggag                         >  2:451786/1‑90 (MQ=11)
                                                gcgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAg                         >  2:1607857/1‑90 (MQ=33)
                                                 cgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGGTATAAGCACGGGTGATGAgc                        >  2:2354897/1‑90 (MQ=11)
                                                  gcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcg                       >  1:1082361/1‑90 (MQ=33)
                                                     cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    <  2:123501/90‑1 (MQ=33)
                                                       aTCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcaccc                  <  1:1117509/90‑1 (MQ=33)
                                                       aTCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcaccc                  <  2:1913237/90‑1 (MQ=33)
                                                          tttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacccacc               <  2:1308978/90‑1 (MQ=33)
                                                          tttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacccacc               <  2:2407065/90‑1 (MQ=33)
                                                           ttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCg              >  2:1070336/1‑90 (MQ=33)
                                                            tACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGGCCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCCCCgg             >  2:2708027/1‑90 (MQ=255)
                                                             aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    <  2:3180863/82‑1 (MQ=33)
                                                             aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                    >  1:3180863/1‑82 (MQ=35)
                                                             aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGt            >  1:961580/1‑90 (MQ=35)
                                                              ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  2:2731551/90‑1 (MQ=33)
                                                              ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  1:1093642/90‑1 (MQ=33)
                                                              ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  1:3176252/90‑1 (MQ=33)
                                                              ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  2:2684829/90‑1 (MQ=33)
                                                               cGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  2:2572670/89‑1 (MQ=33)
                                                               cGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           >  1:2572670/1‑89 (MQ=35)
                                                               cGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAc          >  1:1878991/1‑90 (MQ=35)
                                                                ggTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAcc         <  1:457018/90‑1 (MQ=33)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTc                                                          <  2:2264274/40‑1 (MQ=255)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTc                                                          >  1:2264274/1‑40 (MQ=40)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCTGCACGGATGATGAGCGCACCCCCCGGTACCa        >  2:814375/1‑90 (MQ=255)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   >  2:2646878/1‑63 (MQ=35)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                                   <  1:2646878/63‑1 (MQ=35)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcaccc                  >  2:2719462/1‑80 (MQ=34)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcaccc                  <  1:2719462/80‑1 (MQ=34)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        <  1:2943931/90‑1 (MQ=33)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:2073514/1‑90 (MQ=33)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:3380682/1‑90 (MQ=33)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:857482/1‑90 (MQ=35)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:903965/1‑90 (MQ=33)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:435852/1‑90 (MQ=33)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:1640384/1‑90 (MQ=35)
                                                                 gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGGTATCCGCACGGCTGATGAGCGCACCCACCGGTACCa        >  2:471044/1‑90 (MQ=255)
                                                                 gTATGACCGTCACGGCAGGATGACGGAGAAAACCGACCCCATCCCGGAAGGGGTTATCCTCACGGATGCTGAGCACACCCCCCGGGACCa        >  2:1246334/1‑90 (MQ=11)
                                                                  tATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAt       <  1:3108403/90‑1 (MQ=33)
                                                                  tATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAt       <  2:2259668/90‑1 (MQ=33)
                                                                   aTGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAtt      >  1:401533/1‑90 (MQ=35)
                                                                   aTGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAtt      >  2:522988/1‑90 (MQ=33)
                                                                   aTGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAtt      >  2:2175170/1‑90 (MQ=33)
                                                                     gACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc    >  2:577671/1‑90 (MQ=35)
                                                                     gACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc    <  1:1070336/90‑1 (MQ=33)
                                                                       ccGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa  <  1:1486525/90‑1 (MQ=33)
                                                                       ccGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa  <  1:940561/90‑1 (MQ=33)
                                                                       ccGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa  <  2:1078115/90‑1 (MQ=33)
                                                                       ccGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa  <  1:1820039/90‑1 (MQ=33)
                                                                       |                                                                                         
CCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGA  >  NZ_CP009273/729658‑729818

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