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)95 (2.070) 67/144 0.0 98.9% 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)37 (0.810) 31/144 0.5 97.1% 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 137 total aligned reads displayed.

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524447‑4524597
                                                                           |                                                                           
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCATTATTTTGCAACCGGGGc                                                                             >  1:1355869/1‑76 (MQ=255)
 gaggcatcgtcttaacgagggcccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            <  2:304216‑M2/2‑1 (MQ=255)
  tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           >  1:645337/1‑76 (MQ=255)
     aGTTCAGAAACAACCCTGAATAACGGANCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  1:1524460/1‑76 (MQ=255)
     aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  2:799978/1‑76 (MQ=255)
      gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                       >  2:1327313‑M2/1‑75 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       >  2:912080‑M2/70‑76 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       <  1:398308‑M2/7‑1 (MQ=255)
       tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                      >  1:974595‑M2/69‑76 (MQ=255)
         cANAAAGAACCGTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCCACCGGGGCCGCACccca                                                                    >  2:1439211‑M2/1‑72 (MQ=255)
          aGAAACAACTCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                   >  1:53563‑M2/1‑71 (MQ=255)
            aaaCCACCCTGAATAACTGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAtcccnta                                                                 >  2:1443334/1‑73 (MQ=255)
            aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  1:906319‑M2/1‑69 (MQ=255)
             aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                <  2:1138960‑M2/76‑9 (MQ=255)
             aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                <  2:1021303‑M2/76‑9 (MQ=255)
              acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                               >  1:1239058‑M2/1‑67 (MQ=255)
              aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               >  1:758775‑M2/62‑76 (MQ=255)
               caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:905818‑M2/1‑66 (MQ=255)
               caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  2:637548‑M2/1‑66 (MQ=255)
                 gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                            >  1:446789‑M2/59‑76 (MQ=255)
                    cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                         >  1:328032‑M2/1‑61 (MQ=255)
                     caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        >  1:362430‑M2/55‑76 (MQ=255)
                       ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                      <  1:1172621‑M2/24‑1 (MQ=255)
                       aaTAACGGAACCAATACCNTTCCGTTCCAGGCGCGTTATATTGCAATCGGGGCCGCAAcccataagcgctaactta                                                      <  1:1458508‑M2/76‑19 (MQ=255)
                       aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      <  2:748572‑M2/76‑19 (MQ=255)
                        aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     >  1:984996‑M2/1‑57 (MQ=255)
                        aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:159392‑M2/76‑20 (MQ=255)
                         tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                    <  1:1132479‑M2/76‑21 (MQ=255)
                             cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                <  1:1183497‑M2/30‑1 (MQ=255)
                              gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                               <  1:571394‑M2/31‑1 (MQ=255)
                              gAACCACTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                               <  2:1174856‑M2/76‑26 (MQ=255)
                              gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                               >  2:1196323‑M2/1‑51 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                              >  1:1332957‑M2/1‑50 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                              <  2:1043125‑M2/76‑27 (MQ=255)
                                aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                             >  2:941130‑M2/1‑49 (MQ=255)
                                  cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                           <  2:94853‑M2/76‑30 (MQ=255)
                                   attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATACGACCt                                          <  1:86650‑M2/36‑1 (MQ=255)
                                   aaTATCTTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggga                                          <  1:92504‑M2/76‑31 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          >  2:885835‑M2/1‑46 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          <  1:234042‑M2/76‑31 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        <  1:1064256‑M2/76‑33 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgggatt                                        <  1:863080‑M2/76‑33 (MQ=255)
                                      cttcagatggttcnacccttaagttagcgcttatgggCCGCAACCTCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       <  2:1390968‑M2/39‑1 (MQ=255)
                                      aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                       >  2:518925‑M2/1‑43 (MQ=255)
                                       ttcagatggttcaacccttaagttagcgcttatgggCCGCAACANCGNGTGCTGCTAATGCGNATGCGACCTTCaa                                      >  2:1384448‑M2/37‑76 (MQ=255)
                                       ccATTCCGTTCCAGTCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcgaactgaagggttgtggtattac                                      >  2:232917‑M2/1‑42 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  1:1094659‑M2/1‑42 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  2:1046021‑M2/1‑42 (MQ=255)
                                        tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                     <  2:621321‑M2/41‑1 (MQ=255)
                                        tcagatggttcaacccttaagtnngcgcttatggnNTGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                     <  1:1374483‑M2/41‑1 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                     >  2:442606‑M2/1‑41 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTAGCAGCCNGGGCCGCAAcccatnagcgcnaacttaagggttgtggtattacg                                     >  2:192963‑M2/1‑41 (MQ=255)
                                         aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                    <  1:829956‑M2/76‑37 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                  <  2:1045076‑M2/76‑39 (MQ=255)
                                            ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 >  2:879497‑M2/1‑37 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                 <  2:46723‑M2/45‑1 (MQ=255)
                                             cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                <  1:996792‑M2/76‑41 (MQ=255)
                                             cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                <  1:424298‑M2/76‑41 (MQ=255)
                                              gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                               >  1:606836‑M2/1‑35 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  2:1329666‑M2/76‑43 (MQ=255)
                                               gttcaacgcttaagttagcgttnntgggCCGCAANNGCGGGTGCAGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              <  1:1382600‑M2/48‑1 (MQ=255)
                                                  caacccttaagttngcgcttatgggCCGNAANACCGGGTGCTGCTAATACGGATGCGACCTTCAAGGTTCAGTATc                           <  2:1379975‑M2/51‑1 (MQ=255)
                                                   aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                          <  1:1048001‑M2/52‑1 (MQ=255)
                                                   aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                          <  1:304545‑M2/76‑47 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:571874‑M2/1‑27 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:767600‑M2/1‑27 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgataagtttta                       >  1:87737‑M2/1‑27 (MQ=255)
                                                      ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  2:657760‑M2/55‑1 (MQ=255)
                                                      ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  2:657752‑M2/55‑1 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  2:1135823‑M2/1‑26 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      >  2:417251‑M2/21‑76 (MQ=255)
                                                        cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     <  1:504863‑M2/76‑52 (MQ=255)
                                                         gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                    <  2:246455‑M2/76‑53 (MQ=255)
                                                         gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                    >  1:789310‑M2/1‑24 (MQ=255)
                                                          ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                   >  1:1330601‑M2/1‑23 (MQ=255)
                                                          ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                   >  2:927220‑M2/1‑23 (MQ=255)
                                                          ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                   >  2:526667‑M2/1‑23 (MQ=255)
                                                           tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  >  2:1174417‑M2/1‑22 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  >  2:1104983‑M2/17‑76 (MQ=255)
                                                            gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 <  1:1005567‑M2/61‑1 (MQ=255)
                                                            gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGNTCAGTATCAATAACCTAc                 <  1:1455362‑M2/61‑1 (MQ=255)
                                                            aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                 >  1:1212618‑M2/1‑21 (MQ=255)
                                                             ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                >  2:796268‑M2/15‑76 (MQ=255)
                                                              gttgCATCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg               >  2:1497318‑M2/2‑19 (MQ=255)
                                                               ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatganttaacgtgccga              >  1:21057‑M2/1‑18 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat             >  1:1138961‑M2/1‑17 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat             >  1:1021304‑M2/1‑17 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat             <  2:330071‑M2/76‑60 (MQ=255)
                                                                  gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt           <  1:512237‑M2/67‑1 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:1123652‑M2/1‑15 (MQ=255)
                                                                   aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggaattacgcctgatatgatttaacgtgccgatgaa          <  2:690569‑M2/76‑63 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:585773‑M2/8‑76 (MQ=255)
                                                                     tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        >  2:91357‑M2/7‑76 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:993009‑M2/1‑12 (MQ=255)
                                                                      cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  2:893503‑M2/1‑11 (MQ=255)
                                                                        gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact     >  2:1059542‑M2/1‑9 (MQ=255)
                                                                         ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    <  1:315947‑M2/76‑69 (MQ=255)
                                                                         ggCCGCAACNCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa    >  1:33682‑M2/3‑76 (MQ=255)
                                                                          gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct   >  2:643642‑M2/1‑7 (MQ=255)
                                                                           ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  >  1:1279254‑M2/1‑6 (MQ=255)
                                                                           |                                                                           
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524447‑4524597

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.