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)56 (0.990) 40/144 0.3 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)62 (1.100) 45/144 0.2 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.

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602
                                                                           |                                                                           
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                             <  2:194891‑M2/6‑1 (MQ=255)
aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                             <  1:1625563/76‑1 (MQ=255)
 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCACAAc                                                                            <  1:1515834‑M2/7‑1 (MQ=255)
  ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                           >  1:1725935‑M2/1‑74 (MQ=255)
  tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                           <  2:1356048‑M2/8‑1 (MQ=255)
    gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         >  1:212213‑M2/67‑76 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  2:741793‑M2/1‑70 (MQ=255)
       aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      <  2:776860‑M2/76‑8 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    >  1:852965‑M2/62‑76 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    >  2:1623796‑M2/62‑76 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    <  1:1467318‑M2/15‑1 (MQ=255)
            gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                 <  2:605409‑M2/18‑1 (MQ=255)
            acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                 >  2:512808‑M2/1‑64 (MQ=255)
             aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                                <  1:767443‑M2/19‑1 (MQ=255)
              ggcaccgaggcgtcgcattcttcagatggttcaccccttaagttagcgcttaggggCCGCACCCCCGGGTgctgct                                                               <  2:287552‑M2/20‑1 (MQ=255)
              ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                               <  1:314144‑M2/20‑1 (MQ=255)
              ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                               <  1:1122547‑M2/20‑1 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            <  2:217228‑M2/76‑18 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            <  2:406438‑M2/76‑18 (MQ=255)
                 accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                            <  1:1841331‑M2/23‑1 (MQ=255)
                 accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                            >  2:440223‑M2/54‑76 (MQ=255)
                  ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                           <  2:761467‑M2/24‑1 (MQ=255)
                    gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         >  2:774058‑M2/51‑76 (MQ=255)
                    gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         >  1:1187646‑M2/51‑76 (MQ=255)
                     aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                        <  1:1049729‑M2/27‑1 (MQ=255)
                     aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  2:549718‑M2/76‑22 (MQ=255)
                      ggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa                                                       <  2:686967‑M2/28‑1 (MQ=255)
                       gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:771769‑M2/29‑1 (MQ=255)
                       cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                      <  2:1729576‑M2/76‑24 (MQ=255)
                        cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                     <  1:1089172‑M2/30‑1 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   <  1:979771‑M2/32‑1 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   <  1:1593417‑M2/32‑1 (MQ=255)
                           cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  1:883714‑M2/44‑76 (MQ=255)
                           aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                  >  2:1145951‑M2/1‑49 (MQ=255)
                             cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                                >  1:867771‑M2/1‑47 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                              >  1:1435895‑M2/1‑45 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                              <  1:621675‑M2/76‑32 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                              >  1:366520‑M2/1‑45 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                              >  2:1387909‑M2/1‑45 (MQ=255)
                                 cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                            <  1:1286479‑M2/39‑1 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            <  2:1762500‑M2/76‑34 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            <  1:1133409‑M2/76‑34 (MQ=255)
                                  ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           >  1:1238280‑M2/37‑76 (MQ=255)
                                  ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           >  1:194714‑M2/37‑76 (MQ=255)
                                  ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCAGGTGCTGCTAATGCGGATGCGACCTTCaa                                           >  2:647137‑M2/37‑76 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactttaggggtgtggtattacg                                          <  1:1485746‑M2/76‑36 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          <  1:702704‑M2/76‑36 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  1:516364‑M2/1‑41 (MQ=255)
                                    cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                         >  2:572953‑M2/35‑76 (MQ=255)
                                    cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                         >  2:47710‑M2/35‑76 (MQ=255)
                                      tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                       >  1:1498124‑M2/1‑38 (MQ=255)
                                      tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                       >  1:115644‑M2/1‑38 (MQ=255)
                                      gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  2:1544015‑M2/44‑1 (MQ=255)
                                      gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  2:348056‑M2/44‑1 (MQ=255)
                                      gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  1:249215‑M2/44‑1 (MQ=255)
                                       ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  2:1027362‑M2/1‑37 (MQ=255)
                                       ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  2:207205‑M2/1‑37 (MQ=255)
                                         ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                    >  1:1629555‑M2/30‑76 (MQ=255)
                                          ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  1:672032‑M2/1‑34 (MQ=255)
                                             cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                                >  1:873328‑M2/1‑31 (MQ=255)
                                             cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                                <  2:742594‑M2/76‑46 (MQ=255)
                                              aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               <  2:1545734‑M2/52‑1 (MQ=255)
                                              aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               <  2:159812‑M2/52‑1 (MQ=255)
                                                ncgcgtTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             <  1:561724‑M2/75‑49 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             >  2:1275193‑M2/1‑28 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             >  1:138261‑M2/1‑28 (MQ=255)
                                                 cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  2:1041748‑M2/1‑27 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  1:52056‑M2/55‑1 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  1:1062103‑M2/55‑1 (MQ=255)
                                                  cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           <  2:708664‑M2/56‑1 (MQ=255)
                                                    taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         >  1:1291822‑M2/19‑76 (MQ=255)
                                                    taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         <  1:87369‑M2/58‑1 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  2:818927‑M2/18‑76 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:138547‑M2/1‑22 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  1:148146‑M2/1‑22 (MQ=255)
                                                      agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:1108686‑M2/60‑1 (MQ=255)
                                                       gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      <  1:729862‑M2/61‑1 (MQ=255)
                                                        ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                     <  1:923330‑M2/76‑57 (MQ=255)
                                                        ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                     <  1:706009‑M2/76‑57 (MQ=255)
                                                         tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCTTGGTTCAGTATCAATAACAtactc                    >  2:1153687‑M2/14‑74 (MQ=255)
                                                         tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                    <  1:1794006‑M2/63‑1 (MQ=255)
                                                          ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccga                   >  2:1590879‑M2/1‑18 (MQ=255)
                                                           tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  <  1:676696‑M2/76‑60 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  2:841839‑M2/65‑1 (MQ=255)
                                                            gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 >  2:1082370‑M2/1‑16 (MQ=255)
                                                             cAACCGGGGCCGCAAcccatangcgctaacttaagggtggtggtattacgcctgatatgatttaacgtgccgatga                >  2:293029‑M2/1‑15 (MQ=255)
                                                             cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga                >  2:1548779‑M2/1‑15 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  2:673213‑M2/7‑76 (MQ=255)
                                                                ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             <  1:905836‑M2/76‑65 (MQ=255)
                                                                  tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg           >  2:1084043‑M2/5‑76 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           <  1:926757‑M2/76‑67 (MQ=255)
                                                                   gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  2:622513‑M2/73‑1 (MQ=255)
                                                                   gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  1:237551‑M2/73‑1 (MQ=255)
                                                                    ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  1:609197‑M2/74‑1 (MQ=255)
                                                                    ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  2:1653435‑M2/74‑1 (MQ=255)
                                                                     gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct        >  1:634829‑M2/1‑7 (MQ=255)
                                                                        gCAAcccataagcgctaacttaagggttgtggtattatgcctgatatgacttacngtgccgatgaattactctcac     <  2:294342‑M2/76‑73 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:688932/1‑76 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:478497/1‑76 (MQ=255)
                                                                           acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  <  1:719556‑M2/76‑76 (MQ=255)
                                                                           |                                                                           
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602

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