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)43 (0.830) 36/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)56 (1.080) 46/144 0.1 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

TTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524449‑4524597
                                                                         |                                                                           
aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                           >  1:661734‑M2/74‑76 (MQ=255)
  gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                         <  2:1273900‑M2/5‑1 (MQ=255)
  gcatcgtcttaacgaggcaccgaggcgacgcattcttcagatggttcaacccgtaagttagcgcttatgggCCGCa                                                                         >  2:118779‑M2/72‑76 (MQ=255)
   catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                        <  2:1575148‑M2/6‑1 (MQ=255)
   aGTTCAGCAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  1:1140126/1‑76 (MQ=255)
   aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  2:607854/1‑76 (MQ=255)
        aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                   >  2:1724258‑M2/1‑71 (MQ=255)
         cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  <  2:1608517‑M2/12‑1 (MQ=255)
          aaaCAACCCTGAATAACGGACCCAATCCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCGCAAcccataag                                                                >  2:904023‑M2/1‑68 (MQ=255)
          aaaCAACCCTGAATAACGGAACCAATCCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  1:68849‑M2/1‑69 (MQ=255)
            aacggggcagcgaggcgtcgcattcttcagatggttcaacccttaagtgagcgcttatgggCCGCAACCCCGGGtg                                                               <  2:537363‑M2/15‑1 (MQ=255)
            aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               >  2:395694‑M2/62‑76 (MQ=255)
            aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               <  2:259748‑M2/15‑1 (MQ=255)
            aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               <  2:1422502‑M2/15‑1 (MQ=255)
            aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               >  1:1200023‑M2/62‑76 (MQ=255)
             caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCACcccatacgcg                                                              <  2:1542912‑M2/76‑11 (MQ=255)
             caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              <  2:1171706‑M2/76‑11 (MQ=255)
             acgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAGCCCCGGGtgc                                                              <  1:760087‑M2/16‑1 (MQ=255)
              aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             <  2:1631284‑M2/76‑12 (MQ=255)
               gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                            <  2:1137511‑M2/18‑1 (MQ=255)
                cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           <  1:1234494‑M2/76‑14 (MQ=255)
                aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                           >  2:1634689‑M2/58‑76 (MQ=255)
                 ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                          >  2:933308‑M2/1‑62 (MQ=255)
                  cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                         >  2:1463109‑M2/1‑61 (MQ=255)
                   tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                        >  1:482915‑M2/1‑60 (MQ=255)
                    gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       >  1:1528713‑M2/1‑59 (MQ=255)
                    gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       >  1:164139‑M2/1‑59 (MQ=255)
                    accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       >  2:1574909‑M2/54‑76 (MQ=255)
                     ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                      >  2:7541‑M2/53‑76 (MQ=255)
                     aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  1:672468‑M2/1‑58 (MQ=255)
                     aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  1:1325854‑M2/1‑58 (MQ=255)
                       tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                    <  1:1346779‑M2/76‑21 (MQ=255)
                           cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                <  2:548008‑M2/30‑1 (MQ=255)
                            gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGAGCTGCTAACGCGGATGc                                               >  2:100989‑M2/46‑76 (MQ=255)
                            gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggatg                                               <  2:167928‑M2/76‑26 (MQ=255)
                             tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              <  1:1035‑M2/32‑1 (MQ=255)
                                cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                           <  1:898484‑M2/35‑1 (MQ=255)
                                cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                           <  2:1216948‑M2/76‑30 (MQ=255)
                                 attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          <  1:1223255‑M2/36‑1 (MQ=255)
                                 aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          >  2:475348‑M2/1‑46 (MQ=255)
                                 aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          <  1:1484149‑M2/76‑31 (MQ=255)
                                  ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  2:1069548‑M2/37‑1 (MQ=255)
                                  ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  2:1184846‑M2/37‑1 (MQ=255)
                                  aTACCTTTCCGTTCCAGGCGCGTTATTTTGCACCCGGGGCCGCAAcccataagcgctaacttaagggttttggtat                                         >  2:1076712‑M2/1‑45 (MQ=255)
                                   tgccatTCCTTTCCAGGCGCGTTATTTTGCACCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  2:530598‑M2/3‑44 (MQ=255)
                                   tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                        <  1:1236202‑M2/38‑1 (MQ=255)
                                   tcccaaTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  2:156393‑M2/3‑44 (MQ=255)
                                   tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  1:1158263‑M2/1‑44 (MQ=255)
                                    cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       >  2:72487‑M2/38‑76 (MQ=255)
                                    cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       <  1:626989‑M2/39‑1 (MQ=255)
                                     ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  2:246066‑M2/1‑42 (MQ=255)
                                      tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                     >  2:1125140‑M2/36‑76 (MQ=255)
                                         tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttntggtattacgcct                                  >  1:1183889‑M2/1‑38 (MQ=255)
                                         gatggttcaccccttaactgagcggttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                  <  2:212752‑M2/44‑1 (MQ=255)
                                         gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                  <  2:921008‑M2/44‑1 (MQ=255)
                                         gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                  <  1:334287‑M2/44‑1 (MQ=255)
                                           tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                >  2:555557‑M2/31‑76 (MQ=255)
                                           tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                <  2:1189096‑M2/46‑1 (MQ=255)
                                           tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                <  1:213934‑M2/46‑1 (MQ=255)
                                           tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                >  2:26912‑M2/31‑76 (MQ=255)
                                             gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              <  2:866539‑M2/48‑1 (MQ=255)
                                              ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             >  2:370301‑M2/28‑76 (MQ=255)
                                               tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            <  2:1086466‑M2/50‑1 (MQ=255)
                                                caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           >  2:969782‑M2/26‑76 (MQ=255)
                                                cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                           >  1:684571‑M2/1‑31 (MQ=255)
                                                 aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                          >  1:1373270‑M2/25‑76 (MQ=255)
                                                 aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                          >  1:225372‑M2/1‑30 (MQ=255)
                                                    cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       <  2:1486244‑M2/76‑50 (MQ=255)
                                                    cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:572171‑M2/1‑27 (MQ=255)
                                                    ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  1:1713017‑M2/55‑1 (MQ=255)
                                                     cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      >  1:694995‑M2/21‑76 (MQ=255)
                                                     cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      >  2:485490‑M2/21‑76 (MQ=255)
                                                      ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGGTCAGTATCAATAAc                     >  2:728603‑M2/20‑76 (MQ=255)
                                                      cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     >  2:1740027‑M2/1‑25 (MQ=255)
                                                      cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     >  2:234533‑M2/1‑25 (MQ=255)
                                                       taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                    >  1:1617256‑M2/19‑76 (MQ=255)
                                                        aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                   <  2:1200023‑M2/59‑1 (MQ=255)
                                                        aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                   >  1:1052980‑M2/18‑76 (MQ=255)
                                                         agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:393535‑M2/60‑1 (MQ=255)
                                                          gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 >  1:1726641‑M2/16‑76 (MQ=255)
                                                          gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 <  1:1052541‑M2/61‑1 (MQ=255)
                                                           ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                <  1:212439‑M2/62‑1 (MQ=255)
                                                           ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                >  1:417593‑M2/15‑76 (MQ=255)
                                                            tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg               <  1:442706‑M2/76‑58 (MQ=255)
                                                              gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             <  1:1267217‑M2/65‑1 (MQ=255)
                                                               gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            >  2:642095‑M2/1‑16 (MQ=255)
                                                                cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:1411457‑M2/1‑15 (MQ=255)
                                                                 cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          <  1:548078‑M2/68‑1 (MQ=255)
                                                                 cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          <  1:1429226‑M2/68‑1 (MQ=255)
                                                                 aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa          <  2:1522593‑M2/76‑63 (MQ=255)
                                                                   tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCCATAACCTACCCAGGTTCa        >  2:1703972‑M2/7‑76 (MQ=255)
                                                                   tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        <  1:1430562‑M2/70‑1 (MQ=255)
                                                                   ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  2:1502288‑M2/1‑12 (MQ=255)
                                                                    cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  2:793157‑M2/1‑11 (MQ=255)
                                                                    atgggCCGCACCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCTGGTTCAg       <  1:1123479‑M2/71‑1 (MQ=255)
                                                                       ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    <  2:593916‑M2/76‑69 (MQ=255)
                                                                       ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa    <  2:179469‑M2/74‑1 (MQ=255)
                                                                        gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct   >  1:712431‑M2/1‑7 (MQ=255)
                                                                         ccGCAAcgcataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  <  1:1124599‑M2/76‑71 (MQ=255)
                                                                         ccGCAAcccataagcgctaacttgagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  >  2:892935‑M2/1‑6 (MQ=255)
                                                                         |                                                                           
TTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524449‑4524597

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