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)68 (1.150) 44/144 0.2 98.4% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 1 (0.020)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)78 (1.320) 50/144 0.1 98.6% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 1 (0.020)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 147 total aligned reads displayed.

ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597
                                                                          |                                                                           
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  2:581889‑M2/75‑76 (MQ=255)
 aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                           <  1:1745257‑M2/3‑1 (MQ=255)
 aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                           <  2:898948‑M2/3‑1 (MQ=255)
 aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                           <  2:1226404‑M2/3‑1 (MQ=255)
       tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                     >  1:1713804‑M2/1‑73 (MQ=255)
        gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    <  2:102433‑M2/10‑1 (MQ=255)
        gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    >  2:605822‑M2/67‑76 (MQ=255)
        gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    >  2:605801‑M2/67‑76 (MQ=255)
        gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    >  2:448289‑M2/67‑76 (MQ=255)
        cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    <  2:141555‑M2/76‑5 (MQ=255)
         aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                   >  1:109103‑M2/1‑71 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccccta                                                                  >  1:796428/1‑74 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  <  2:1885625‑M2/76‑7 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  >  2:1798660‑M2/1‑70 (MQ=255)
            taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                >  2:1227706‑M2/63‑76 (MQ=255)
             acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                               >  2:1690210‑M2/1‑67 (MQ=255)
              caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:392468‑M2/1‑66 (MQ=255)
              caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:431032‑M2/1‑66 (MQ=255)
                gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                            >  2:340333‑M2/59‑76 (MQ=255)
                gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                            >  1:531785‑M2/59‑76 (MQ=255)
                 aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                           >  1:1064306‑M2/58‑76 (MQ=255)
                 aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                           >  2:1115450‑M2/58‑76 (MQ=255)
                 aggcaccgaggcgtcccattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                           >  1:1443507‑M2/58‑76 (MQ=255)
                   gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                         <  1:663791‑M2/21‑1 (MQ=255)
                    caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        >  1:1522769‑M2/55‑76 (MQ=255)
                      ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                      >  1:480412‑M2/53‑76 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:1019532‑M2/76‑20 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:1822579‑M2/76‑20 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:756652‑M2/76‑20 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     >  2:1582751‑M2/1‑57 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  1:842335‑M2/29‑1 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  1:1909332‑M2/29‑1 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:1590880‑M2/29‑1 (MQ=255)
                             gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                               <  1:882348‑M2/31‑1 (MQ=255)
                                ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                            >  2:1279657‑M2/1‑48 (MQ=255)
                                 cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccatacgcgctaacttaagggtcgtggt                                           <  2:798768‑M2/76‑30 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          >  2:192191‑M2/41‑76 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          >  1:1198801‑M2/41‑76 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          <  1:598177‑M2/36‑1 (MQ=255)
                                  atgcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          <  2:1763908‑M2/36‑1 (MQ=255)
                                    tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  1:156003‑M2/1‑44 (MQ=255)
                                     cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       >  2:240517‑M2/38‑76 (MQ=255)
                                     cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       <  2:1448561‑M2/39‑1 (MQ=255)
                                      ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  2:1321728‑M2/1‑42 (MQ=255)
                                       cATTCCGTTCCAGGCGCGTTATGTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                     >  1:1497431‑M2/1‑41 (MQ=255)
                                          gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                  <  1:1275027‑M2/44‑1 (MQ=255)
                                           ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 >  1:1660100‑M2/1‑37 (MQ=255)
                                           ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 >  2:802300‑M2/1‑37 (MQ=255)
                                           atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                 <  1:1682395‑M2/45‑1 (MQ=255)
                                             gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                               <  1:619321‑M2/76‑42 (MQ=255)
                                             gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                               >  1:606825‑M2/1‑35 (MQ=255)
                                              ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  1:1736920‑M2/1‑34 (MQ=255)
                                              ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  1:2811‑M2/1‑34 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              >  1:1061401‑M2/29‑76 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              >  2:923362‑M2/29‑76 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              >  1:1660673‑M2/29‑76 (MQ=255)
                                               tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                             >  1:772566‑M2/1‑33 (MQ=255)
                                                tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            <  1:1227706‑M2/50‑1 (MQ=255)
                                                ccAGGCGCGTTATTTTGCAACCNGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            <  2:267976‑M2/76‑45 (MQ=255)
                                                 caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           <  1:207129‑M2/51‑1 (MQ=255)
                                                  aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                          <  1:760417‑M2/52‑1 (MQ=255)
                                                   acccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                         >  1:1268333‑M2/24‑76 (MQ=255)
                                                    gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  1:1225385‑M2/1‑28 (MQ=255)
                                                    gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  2:1601586‑M2/1‑28 (MQ=255)
                                                     cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:1869331‑M2/1‑27 (MQ=255)
                                                     cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:965495‑M2/1‑27 (MQ=255)
                                                     ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       >  2:467097‑M2/22‑76 (MQ=255)
                                                     ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  1:439471‑M2/55‑1 (MQ=255)
                                                      gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  2:378131‑M2/1‑26 (MQ=255)
                                                      gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  1:1669723‑M2/1‑26 (MQ=255)
                                                      cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      <  1:1692044‑M2/56‑1 (MQ=255)
                                                       ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     <  1:204896‑M2/57‑1 (MQ=255)
                                                        taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                    >  2:1560239‑M2/19‑76 (MQ=255)
                                                        taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                    >  1:1803644‑M2/19‑76 (MQ=255)
                                                         ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                   >  2:1933010‑M2/1‑23 (MQ=255)
                                                         aagttagcgcttatgggCCGCAACCCCGGGTGCTCGTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                   <  2:1741044‑M2/59‑1 (MQ=255)
                                                          tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  2:1679378‑M2/76‑55 (MQ=255)
                                                          tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  2:540998‑M2/76‑55 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:1220149‑M2/60‑1 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:294961‑M2/60‑1 (MQ=255)
                                                           gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 <  1:1905627‑M2/61‑1 (MQ=255)
                                                              agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa              >  2:970473‑M2/13‑76 (MQ=255)
                                                               tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat             <  2:506571‑M2/76‑60 (MQ=255)
                                                               gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             <  2:1808674‑M2/65‑1 (MQ=255)
                                                               gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             <  1:133182‑M2/65‑1 (MQ=255)
                                                                gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            <  1:1031115‑M2/76‑61 (MQ=255)
                                                                  aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa          <  1:1260305‑M2/76‑63 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  2:1610805‑M2/8‑76 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:1042322‑M2/8‑76 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:1582554‑M2/8‑76 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  2:682582‑M2/8‑76 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         <  1:395683‑M2/69‑1 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  2:52693‑M2/8‑76 (MQ=255)
                                                                    tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        >  1:1089784‑M2/7‑76 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:317145‑M2/1‑12 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        <  1:89282‑M2/76‑65 (MQ=255)
                                                                      tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      >  2:1139284‑M2/5‑76 (MQ=255)
                                                                        ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa    <  2:740112‑M2/74‑1 (MQ=255)
                                                                        ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa    <  2:1014524‑M2/74‑1 (MQ=255)
                                                                          ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  >  2:522680‑M2/1‑6 (MQ=255)
                                                                          |                                                                           
ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597

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