New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 729706NA (NA)5 (0.100) 5/150 3.2 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.

GCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCA  >  NZ_CP009273/729665‑729812
                                                                |                                                                                   
gCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaaaa                                                            >  1:216751/1‑90 (MQ=33)
      tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                      <  1:1576483/90‑1 (MQ=35)
      tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                      <  2:58554/90‑1 (MQ=33)
           aTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATc                                                 <  1:1343171/90‑1 (MQ=35)
             gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                                               >  2:1385956/1‑90 (MQ=33)
             gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                                               >  1:117546/1‑90 (MQ=33)
             gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                                               >  2:313480/1‑90 (MQ=33)
              tgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCg                                              <  1:1230001/90‑1 (MQ=35)
               ggCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCgg                                             >  2:1104382/1‑90 (MQ=33)
               ggCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCgg                                             >  2:831011/1‑90 (MQ=33)
                 ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                                           <  1:1104382/90‑1 (MQ=35)
                           gTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc                                 <  1:1389342/90‑1 (MQ=33)
                             aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                            >  1:1288587/1‑77 (MQ=34)
                             aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                                            <  2:1288587/77‑1 (MQ=35)
                             aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGc                               <  1:1048218/90‑1 (MQ=33)
                              tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                      <  2:405464/66‑1 (MQ=35)
                              tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                                      >  1:405464/1‑66 (MQ=35)
                               cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                             <  2:768260/90‑1 (MQ=33)
                               cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                             <  2:95361/90‑1 (MQ=33)
                               cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                             <  1:981081/90‑1 (MQ=33)
                               cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                             <  1:1503964/90‑1 (MQ=33)
                                  ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGga                          >  1:669584/1‑90 (MQ=33)
                                    gTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCACCCCGGAAGGGGTTATCCGCACGgatg                        >  1:1360792/1‑90 (MQ=18)
                                          cgcACTATCTTTACCGGTATGACCGTCCCGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgc                  <  1:511590/90‑1 (MQ=18)
                                              cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac              <  2:1293227/90‑1 (MQ=33)
                                                    tacccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCg        <  2:113421/87‑1 (MQ=18)
                                                      aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                             <  1:806289/67‑1 (MQ=35)
                                                      aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc                             >  2:806289/1‑67 (MQ=35)
                                                      aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGt      >  2:1586362/1‑90 (MQ=35)
                                                      aCCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGt      <  2:216751/90‑1 (MQ=33)
                                                       ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa     <  2:1268992/90‑1 (MQ=33)
                                                       ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa     <  2:170646/90‑1 (MQ=33)
                                                       ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa     <  1:1586362/90‑1 (MQ=33)
                                                         ggTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAcc   <  1:1253114/90‑1 (MQ=33)
                                                          gTATGACCGTCACGGCGGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGGTATCCGCACGGGGGATGAGCGCACCCCCCGGGACCa  >  1:565891/1‑90 (MQ=11)
                                                          gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac              <  1:278802/78‑1 (MQ=34)
                                                          gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac              >  2:278802/1‑78 (MQ=35)
                                                          gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa  >  1:1111883/1‑90 (MQ=33)
                                                          gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa  >  1:441975/1‑90 (MQ=33)
                                                          gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa  >  2:1096273/1‑90 (MQ=35)
                                                            aTGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa     <  1:158735/85‑1 (MQ=33)
                                                            aTGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa     >  2:158735/1‑85 (MQ=35)
                                                                |                                                                                   
GCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCA  >  NZ_CP009273/729665‑729812

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