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)62 (1.260) 44/144 0.1 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)40 (0.810) 34/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 104 total aligned reads displayed.

GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602
                                                                          |                                                                           
atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                            <  1:293819‑M2/7‑1 (MQ=255)
 tcgtcttaatgatgcaccgaggcttcgcattcttcagatggttcaacccttaagttagngcttatgggCCGCAAcc                                                                           <  1:208110‑M2/8‑1 (MQ=255)
  tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                          >  2:1585271‑M2/1‑73 (MQ=255)
  cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                          >  1:455664‑M2/68‑76 (MQ=255)
  cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                          >  1:464822‑M2/68‑76 (MQ=255)
   gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         >  1:1627035‑M2/67‑76 (MQ=255)
   cAGAAACAACCCTGAATAACGGAACCANTACCATTCTGTGCCAGGCGCCTTATTTTGCAACCGGGGCCGCAAccca                                                                         <  1:1499165‑M2/76‑5 (MQ=255)
    tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                        <  1:640748‑M2/11‑1 (MQ=255)
    aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                        >  1:694891‑M2/1‑71 (MQ=255)
    aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                        >  2:161366‑M2/1‑71 (MQ=255)
     gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  2:479064‑M2/1‑70 (MQ=255)
     gAAACAACCCTGAATAACGCAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       <  1:1447948‑M2/76‑7 (MQ=255)
       taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgtgCCGCAACCCCGGGt                                                                     >  1:1413811‑M2/63‑76 (MQ=255)
       taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     >  1:669513‑M2/63‑76 (MQ=255)
       aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     <  1:1559405‑M2/76‑9 (MQ=255)
        aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    >  2:314822‑M2/62‑76 (MQ=255)
          cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  2:791608‑M2/60‑76 (MQ=255)
          cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  2:334179‑M2/60‑76 (MQ=255)
           gaggcgccgaggcgtcgcatacttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                 <  2:1485302‑M2/18‑1 (MQ=255)
             ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                               >  2:1178209‑M2/57‑76 (MQ=255)
              gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                              <  2:1418739‑M2/21‑1 (MQ=255)
              cTGAATAACGGAACCTATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                              >  1:1415402‑M2/1‑61 (MQ=255)
              cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                              >  1:1415324‑M2/1‑61 (MQ=255)
               tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             >  2:200187‑M2/1‑60 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  1:473585‑M2/76‑19 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  1:418696‑M2/76‑19 (MQ=255)
                  aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                          >  2:837943‑M2/1‑57 (MQ=255)
                   tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                         >  1:178604‑M2/1‑56 (MQ=255)
                   gaggcgtcgcattcntcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCTGGGGCTGCTAATGCg                                                         <  1:1469261‑M2/26‑1 (MQ=255)
                    aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTGGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        >  2:712865‑M2/1‑55 (MQ=255)
                     aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                       <  2:1493035‑M2/76‑23 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:512307‑M2/29‑1 (MQ=255)
                       ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  1:869610‑M2/76‑25 (MQ=255)
                       cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                     >  2:531891‑M2/47‑76 (MQ=255)
                          cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  2:1584167‑M2/44‑76 (MQ=255)
                          aCCAATACCATTCCGTTCCATGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagctctaacttaagggttgtg                                                  >  1:1698756‑M2/1‑49 (MQ=255)
                           gcattcttcagatggttcaatccttaagttagcgcgtataggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                 <  2:1365659‑M2/34‑1 (MQ=255)
                           ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                                 >  2:379939‑M2/1‑48 (MQ=255)
                           ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                                 >  2:1013353‑M2/1‑48 (MQ=255)
                           ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                                 >  2:1441724‑M2/1‑48 (MQ=255)
                            cAACACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                                >  1:774957‑M2/1‑47 (MQ=255)
                             aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               <  1:987033‑M2/76‑31 (MQ=255)
                              attccattccTNTCNAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctatcttaagggttgttgtat                                              >  2:270412‑M2/1‑45 (MQ=255)
                              aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaggggttgtggtat                                              >  2:1248227‑M2/1‑45 (MQ=255)
                              aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                              >  1:624441‑M2/1‑45 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  1:364734‑M2/1‑44 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  2:514343‑M2/1‑44 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            <  2:1202826‑M2/76‑34 (MQ=255)
                                 ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           >  2:1069059‑M2/37‑76 (MQ=255)
                                 ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  1:1198066‑M2/1‑42 (MQ=255)
                                 ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  2:1314918‑M2/1‑42 (MQ=255)
                                  tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                          <  1:925813‑M2/41‑1 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          <  1:1154487‑M2/76‑36 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  2:1183541‑M2/1‑41 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTAGCANACGGGGCCGCGAcccataagcnctcacttaagggttgtggtattacg                                          <  1:255640‑M2/76‑36 (MQ=255)
                                    ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                        >  2:1161319‑M2/1‑39 (MQ=255)
                                      ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      <  1:466735‑M2/76‑40 (MQ=255)
                                      ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  2:1470581‑M2/1‑37 (MQ=255)
                                        ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                    <  2:707394‑M2/47‑1 (MQ=255)
                                        gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                    <  1:1293685‑M2/76‑42 (MQ=255)
                                        gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                    <  1:438906‑M2/76‑42 (MQ=255)
                                          tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  <  1:349434‑M2/76‑44 (MQ=255)
                                          tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  <  1:553785‑M2/76‑44 (MQ=255)
                                          tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  >  2:1129919‑M2/1‑33 (MQ=255)
                                            caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                <  1:69687‑M2/51‑1 (MQ=255)
                                                ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  2:642103‑M2/55‑1 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           >  2:1077821‑M2/21‑76 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           <  2:1627033‑M2/56‑1 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           <  1:1185463‑M2/56‑1 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           >  1:444248‑M2/21‑76 (MQ=255)
                                                    ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        >  1:1304709‑M2/1‑23 (MQ=255)
                                                    aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  2:463817‑M2/18‑76 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:623326‑M2/1‑22 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:638438‑M2/1‑22 (MQ=255)
                                                      gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      >  2:804816‑M2/16‑76 (MQ=255)
                                                       ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                     <  1:1211784‑M2/76‑57 (MQ=255)
                                                       ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                     <  1:972902‑M2/76‑57 (MQ=255)
                                                          tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  >  2:602236‑M2/1‑17 (MQ=255)
                                                          gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTNCAGCTTCAGTATCAATAACCTACCCAg                  <  2:5134‑M2/65‑1 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 >  2:692030‑M2/1‑16 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 >  2:362198‑M2/1‑16 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 <  2:312291‑M2/76‑61 (MQ=255)
                                                            gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                >  2:1467082‑M2/10‑76 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  2:1615384‑M2/70‑1 (MQ=255)
                                                               ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             <  1:573757‑M2/76‑65 (MQ=255)
                                                               ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             <  1:663829‑M2/76‑65 (MQ=255)
                                                                cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta            <  1:934511‑M2/76‑66 (MQ=255)
                                                                cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta            >  2:439911‑M2/1‑11 (MQ=255)
                                                                  gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          >  2:175623‑M2/4‑76 (MQ=255)
                                                                  gggCCGCAACCCCGGGTCGTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  2:205293‑M2/73‑1 (MQ=255)
                                                                    gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct        >  2:384183‑M2/1‑7 (MQ=255)
                                                                     ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  2:603270‑M2/1‑6 (MQ=255)
                                                                     ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  1:1530613‑M2/1‑6 (MQ=255)
                                                                     ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg       >  1:913749/1‑76 (MQ=255)
                                                                     ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg       <  1:1081505/76‑1 (MQ=255)
                                                                       gCNACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCCGGGAANTc     <  2:1399780/76‑1 (MQ=38)
                                                                        cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    <  1:674905‑M2/76‑74 (MQ=255)
                                                                         aacccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacga   <  1:1277277‑M2/76‑75 (MQ=255)
                                                                          acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  2:1463110‑M2/1‑1 (MQ=255)
                                                                          acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  1:1183652‑M2/1‑1 (MQ=255)
                                                                          |                                                                           
GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602

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