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)55 (1.180) 45/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)43 (0.920) 35/144 0.4 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602
                                                                           |                                                                           
aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                             <  2:1441822/76‑1 (MQ=255)
 gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                            >  1:174001‑M2/1‑75 (MQ=255)
 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                            <  2:1191424‑M2/7‑1 (MQ=255)
    gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         <  1:1331898‑M2/10‑1 (MQ=255)
    gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         >  2:1103240‑M2/67‑76 (MQ=255)
     gcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                        <  2:145969‑M2/11‑1 (MQ=255)
     aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                        >  1:363130‑M2/1‑71 (MQ=255)
       aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      <  2:511588‑M2/76‑8 (MQ=255)
        taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     >  1:609192‑M2/63‑76 (MQ=255)
         acaacCCTGAATAACGGAACCAATACCATTCCGTTCCACGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                    <  2:966476‑M2/76‑10 (MQ=255)
          caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                   <  1:443666‑M2/76‑11 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                  <  2:196613‑M2/76‑12 (MQ=255)
            gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                 <  2:87372‑M2/18‑1 (MQ=255)
             cccTGAAAAACGGAACCAATACCATTCCGTGCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                >  2:602381‑M2/1‑63 (MQ=255)
              ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                               >  2:930750‑M2/57‑76 (MQ=255)
              ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               >  1:1015588‑M2/1‑62 (MQ=255)
              ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               >  2:842692‑M2/1‑62 (MQ=255)
               cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                              <  1:1314713‑M2/76‑16 (MQ=255)
                tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             >  1:1260528‑M2/1‑60 (MQ=255)
                tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             <  1:640497‑M2/76‑17 (MQ=255)
                 gccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                            <  1:885891‑M2/23‑1 (MQ=255)
                  ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                           >  2:90440‑M2/53‑76 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  1:347092‑M2/76‑19 (MQ=255)
                    gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         >  2:1057290‑M2/51‑76 (MQ=255)
                    gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         >  1:108154‑M2/51‑76 (MQ=255)
                     aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                        <  1:1234231‑M2/27‑1 (MQ=255)
                     aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        >  2:480623‑M2/1‑55 (MQ=255)
                     aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  2:42978‑M2/76‑22 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     >  1:1018546‑M2/1‑52 (MQ=255)
                        cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                     <  1:369736‑M2/30‑1 (MQ=255)
                         gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccgatatgcgttcacttaaggggtg                                                    <  2:860651‑M2/76‑26 (MQ=255)
                         gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                    >  2:1308674‑M2/1‑51 (MQ=255)
                         gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                    >  1:592395‑M2/1‑51 (MQ=255)
                            gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                 <  2:750863‑M2/34‑1 (MQ=255)
                            acaATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                                 <  1:848073‑M2/75‑29 (MQ=255)
                             cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                >  1:260854‑M2/42‑76 (MQ=255)
                             cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                >  2:1290001‑M2/42‑76 (MQ=255)
                             cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                <  1:1084875‑M2/35‑1 (MQ=255)
                             cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                                <  1:842262‑M2/76‑30 (MQ=255)
                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               >  2:65444‑M2/1‑46 (MQ=255)
                                tcttcagatggttcaacccttaagttagcgcttatgggCGGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                             <  1:657800‑M2/38‑1 (MQ=255)
                                tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                             <  1:931029‑M2/38‑1 (MQ=255)
                                tACCATGCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  1:862598‑M2/1‑44 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:233778‑M2/1‑43 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTACAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  1:1337591‑M2/1‑43 (MQ=255)
                                  ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  1:937545‑M2/1‑42 (MQ=255)
                                  ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  1:598475‑M2/1‑42 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          <  2:1039899‑M2/76‑36 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  2:76135‑M2/1‑41 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         <  1:340986‑M2/76‑37 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         >  2:1259308‑M2/1‑40 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         <  1:1523208‑M2/76‑37 (MQ=255)
                                      gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       >  2:292945‑M2/33‑76 (MQ=255)
                                       atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                      >  2:765732‑M2/32‑76 (MQ=255)
                                        cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                     <  1:1001607‑M2/76‑41 (MQ=255)
                                         gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                    >  2:85235‑M2/1‑35 (MQ=255)
                                          gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   >  1:22939‑M2/29‑76 (MQ=255)
                                             naacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                >  1:83767‑M2/26‑76 (MQ=255)
                                              aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               <  1:234405‑M2/52‑1 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             >  1:152719‑M2/1‑28 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             >  1:44974‑M2/1‑28 (MQ=255)
                                                cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATTCGGATGCGACCTTCAAGGTTCAGTATCAAt                             <  2:1247977‑M2/54‑1 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  2:1104818‑M2/55‑1 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  2:355042‑M2/55‑1 (MQ=255)
                                                  gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                           >  1:142982‑M2/1‑26 (MQ=255)
                                                  gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                           <  1:613397‑M2/76‑51 (MQ=255)
                                                  cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           >  2:1488268‑M2/21‑76 (MQ=255)
                                                   cgTTATTTTGCAANCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                          <  1:1033498‑M2/76‑52 (MQ=255)
                                                   cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactaaagggttgtggcattacgcctgatatactttaac                          <  1:1412508‑M2/76‑52 (MQ=255)
                                                    taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         >  1:887240‑M2/19‑76 (MQ=255)
                                                    taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         >  2:1328444‑M2/19‑76 (MQ=255)
                                                     ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        >  1:574545‑M2/1‑23 (MQ=255)
                                                     ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        <  2:158793‑M2/76‑54 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        <  1:119608‑M2/59‑1 (MQ=255)
                                                      tATTTTGCAACCGGGGCCTCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:417079‑M2/1‑22 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       <  1:1165685‑M2/76‑55 (MQ=255)
                                                      agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       >  1:389904‑M2/17‑76 (MQ=255)
                                                       gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      <  1:491065‑M2/61‑1 (MQ=255)
                                                          agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAGCTACCCa                   <  1:966290‑M2/64‑1 (MQ=255)
                                                          agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                   >  1:1441542‑M2/13‑76 (MQ=255)
                                                          agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                   <  2:1122652‑M2/64‑1 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  2:782457‑M2/65‑1 (MQ=255)
                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 >  2:1512268‑M2/11‑76 (MQ=255)
                                                               aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat              >  2:1128145‑M2/1‑13 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  2:914393‑M2/70‑1 (MQ=255)
                                                                ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             <  2:417713‑M2/76‑65 (MQ=255)
                                                                 atgggCCGCAACCCTGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            <  2:1451201‑M2/71‑1 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           >  1:796034‑M2/1‑10 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           <  1:1149309‑M2/76‑67 (MQ=255)
                                                                   gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact          >  2:613239‑M2/1‑9 (MQ=255)
                                                                     gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct        >  2:1165684‑M2/1‑7 (MQ=255)
                                                                      ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  2:347091‑M2/1‑6 (MQ=255)
                                                                        gCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcac     >  1:980674‑M2/1‑4 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     >  2:321292/1‑76 (MQ=255)
                                                                         cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    >  1:1233824‑M2/1‑3 (MQ=255)
                                                                         cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa    >  1:996467/1‑76 (MQ=255)
                                                                          aacccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacga   <  2:1086387‑M2/76‑75 (MQ=255)
                                                                          aacccataagcgctaacttaagggttgaggtattaagcatgatatgatttaacgtgccgataaattcctatcacga   <  2:1451168‑M2/76‑75 (MQ=255)
                                                                           acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  2:1001606‑M2/1‑1 (MQ=255)
                                                                           |                                                                           
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602

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