Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,542,223 IS1 (–) +9 bp 60.7% coding (187‑195/597 nt) fimE → regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3583428 =NA (NA)55 (0.580) 49/176 NT 57.8% noncoding (1/768 nt) IS1 repeat region
?NC_000913 4542223 = 36 (0.420)coding (187/597 nt) fimE regulator for fimA
* ? NC_000913 = 3584195NA (NA)56 (0.590) 46/176 NT 58.2% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 4542231 36 (0.420)coding (195/597 nt) fimE regulator for fimA

Only 100 of 175 total aligned reads displayed.

CTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCG  >  NC_000913/4542133‑4542329
                                                                                                  |                                                                                                  
cTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTg                                                                                                    <  2:1720916/99‑1 (MQ=255)
  gTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaa                                                                                                  <  2:2328834/99‑1 (MQ=255)
  gTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGgg                                                                                                  >  1:2117848‑M2/1‑97 (MQ=255)
  acaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTgaa                                                                                                  <  2:1881205‑M2/11‑1 (MQ=255)
  acaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTgaa                                                                                                  >  1:76167‑M2/89‑99 (MQ=255)
   ttGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggt                                                                                                 >  2:698380‑M2/1‑96 (MQ=255)
   caggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTgaag                                                                                                 >  1:739245‑M2/88‑99 (MQ=255)
     tcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaagaa                                                                                               <  1:2246107‑M2/76‑1 (MQ=255)
      gCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAAc                                                                                              <  2:313892/99‑1 (MQ=255)
      gCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaat                                                                                              >  1:162121‑M2/1‑93 (MQ=255)
        gggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACgg                                                                                            <  2:1368343‑M2/17‑1 (MQ=255)
          gacagctgatagaaacagaagccactggagcacctcaaaaacaccgtcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGAt                                                                                          <  2:75496‑M2/19‑1 (MQ=255)
            cagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGAttt                                                                                        >  2:2248847‑M2/79‑99 (MQ=255)
             agctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGAtttt                                                                                       <  2:279069‑M2/22‑1 (MQ=255)
              tgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTc                                                                                      >  1:1875899‑M2/15‑99 (MQ=255)
              gctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTc                                                                                      <  1:633368‑M2/23‑1 (MQ=255)
               cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCt                                                                                     <  2:725856/99‑1 (MQ=255)
               cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaac                                                                                     >  2:1693366‑M2/1‑84 (MQ=255)
               cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCt                                                                                     >  1:1944505/1‑99 (MQ=255)
                aTGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaact                                                                                    >  2:284318‑M2/1‑83 (MQ=255)
                 tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaactt                                                                                   <  1:1158260‑M2/99‑18 (MQ=255)
                 gatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTAc                                                                                   >  1:2174812‑M2/74‑99 (MQ=255)
                  gggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaactta                                                                                  >  2:405300‑M2/1‑81 (MQ=255)
                   tagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCg                                                                                 >  1:1362338‑M2/72‑99 (MQ=255)
                      tGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattga                                                                              >  2:215674‑M2/1‑77 (MQ=255)
                      tGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattga                                                                              <  1:282419‑M2/99‑23 (MQ=255)
                       gCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCa                                                                             <  2:637161/99‑1 (MQ=255)
                        acagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCAc                                                                            >  1:1515663‑M2/67‑99 (MQ=255)
                         gTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatag                                                                           <  2:1812186‑M2/99‑26 (MQ=255)
                           gaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCg                                                                         >  1:2199229‑M2/64‑99 (MQ=255)
                            ttAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgt                                                                        >  1:2018093‑M2/1‑71 (MQ=255)
                            aagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGt                                                                        <  2:934330‑M2/37‑1 (MQ=255)
                                  aaCTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatg                                                                  <  2:1057966‑M2/99‑35 (MQ=255)
                                   tggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGttt                                                                 <  2:1600052‑M2/44‑1 (MQ=255)
                                   aCTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGttt                                                                 >  2:506198/1‑99 (MQ=255)
                                    ggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGtttt                                                                >  2:2108904‑M2/55‑99 (MQ=255)
                                      gCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtga                                                              >  2:1275397/1‑99 (MQ=255)
                                        ttGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcaga                                                            >  2:1666686‑M2/1‑59 (MQ=255)
                                         tGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagat                                                           >  2:1210384‑M2/1‑58 (MQ=255)
                                          cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAg                                                          <  1:2092354‑M2/51‑1 (MQ=255)
                                          cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAg                                                          <  1:304787‑M2/51‑1 (MQ=255)
                                            tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg                                                        >  1:301262/1‑99 (MQ=255)
                                            tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg                                                        <  1:1836101/99‑1 (MQ=255)
                                             caaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGt                                                       >  1:530368‑M2/46‑99 (MQ=255)
                                              tGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTg                                                      <  2:1620090/99‑1 (MQ=255)
                                               aaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGa                                                     >  1:23233‑M2/44‑99 (MQ=255)
                                               aaaacaccatcatacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGa                                                     >  1:561590‑M2/44‑99 (MQ=255)
                                                cATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgccc                                                    <  2:2482968‑M2/99‑49 (MQ=255)
                                                 aTTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAg                                                   <  1:100852/99‑1 (MQ=255)
                                                      cgagaccttgaccttaatGGAGGTAGAATAAATTCTTCCCGAATAggtaatgactccaacttattgatagtgttttatgtgcagataatgcccgatgac                                              >  2:1336215‑M2/4‑45 (MQ=255)
                                                        ggACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGa                                            >  2:2415900/1‑99 (MQ=255)
                                                         gaccttgacctcAATGAAGGTAGACTGAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGaa                                           >  2:2272390/1‑99 (MQ=255)
                                                         catacactaaatcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGaa                                           <  2:1688275‑M2/66‑1 (MQ=255)
                                                          accttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                          <  1:475336‑M2/99‑59 (MQ=255)
                                                           ccttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgt                                         >  1:1454363‑M2/1‑40 (MQ=255)
                                                             ttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtca                                       >  2:1835786‑M2/1‑38 (MQ=255)
                                                             ttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCt                                       <  1:1730935/99‑1 (MQ=255)
                                                               gaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatg                                     <  2:746493‑M2/99‑64 (MQ=255)
                                                                accttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgc                                    <  2:1013228‑M2/99‑65 (MQ=255)
                                                                   ttAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGAccc                                 >  1:1306077/1‑99 (MQ=255)
                                                                   atcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGAccc                                 <  1:1648383‑M2/76‑1 (MQ=255)
                                                                   atcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGAccc                                 >  2:1955104‑M2/24‑99 (MQ=255)
                                                                    tcagtaagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCa                                >  2:1140311‑M2/23‑99 (MQ=255)
                                                                     aaTGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctc                               >  1:1964818‑M2/1‑30 (MQ=255)
                                                                       tGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcca                             <  2:130325‑M2/99‑72 (MQ=255)
                                                                        taagttggcagcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGaa                            >  1:1843133‑M2/19‑99 (MQ=255)
                                                                        gAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccac                            <  1:1420019‑M2/99‑73 (MQ=255)
                                                                         aaGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc                           <  2:68797‑M2/99‑74 (MQ=255)
                                                                         aaGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc                           >  2:863691‑M2/1‑26 (MQ=255)
                                                                          aGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACg                          >  2:2042856/1‑99 (MQ=255)
                                                                           ggTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccga                         <  2:1659445‑M2/99‑76 (MQ=255)
                                                                           ggTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccga                         >  2:1231918‑M2/1‑24 (MQ=255)
                                                                            gTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgat                        <  1:1913322‑M2/99‑77 (MQ=255)
                                                                             tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt                       >  1:2510714‑M2/1‑22 (MQ=255)
                                                                             tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt                       >  2:2229521‑M2/1‑22 (MQ=255)
                                                                             tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt                       <  1:464800‑M2/99‑78 (MQ=255)
                                                                             tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt                       >  1:671070‑M2/1‑22 (MQ=255)
                                                                             tAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGc                       >  1:1699070/1‑99 (MQ=255)
                                                                               gAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTa                     >  2:196526/1‑99 (MQ=255)
                                                                                 agcatcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAAc                   <  2:352421‑M2/90‑1 (MQ=255)
                                                                                   aaataTTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgaga                 >  1:746557‑M2/1‑16 (MQ=255)
                                                                                    atcaccCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTgg                >  2:333952‑M2/7‑99 (MQ=255)
                                                                                    aataTTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaa                <  1:1060028‑M2/99‑85 (MQ=255)
                                                                                        ttCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgac            <  1:435871‑M2/99‑89 (MQ=255)
                                                                                        ttCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAg            >  2:1850734/1‑99 (MQ=255)
                                                                                        ttCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAg            >  1:1801085/1‑99 (MQ=255)
                                                                                         tCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAgg           <  2:27072/99‑1 (MQ=255)
                                                                                          cGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGgc          <  2:457482/99‑1 (MQ=255)
                                                                                            ccGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGgcgc        <  2:2470820/99‑1 (MQ=255)
                                                                                             cGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCt       <  2:778985/99‑1 (MQ=255)
                                                                                              gACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTg      >  1:2285936/1‑99 (MQ=255)
                                                                                               aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGTAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGa     >  1:1866155/1‑99 (MQ=255)
                                                                                               aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGa     <  1:237917/99‑1 (MQ=255)
                                                                                               aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGa     <  1:807420/99‑1 (MQ=255)
                                                                                                cTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGAc    >  1:1904710/1‑99 (MQ=255)
                                                                                                 tGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttc   <  2:1908423‑M2/99‑98 (MQ=255)
                                                                                                 tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGAcc   >  1:934923/1‑99 (MQ=255)
                                                                                                 tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGAcc   >  2:2335686/1‑99 (MQ=255)
                                                                                                  gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCg  <  2:1003249/99‑1 (MQ=255)
                                                                                                  gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCg  >  2:251961/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
CTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCG  >  NC_000913/4542133‑4542329

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