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)

Only 100 of 107 total aligned reads displayed.

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602
                                                                           |                                                                           
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                             <  2:1575148‑M2/6‑1 (MQ=255)
aGTTCAGCAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                             >  1:1140126/1‑76 (MQ=255)
aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                             >  2:607854/1‑76 (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                                                                    >  1:1200023‑M2/62‑76 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    <  2:259748‑M2/15‑1 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    >  2:395694‑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:1325854‑M2/1‑58 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  1:672468‑M2/1‑58 (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:1184846‑M2/37‑1 (MQ=255)
                               ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                              <  2:1069548‑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                                     >  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)
                                             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:572171‑M2/1‑27 (MQ=255)
                                                 cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            <  2:1486244‑M2/76‑50 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  1:1713017‑M2/55‑1 (MQ=255)
                                                  cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           >  2:485490‑M2/21‑76 (MQ=255)
                                                  cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           >  1:694995‑M2/21‑76 (MQ=255)
                                                   ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGGTCAGTATCAATAAc                          >  2:728603‑M2/20‑76 (MQ=255)
                                                   cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                          >  2:234533‑M2/1‑25 (MQ=255)
                                                   cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                          >  2:1740027‑M2/1‑25 (MQ=255)
                                                    taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         >  1:1617256‑M2/19‑76 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  1:1052980‑M2/18‑76 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        <  2:1200023‑M2/59‑1 (MQ=255)
                                                      agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:393535‑M2/60‑1 (MQ=255)
                                                       gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      <  1:1052541‑M2/61‑1 (MQ=255)
                                                       gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      >  1:1726641‑M2/16‑76 (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:1429226‑M2/68‑1 (MQ=255)
                                                              cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               <  1:548078‑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)
                                                                      ccGCAAcgcataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       <  1:1124599‑M2/76‑71 (MQ=255)
                                                                      ccGCAAcccataagcgctaacttgagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  2:892935‑M2/1‑6 (MQ=255)
                                                                       cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGATCAGTATCAATAACCTACCCAGGGTCAGGGACga      >  1:482614/1‑75 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  1:173662/76‑1 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  2:479268/76‑1 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  2:1468453/76‑1 (MQ=255)
                                                                         cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    >  1:1586302‑M2/1‑3 (MQ=255)
                                                                         cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa    <  1:180227/76‑1 (MQ=255)
                                                                         cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa    <  1:814010/76‑1 (MQ=255)
                                                                          aaCCCCGGGCGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCAAGGTTCAGGGACGCCAt   >  1:1377409/1‑76 (MQ=38)
                                                                           acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  1:507397‑M2/1‑1 (MQ=255)
                                                                           aCCCCGGGTGCTGCTAATGCGGATGCGACCTTTAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  <  1:1122456/76‑1 (MQ=255)
                                                                           |                                                                           
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602

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