Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,542,020 IS5 (–) +4 bp 65.7% intergenic (+461/‑14) fimB → / → fimE regulator for fimA/regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)66 (0.580) 51/174 NT 61.7% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 4542023 39 (0.360)intergenic (+464/‑14) fimB/fimE regulator for fimA/regulator for fimA
* ? NC_000913 = 2290113NA (NA)87 (0.750) 64/176 NT 66.9% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4542020 = 41 (0.370)intergenic (+461/‑17) fimB/fimE regulator for fimA/regulator for fimA

Only 100 of 258 total aligned reads displayed.

GCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATT  >  NC_000913/4541926‑4542121
                                                                                                 |                                                                                                  
tcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAgg                                                                                                   >  1:1262326‑M2/95‑99 (MQ=255)
  cacaTCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGgga                                                                                                 <  1:1102264‑M2/99‑4 (MQ=255)
  cacaTCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGgga                                                                                                 <  1:1231409‑M2/99‑4 (MQ=255)
    caTCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGta                                                                                               <  1:2441644/99‑1 (MQ=255)
    caTCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGta                                                                                               <  2:3048617/99‑1 (MQ=255)
     aTCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGtac                                                                                              <  2:1981662/99‑1 (MQ=255)
       cACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtg                                                                                            >  1:2075608‑M2/1‑91 (MQ=255)
       cACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGtacta                                                                                            <  1:2805923/99‑1 (MQ=255)
          tCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTActacta                                                                                         <  2:476904/99‑1 (MQ=255)
           ccGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaa                                                                                        >  2:585472‑M2/1‑87 (MQ=255)
             gtgtgtGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaata                                                                                      <  1:1793157‑M2/99‑15 (MQ=255)
              ttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATtgt                                                                                     >  1:22328‑M2/81‑99 (MQ=255)
              tgtgtgGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATtgt                                                                                     >  1:2852122/1‑99 (MQ=255)
               gtgtgGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATtgtg                                                                                    >  1:767253/1‑99 (MQ=255)
                    gTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcgggg                                                                               <  1:2007075‑M2/99‑22 (MQ=255)
                      tcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAAc                                                                             <  1:1651768‑M2/27‑1 (MQ=255)
                      tcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAAc                                                                             >  1:884330‑M2/73‑99 (MQ=255)
                        tCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaat                                                                           >  1:1416757‑M2/1‑74 (MQ=255)
                           atgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAAcgtcgt                                                                        <  1:347693‑M2/32‑1 (MQ=255)
                              aaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTAt                                                                     <  1:2530362‑M2/35‑1 (MQ=255)
                               aTCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATc                                                                    <  2:610996/99‑1 (MQ=255)
                                tCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcg                                                                   <  2:2924800‑M2/99‑34 (MQ=255)
                                ccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCt                                                                   <  1:953135‑M2/37‑1 (MQ=255)
                                 cTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCtt                                                                  <  1:57567/99‑1 (MQ=255)
                                  tATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggc                                                                 <  2:2105031‑M2/99‑36 (MQ=255)
                                  ctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTa                                                                 <  1:1242101‑M2/39‑1 (MQ=255)
                                    ttGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctg                                                               <  2:303151‑M2/99‑38 (MQ=255)
                                     tGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctga                                                              <  1:1025508‑M2/99‑39 (MQ=255)
                                     tGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctga                                                              >  2:2124144‑M2/1‑61 (MQ=255)
                                      gggCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgac                                                             >  2:1205010‑M2/1‑60 (MQ=255)
                                       ggCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgact                                                            <  1:1320123‑M2/99‑41 (MQ=255)
                                       ggCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGt                                                            <  2:1070065/99‑1 (MQ=255)
                                       ctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGt                                                            >  2:1380037‑M2/56‑99 (MQ=255)
                                       ctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGt                                                            <  1:1393946‑M2/44‑1 (MQ=255)
                                         cTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTaa                                                          >  1:446691/1‑99 (MQ=255)
                                          tggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTaaa                                                         >  2:2838031‑M2/53‑99 (MQ=255)
                                          tAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcag                                                         <  1:211692‑M2/99‑44 (MQ=255)
                                          tAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcag                                                         >  1:1622793‑M2/1‑56 (MQ=255)
                                            gctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAaaga                                                       >  1:2402663‑M2/51‑99 (MQ=255)
                                             ttttGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtca                                                      <  2:670199‑M2/99‑47 (MQ=255)
                                             ttttGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAaagaa                                                      >  2:2569903/1‑99 (MQ=255)
                                              tttGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAaagaag                                                     <  2:21753/99‑1 (MQ=255)
                                              tttGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAaagaag                                                     >  1:2695795/1‑99 (MQ=255)
                                               ttGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGt                                                    <  1:565480/99‑1 (MQ=255)
                                               ccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGt                                                    <  2:2064070‑M2/52‑1 (MQ=255)
                                                cagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGtt                                                   >  2:1368792‑M2/47‑99 (MQ=255)
                                                  aCCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCa                                                 <  1:2535680/99‑1 (MQ=255)
                                                      aTTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGcc                                             >  1:488912/1‑99 (MQ=255)
                                                          acatgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCatga                                         <  1:271167‑M2/63‑1 (MQ=255)
                                                            gTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatg                                       >  1:1576117‑M2/1‑38 (MQ=255)
                                                            atgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCatgatg                                       >  1:2325895‑M2/35‑99 (MQ=255)
                                                            atgatctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCatgatg                                       <  2:1872132‑M2/65‑1 (MQ=255)
                                                             tttCCTATAGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgt                                      >  2:1488544‑M2/1‑37 (MQ=255)
                                                                 ctcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGc                                  >  1:1720435‑M2/30‑99 (MQ=255)
                                                                  tcatatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCg                                 >  2:580006‑M2/29‑99 (MQ=255)
                                                                    taGGTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagcc                               <  1:1948522‑M2/99‑70 (MQ=255)
                                                                    atatcagggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGt                               >  1:53442‑M2/27‑99 (MQ=255)
                                                                       gTATTCATTCAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgat                            >  1:2922509‑M2/1‑27 (MQ=255)
                                                                         attcattcAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattt                          <  1:2768161‑M2/99‑75 (MQ=255)
                                                                         attcattcAAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattt                          <  2:3004429‑M2/99‑75 (MQ=255)
                                                                         attcattcAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGtt                          <  2:169322/99‑1 (MQ=255)
                                                                          ttcattcaAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTa                         <  2:428325/99‑1 (MQ=255)
                                                                          gggacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTa                         >  1:1351019‑M2/21‑99 (MQ=255)
                                                                           tcattcaAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttt                        >  1:2950338‑M2/1‑23 (MQ=255)
                                                                            gacttgttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACg                       >  2:744743‑M2/19‑99 (MQ=255)
                                                                            cattcaAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgatttttt                       <  2:903298‑M2/99‑78 (MQ=255)
                                                                              ttcaAATATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttct                     <  2:1565760‑M2/99‑80 (MQ=255)
                                                                                  ttcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCaa                 >  1:1542002‑M2/13‑99 (MQ=255)
                                                                                  aaTATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccg                 >  1:186172‑M2/1‑16 (MQ=255)
                                                                                  aaTATATCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccg                 <  2:752053‑M2/99‑84 (MQ=255)
                                                                                   tcgcaccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAAc                <  1:246679‑M2/88‑1 (MQ=255)
                                                                                    tatatCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACg               >  2:2955562/1‑99 (MQ=255)
                                                                                      tatCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaa             <  2:2742907‑M2/99‑88 (MQ=255)
                                                                                      tatCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaa             <  2:2339269‑M2/99‑88 (MQ=255)
                                                                                      tatCTCAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaa             <  2:1292106‑M2/99‑88 (MQ=255)
                                                                                       accttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGa            >  2:712801‑M2/8‑99 (MQ=255)
                                                                                        tctcAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatg           >  2:2252422‑M2/1‑10 (MQ=255)
                                                                                        ccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAg           <  2:2784653‑M2/93‑1 (MQ=255)
                                                                                        ccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAg           <  2:1348620‑M2/93‑1 (MQ=255)
                                                                                        ccttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAg           <  1:822483‑M2/93‑1 (MQ=255)
                                                                                         ctcAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGc          <  2:1756564/99‑1 (MQ=255)
                                                                                          ttccTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGcc         >  2:1962166‑M2/5‑99 (MQ=255)
                                                                                          tcAGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgcc         <  1:1495407‑M2/99‑92 (MQ=255)
                                                                                            aGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgcctt       >  2:842263‑M2/1‑6 (MQ=255)
                                                                                             cTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCaga      <  2:1085683‑M2/98‑1 (MQ=255)
                                                                                              ttAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttga     <  2:2601340‑M2/99‑96 (MQ=255)
                                                                                              ttAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttga     <  2:2526040‑M2/99‑96 (MQ=255)
                                                                                              ttAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCagag     <  2:2466697/99‑1 (MQ=255)
                                                                                              ttAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCagag     >  1:636377/1‑99 (MQ=255)
                                                                                               tAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaa    >  1:1008744‑M2/1‑3 (MQ=255)
                                                                                               tAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaa    >  2:1851211‑M2/1‑3 (MQ=255)
                                                                                               tAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAgaga    <  1:2200954/99‑1 (MQ=255)
                                                                                                aGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaat   <  2:319987‑M2/99‑98 (MQ=255)
                                                                                                aGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaat   >  2:160262‑M2/1‑2 (MQ=255)
                                                                                                aGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaat   >  2:91968‑M2/1‑2 (MQ=255)
                                                                                                aGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatatcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaat   <  1:3010366‑M2/99‑98 (MQ=255)
                                                                                                 gggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatc  <  2:513162‑M2/99‑99 (MQ=255)
                                                                                                 gggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatc  <  1:1885457‑M2/99‑99 (MQ=255)
                                                                                                 gggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatc  <  1:2348705‑M2/99‑99 (MQ=255)
                                                                                                 ggAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGatt  >  1:2661755/1‑99 (MQ=255)
                                                                                                 |                                                                                                  
GCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATT  >  NC_000913/4541926‑4542121

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