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)65 (1.150) 43/144 0.2 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)50 (0.890) 41/144 0.3 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 124 total aligned reads displayed.

TTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524454‑4524602
                                                                         |                                                                           
ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                           <  2:1697537‑M2/76‑3 (MQ=255)
ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                           >  1:450880‑M2/1‑74 (MQ=255)
  gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         >  1:770901‑M2/67‑76 (MQ=255)
   aGAAACAACCCTGAATAACGGAACCAATACCATCCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                        >  1:603875‑M2/1‑71 (MQ=255)
    cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                       <  1:673136‑M2/12‑1 (MQ=255)
    cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                       >  1:215817‑M2/65‑76 (MQ=255)
     ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCggg                                                                      >  2:1862109‑M2/64‑76 (MQ=255)
     aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      >  1:1005908‑M2/1‑69 (MQ=255)
     aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      >  1:271903‑M2/1‑69 (MQ=255)
      taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     >  1:460119‑M2/63‑76 (MQ=255)
      aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     >  2:345294‑M2/1‑68 (MQ=255)
       acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                    >  2:303070‑M2/1‑67 (MQ=255)
       acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                    >  2:526480‑M2/1‑67 (MQ=255)
       acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                    <  1:1876702‑M2/76‑10 (MQ=255)
         cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  1:894148‑M2/60‑76 (MQ=255)
         cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  1:998975‑M2/60‑76 (MQ=255)
         cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  2:1045836‑M2/60‑76 (MQ=255)
          acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                 <  2:807434‑M2/76‑13 (MQ=255)
              tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             >  2:1803833‑M2/1‑60 (MQ=255)
              tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             >  2:700133‑M2/1‑60 (MQ=255)
              caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                             <  1:344704‑M2/22‑1 (MQ=255)
               gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            >  1:14582‑M2/1‑59 (MQ=255)
               cccggggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                            <  1:716436‑M2/23‑1 (MQ=255)
                aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  2:682446‑M2/76‑19 (MQ=255)
                 cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                          <  1:964574‑M2/25‑1 (MQ=255)
                 aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                          >  1:1080213‑M2/1‑57 (MQ=255)
                  tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                         >  1:438970‑M2/1‑56 (MQ=255)
                   aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  1:154343‑M2/76‑22 (MQ=255)
                   aaCGGAACCAATACCATTCCGTTCCAGGCGCGGTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  1:1597947‑M2/76‑22 (MQ=255)
                   aaCGGAACCAATACCATCCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAActcataagcgctaacttaagg                                                        >  1:911960‑M2/1‑55 (MQ=255)
                     cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                      <  2:51871‑M2/76‑24 (MQ=255)
                       gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    <  1:213932‑M2/31‑1 (MQ=255)
                       gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    >  1:1724078‑M2/46‑76 (MQ=255)
                       gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    <  1:1679738‑M2/31‑1 (MQ=255)
                         cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  2:204264‑M2/44‑76 (MQ=255)
                         cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  1:1517084‑M2/44‑76 (MQ=255)
                            attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                               <  2:599892‑M2/36‑1 (MQ=255)
                            attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                               >  2:92230‑M2/41‑76 (MQ=255)
                            aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               >  1:796707‑M2/1‑46 (MQ=255)
                               aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:694680‑M2/1‑43 (MQ=255)
                               aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  1:709442‑M2/1‑43 (MQ=255)
                               aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:634570‑M2/1‑43 (MQ=255)
                                ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           <  1:953587‑M2/40‑1 (MQ=255)
                                 cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  2:159576‑M2/1‑41 (MQ=255)
                                  aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         <  1:1287741‑M2/76‑37 (MQ=255)
                                    tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                       >  2:121382‑M2/1‑38 (MQ=255)
                                    tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgggattacgcct                                       <  1:1185800‑M2/76‑39 (MQ=255)
                                    gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  1:1230911‑M2/44‑1 (MQ=255)
                                     atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                      <  2:1481512‑M2/45‑1 (MQ=255)
                                        ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   <  1:1750553‑M2/76‑43 (MQ=255)
                                        gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   <  1:124436‑M2/48‑1 (MQ=255)
                                        gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   <  1:1219893‑M2/48‑1 (MQ=255)
                                        gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   <  1:449245‑M2/48‑1 (MQ=255)
                                         ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                  <  1:1579750‑M2/49‑1 (MQ=255)
                                         tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  <  1:535033‑M2/76‑44 (MQ=255)
                                         tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  <  1:1402967‑M2/76‑44 (MQ=255)
                                          tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                 <  2:1676888‑M2/50‑1 (MQ=255)
                                          tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                 >  2:858317‑M2/27‑76 (MQ=255)
                                          tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                 <  2:446017‑M2/50‑1 (MQ=255)
                                          ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  1:1733268‑M2/1‑32 (MQ=255)
                                          ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  2:1180789‑M2/1‑32 (MQ=255)
                                          ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 <  2:318690‑M2/76‑45 (MQ=255)
                                            aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                               <  1:682851‑M2/76‑47 (MQ=255)
                                              cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                             <  2:1356354‑M2/54‑1 (MQ=255)
                                               cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  1:1525332‑M2/1‑27 (MQ=255)
                                               cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            <  2:1587862‑M2/76‑50 (MQ=255)
                                               cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  1:1001926‑M2/1‑27 (MQ=255)
                                                gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                           >  1:752613‑M2/1‑26 (MQ=255)
                                                  taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         <  1:579626‑M2/58‑1 (MQ=255)
                                                   nagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  1:544602‑M2/18‑76 (MQ=255)
                                                   aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  2:1226858‑M2/18‑76 (MQ=255)
                                                    tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:1801223‑M2/1‑22 (MQ=255)
                                                    agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:528049‑M2/60‑1 (MQ=255)
                                                     aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                      <  1:52166‑M2/76‑56 (MQ=255)
                                                       tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg                    >  1:1571429‑M2/1‑19 (MQ=255)
                                                       tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg                    >  2:1049645‑M2/1‑19 (MQ=255)
                                                       tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg                    <  1:1605127‑M2/76‑58 (MQ=255)
                                                        agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                   <  2:1309319‑M2/64‑1 (MQ=255)
                                                         tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  <  1:1790889‑M2/76‑60 (MQ=255)
                                                         tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgan                  >  1:1003618‑M2/1‑17 (MQ=255)
                                                          gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 >  2:1752823‑M2/1‑16 (MQ=255)
                                                            cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               >  1:18873‑M2/9‑76 (MQ=255)
                                                              tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  2:789052‑M2/70‑1 (MQ=255)
                                                               cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta            >  2:608171‑M2/1‑11 (MQ=255)
                                                               cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta            >  2:1604242‑M2/1‑11 (MQ=255)
                                                               atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            <  1:1305219‑M2/71‑1 (MQ=255)
                                                                 gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  1:371092‑M2/73‑1 (MQ=255)
                                                                 gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  1:381821‑M2/73‑1 (MQ=255)
                                                                  ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc         <  1:1259393‑M2/76‑69 (MQ=255)
                                                                  ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc         <  2:886138‑M2/76‑69 (MQ=255)
                                                                  ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  1:1239405‑M2/74‑1 (MQ=255)
                                                                  ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  2:1374709‑M2/74‑1 (MQ=255)
                                                                  ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  2:7597‑M2/74‑1 (MQ=255)
                                                                  ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  2:1702997‑M2/74‑1 (MQ=255)
                                                                   gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct        >  2:742150‑M2/1‑7 (MQ=255)
                                                                      gCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcac     >  1:753639‑M2/1‑4 (MQ=255)
                                                                      gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  2:1517084/76‑1 (MQ=255)
                                                                        aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:885474/1‑76 (MQ=255)
                                                                        aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:1849840/1‑76 (MQ=255)
                                                                         acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  <  2:1164365‑M2/76‑76 (MQ=255)
                                                                         |                                                                           
TTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524454‑4524602

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