Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp coding (494‑499/549 nt) fimA → major type 1 subunit fimbrin (pilin)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2772854 =NA (NA)179 (0.950) 84/144 0.4 100% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 0 (0.000)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)210 (1.120) 93/144 0.2 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 410 total aligned reads displayed.

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602
                                                                           |                                                                           
aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                             >  2:4373889/1‑76 (MQ=255)
 gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                            <  2:3975768‑M2/76‑2 (MQ=255)
 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                            <  1:5933167‑M2/7‑1 (MQ=255)
  ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                           <  2:2598913‑M2/76‑3 (MQ=255)
  tcgtcttaacgaggcagcgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                           <  2:4600228‑M2/8‑1 (MQ=255)
   tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                          >  2:1538653‑M2/1‑73 (MQ=255)
   tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                          >  1:1611372‑M2/1‑73 (MQ=255)
   tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                          >  1:6143956‑M2/1‑73 (MQ=255)
    gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         <  1:1900240‑M2/10‑1 (MQ=255)
    cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         <  2:5003785‑M2/76‑5 (MQ=255)
    cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         <  2:5188405‑M2/76‑5 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  1:6314936‑M2/1‑70 (MQ=255)
      cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                       >  2:2307205‑M2/65‑76 (MQ=255)
        taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     <  2:3792561‑M2/14‑1 (MQ=255)
        taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     >  2:369062‑M2/63‑76 (MQ=255)
        taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     >  1:5593833‑M2/63‑76 (MQ=255)
        aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     <  1:6092883‑M2/76‑9 (MQ=255)
        aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     <  1:1523346‑M2/76‑9 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    <  2:5371880‑M2/15‑1 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                  <  1:4248054‑M2/76‑12 (MQ=255)
            gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                 <  1:3370848‑M2/18‑1 (MQ=255)
            gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                 >  2:1868568‑M2/59‑76 (MQ=255)
            acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                 >  1:5015133‑M2/1‑64 (MQ=255)
            acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                 >  2:4207426‑M2/1‑64 (MQ=255)
             cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                >  1:2193515‑M2/1‑63 (MQ=255)
             aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                                <  1:4370398‑M2/19‑1 (MQ=255)
              ccTGAATAACGGAACCAATCCCATGCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               >  1:784868‑M2/1‑62 (MQ=255)
               cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                              >  2:5047261‑M2/1‑61 (MQ=255)
                tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             >  1:5726005‑M2/1‑60 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            >  2:6289136‑M2/1‑59 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            <  1:2269738‑M2/76‑18 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            <  1:6036814‑M2/76‑18 (MQ=255)
                  ccgaggcgtcgcattcttcagatggtttaacccttaagttagcgcttatgggCCGCAACCCCGGATGCTGCTAATg                                                           >  2:6348175‑M2/53‑76 (MQ=255)
                  ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                           <  2:5459883‑M2/24‑1 (MQ=255)
                  aaTAACGGAACCAATCCCATTCCGTTCAAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  2:3183488‑M2/1‑58 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  2:5501922‑M2/1‑58 (MQ=255)
                   cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                          <  1:401776‑M2/25‑1 (MQ=255)
                   cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                          <  1:3542311‑M2/25‑1 (MQ=255)
                    gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         >  2:3238507‑M2/51‑76 (MQ=255)
                     aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                        <  1:4818202‑M2/27‑1 (MQ=255)
                     aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                        <  1:29734‑M2/27‑1 (MQ=255)
                     aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                        <  1:901268‑M2/27‑1 (MQ=255)
                     aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  1:2667899‑M2/76‑22 (MQ=255)
                       cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                      <  1:1730102‑M2/76‑24 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  1:2365774‑M2/76‑25 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     >  1:1867121‑M2/1‑52 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    >  1:3885925‑M2/46‑76 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   <  1:1027315‑M2/32‑1 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   >  1:4911659‑M2/45‑76 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   >  1:1230514‑M2/45‑76 (MQ=255)
                           cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  2:2928291‑M2/44‑76 (MQ=255)
                            ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                                 <  2:3503‑M2/76‑29 (MQ=255)
                             cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                                >  1:5473051‑M2/1‑47 (MQ=255)
                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               <  2:4101838‑M2/76‑31 (MQ=255)
                                tcccatTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  1:108982‑M2/3‑44 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:6207455‑M2/1‑43 (MQ=255)
                                  ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           <  2:1160245‑M2/40‑1 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  1:2795030‑M2/1‑41 (MQ=255)
                                       ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  2:5957969‑M2/1‑37 (MQ=255)
                                       atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                      <  2:1323932‑M2/45‑1 (MQ=255)
                                        tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAGGCGACCTTCAAGGTTCa                                     <  2:4932315‑M2/46‑1 (MQ=255)
                                         ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                    >  1:3233385‑M2/30‑76 (MQ=255)
                                            ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  2:6044918‑M2/1‑32 (MQ=255)
                                             caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                <  1:4135629‑M2/51‑1 (MQ=255)
                                             caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                >  1:468014‑M2/26‑76 (MQ=255)
                                             cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                                <  2:3916724‑M2/76‑46 (MQ=255)
                                              aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               >  2:4743029‑M2/25‑76 (MQ=255)
                                              aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               >  1:2944303‑M2/25‑76 (MQ=255)
                                                cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                             <  2:1098853‑M2/54‑1 (MQ=255)
                                                 ggtgtttTTATGCCACAGGGTCCGCAAcccataagcgctaacttaagggttgttgtattacgcctgatatgattta                            >  2:4847846‑M2/4‑27 (MQ=255)
                                                 cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  2:751643‑M2/1‑27 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  1:5326754‑M2/55‑1 (MQ=255)
                                                  gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                           >  1:365171‑M2/1‑26 (MQ=255)
                                                       aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                      <  2:2653378‑M2/76‑56 (MQ=255)
                                                          agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                   <  1:2106031‑M2/64‑1 (MQ=255)
                                                           tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  <  2:1624195‑M2/76‑60 (MQ=255)
                                                           tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  <  2:2838818‑M2/76‑60 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  1:4457012‑M2/65‑1 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  2:5285960‑M2/65‑1 (MQ=255)
                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 <  1:1386897‑M2/66‑1 (MQ=255)
                                                              cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               >  1:3034298‑M2/9‑76 (MQ=255)
                                                              cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGGTCAGTATCAATAACCTACCCAGGtt               >  2:307584‑M2/9‑76 (MQ=255)
                                                               ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc              >  2:401800‑M2/8‑76 (MQ=255)
                                                               gtatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc              <  1:2586459‑M2/69‑1 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  2:2048521‑M2/7‑76 (MQ=255)
                                                                ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             <  1:619909‑M2/76‑65 (MQ=255)
                                                                 atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            >  2:2549103‑M2/6‑76 (MQ=255)
                                                                 atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            <  2:5621787‑M2/71‑1 (MQ=255)
                                                                  tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg           >  2:2198681‑M2/5‑76 (MQ=255)
                                                                   gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          >  2:3318451‑M2/4‑76 (MQ=255)
                                                                     gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct        >  2:10142‑M2/1‑7 (MQ=255)
                                                                     gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc        <  1:3776891‑M2/75‑1 (MQ=255)
                                                                      ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  2:4536715‑M2/1‑6 (MQ=255)
                                                                        gCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcac     >  1:4566094‑M2/1‑4 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  2:5852818/76‑1 (MQ=255)
                                                                         cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    >  1:5123767‑M2/1‑3 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:2422882/1‑76 (MQ=255)
                                                                           acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  2:1221590‑M2/1‑1 (MQ=255)
                                                                           aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  >  2:5236942/1‑76 (MQ=255)
                                                                           aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  <  1:1466083/76‑1 (MQ=255)
                                                                           |                                                                           
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602

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