New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 288457460 (0.810)5 (0.070) 5/276 10.1 8.1% coding (2644/2667 nt) ygcB Cascade complex anti‑viral R‑loop helicase‑annealase Cas3
?NC_000913 = 2884578 57 (0.820)coding (2640/2667 nt) ygcB Cascade complex anti‑viral R‑loop helicase‑annealase Cas3
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGA  >  NC_000913/2884430‑2884726
                                                                                                                                                    |                                                                                                                                                    
tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                      <  1:1121613/149‑1 (MQ=255)
tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                      <  2:252742/149‑1 (MQ=255)
tgAATAGAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                      <  2:334270/149‑1 (MQ=255)
tgAAGATAAAATCGCATAACTGGAAGTTTAAAGAAGATGTACATTGTGCGCCTTCCCTACTTAAGTAGGGATAAACCGTTCTTGGTGTTATTATCGTCATTGATAACAATCATTCCCTGAGTTAGTTGGGATTTTCTGGGATGACTCTg                                                                                                                                                      <  2:234532/149‑1 (MQ=255)
  aaTATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGt                                                                                                                                                    >  1:1119004/1‑149 (MQ=255)
       aaaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTAATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATcc                                                                                                                                               >  1:20185/1‑149 (MQ=255)
            ctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCa                                                                                                                                          >  2:1123793/1‑149 (MQ=255)
                tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATccc                                                                                                                                      <  1:207383/149‑1 (MQ=255)
                  ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCt                                                                                                                                    <  1:1056490/149‑1 (MQ=255)
                    gggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGt                                                                                                                                  >  1:174329/1‑149 (MQ=255)
                        aaTTTAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGAAAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAg                                                                                                                              <  2:126518/149‑1 (MQ=255)
                          ttAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTGTTATTATCGTCATTGATAACAATCGTTCCCGAAGTTATTTGGGATTTGCATGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGt                                                                                                                            <  1:806611/149‑1 (MQ=255)
                           tAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACAGTTATTGGTCTTCTTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCTTCCCTTCATCCCCTGTATAGGta                                                                                                                           <  2:853840/149‑1 (MQ=255)
                            aaaaGAAGATGTACATGTGGCACCTTCCCGACTTAAGTAGGGATAAACCGTTACTGGTCTTATTATCGTCATTGATAACAATCAGTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaa                                                                                                                          <  1:1096830/149‑1 (MQ=255)
                             aaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaat                                                                                                                         <  2:662888/149‑1 (MQ=255)
                               agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaataa                                                                                                                       <  1:339217/149‑1 (MQ=255)
                                gaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAAc                                                                                                                      <  1:56288/149‑1 (MQ=255)
                                gaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCGTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAAc                                                                                                                      >  2:797502/1‑149 (MQ=255)
                                  agaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACaa                                                                                                                    <  2:403024/149‑1 (MQ=255)
                                                 aCCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGGATAGGTAATAACAATACTGTTACCCTGcc                                                                                                     >  1:996263/1‑149 (MQ=255)
                                                  ccTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACGACATTGGTACCCTGCCa                                                                                                    <  2:767541/149‑1 (MQ=255)
                                                               aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGCCCCatccatcc                                                                                       >  1:448713/1‑149 (MQ=255)
                                                                   gggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAt                                                                                   <  1:298987/149‑1 (MQ=255)
                                                                    ggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAtt                                                                                  >  1:1191792/1‑149 (MQ=255)
                                                                     gATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTCTAGGTAATAACAATACTGTTACTCTGCCAGACCCATCCATCCAGAtgc                                                                                 >  2:289178/1‑147 (MQ=255)
                                                                        aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGt                                                                              <  1:721320/149‑1 (MQ=255)
                                                                            cGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACCGTTGTCTGAGCCTGCCAGACCCATCCATCCAGATTCTGTTTcc                                                                          <  2:1136768/149‑1 (MQ=255)
                                                                                 ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTcc                                                                     >  2:408272/1‑149 (MQ=255)
                                                                                   ggTCCTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAg                                                                   <  1:770382/149‑1 (MQ=255)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGt                                                                              >  1:8325/1‑137 (MQ=255)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGt                                                                              <  2:8325/137‑1 (MQ=255)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCCGCATCCCCTGGATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCAATCCAAATTCGGTTTCCCTTCCAGc                                                                  >  1:199973/1‑149 (MQ=255)
                                                                                     tCTTATTATCGTCATTGATAACACTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATGGGTAATAACAATCCTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCTCTTCCCGcc                                                                 >  2:856409/1‑149 (MQ=255)
                                                                                     tCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGcc                                                                 >  2:988040/1‑149 (MQ=255)
                                                                                        tattatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCaaa                                                              >  2:560571/1‑149 (MQ=255)
                                                                                         attatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACCGTTACCCTGCCACACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAg                                                             >  2:1089051/1‑149 (MQ=255)
                                                                                            atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAACGTaa                                                          >  2:584544/1‑149 (MQ=255)
                                                                                            atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                          <  1:337171/149‑1 (MQ=255)
                                                                                               gTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCCGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAacc                                                       >  2:22690/1‑149 (MQ=255)
                                                                                                tCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                      >  1:721306/1‑149 (MQ=255)
                                                                                                tCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                      <  1:1126048/149‑1 (MQ=255)
                                                                                                 cATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCAt                                                     >  1:126259/1‑149 (MQ=255)
                                                                                                     gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTc                                                 <  2:174329/149‑1 (MQ=255)
                                                                                                      aTAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                                <  2:1113328/149‑1 (MQ=255)
                                                                                                           aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                           <  2:95747/149‑1 (MQ=255)
                                                                                                            aTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctact                                          <  2:975092/149‑1 (MQ=255)
                                                                                                                  cccGACGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                    <  2:764143/149‑1 (MQ=255)
                                                                                                                    cGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGa                                  <  2:1111592/149‑1 (MQ=255)
                                                                                                                    cGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGa                                  <  1:705059/149‑1 (MQ=255)
                                                                                                                        gTTATTTGGAATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaaa                              <  2:955221/149‑1 (MQ=255)
                                                                                                                         ttATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                             >  1:1229900/1‑149 (MQ=255)
                                                                                                                           gtttGGGATTGGCAGGGATGACTCTGGGCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAGCCATCCAGATTCTGGTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTa                           <  2:101994/148‑1 (MQ=255)
                                                                                                                           aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCCGCCAAAGTAACCCATCTTCATCTACTAGTTCAGAAAAACTa                           <  1:891755/149‑1 (MQ=255)
                                                                                                                               gggATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGtt                       >  1:389079/1‑149 (MQ=255)
                                                                                                                                ggATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGttt                      <  2:1191792/149‑1 (MQ=255)
                                                                                                                                 gATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTc                     >  1:614552/1‑149 (MQ=255)
                                                                                                                                    ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg                  >  2:724382/1‑149 (MQ=255)
                                                                                                                                           gATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg           <  2:970235/149‑1 (MQ=255)
                                                                                                                                           gATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACACATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAATACGTTTCCAGGGGaatg           >  1:852249/1‑147 (MQ=255)
                                                                                                                                              gACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa        >  1:1087575/1‑149 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTACCGGCACCAATAATACTGTTACCCTGCCAGACCCCACCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACACCTCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGGGCACGGGACAt     >  2:885959/1‑149 (MQ=255)
                                                                                                                                                    ggTCATTCCTGCTGCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCACCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGa  <  2:432527/149‑1 (MQ=255)
                                                                                                                                                    ggTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGa  >  2:502894/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGA  >  NC_000913/2884430‑2884726

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