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)66 (1.080) 52/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)61 (0.990) 50/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 133 total aligned reads displayed.

ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597
                                                                          |                                                                           
gaggcatcgtcttaacgaggctccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            <  1:834654‑M2/2‑1 (MQ=255)
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  2:1575340‑M2/75‑76 (MQ=255)
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  2:1291977‑M2/75‑76 (MQ=255)
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            <  2:103700‑M2/2‑1 (MQ=255)
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  2:597016‑M2/75‑76 (MQ=255)
aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                            <  2:1573556/76‑1 (MQ=255)
    catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                        <  2:926044‑M2/6‑1 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  2:834871/1‑76 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  2:1852937/1‑76 (MQ=255)
     atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       >  1:1136053‑M2/70‑76 (MQ=255)
      ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                      >  2:418999‑M2/1‑74 (MQ=255)
      tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                      >  2:328457‑M2/69‑76 (MQ=255)
      tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                      <  2:1145774‑M2/8‑1 (MQ=255)
        cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    >  1:306518‑M2/1‑72 (MQ=255)
         aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                   >  1:1599251‑M2/1‑71 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  <  1:393085‑M2/76‑7 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  <  1:39485‑M2/76‑7 (MQ=255)
          cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  <  2:1510005‑M2/12‑1 (MQ=255)
           ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCggg                                                                 <  1:991580‑M2/13‑1 (MQ=255)
            taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCcgggg                                                                >  1:975713‑M2/63‑75 (MQ=255)
               aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             <  2:1804865‑M2/76‑12 (MQ=255)
               aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             >  2:1464458‑M2/1‑65 (MQ=255)
                   gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                         >  2:1960058‑M2/56‑76 (MQ=255)
                    caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        <  1:1988959‑M2/22‑1 (MQ=255)
                     gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCCAcccataaacgccaactt                                                       <  1:1362359‑M2/76‑18 (MQ=255)
                     accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       >  2:1910234‑M2/54‑76 (MQ=255)
                     accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       <  1:818697‑M2/23‑1 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      <  1:1555288‑M2/76‑19 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      <  2:264035‑M2/76‑19 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  1:834871‑M2/1‑58 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      <  2:749549‑M2/76‑19 (MQ=255)
                       cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     <  1:1248479‑M2/25‑1 (MQ=255)
                        tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactcaag                                                    >  1:1016352‑M2/1‑56 (MQ=255)
                        gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                    >  1:408545‑M2/51‑76 (MQ=255)
                        gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                    >  2:173551‑M2/51‑76 (MQ=255)
                        aaacGGAACCAATCCCATTCCGCTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                    >  1:1016331‑M2/2‑56 (MQ=255)
                         aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                   >  1:1206171‑M2/50‑76 (MQ=255)
                            ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                <  2:908942‑M2/76‑25 (MQ=255)
                            aaaaccAATACCATGCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                >  2:386037‑M2/3‑52 (MQ=255)
                              aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                              <  1:59803‑M2/76‑27 (MQ=255)
                              aaCCAATACCATTCCGTTCCAGGCGAGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                              <  2:393675‑M2/76‑27 (MQ=255)
                               cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                             >  1:502677‑M2/44‑76 (MQ=255)
                                 cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                           <  2:88632‑M2/35‑1 (MQ=255)
                                 cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                           >  2:850591‑M2/42‑76 (MQ=255)
                                 cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                           >  1:118772‑M2/1‑47 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          >  1:1711843‑M2/41‑76 (MQ=255)
                                  aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          >  2:340083‑M2/1‑46 (MQ=255)
                                  aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          >  2:1314539‑M2/1‑46 (MQ=255)
                                   aTCCCAGTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                         >  2:666251‑M2/1‑45 (MQ=255)
                                    tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                        <  1:47902‑M2/38‑1 (MQ=255)
                                     cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       >  2:1417411‑M2/38‑76 (MQ=255)
                                     aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                       >  2:1586549‑M2/1‑43 (MQ=255)
                                      ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                      <  1:1598521‑M2/40‑1 (MQ=255)
                                        cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                    >  2:1930474‑M2/35‑76 (MQ=255)
                                           ccGTTCCAGGCGCGTTATTTCGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 >  2:1556247‑M2/1‑37 (MQ=255)
                                           atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                 >  1:1961967‑M2/32‑76 (MQ=255)
                                            cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                <  1:728877‑M2/76‑41 (MQ=255)
                                              ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  1:1503059‑M2/76‑43 (MQ=255)
                                              ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  2:1642295‑M2/1‑34 (MQ=255)
                                              ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  1:9304‑M2/76‑43 (MQ=255)
                                              ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggggtgtggtattacgcctggta                              <  1:760699‑M2/76‑43 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              >  1:1575652‑M2/29‑76 (MQ=255)
                                               ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             >  1:658775‑M2/28‑76 (MQ=255)
                                               ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             <  2:673631‑M2/49‑1 (MQ=255)
                                               ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             <  2:591972‑M2/49‑1 (MQ=255)
                                                tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            <  1:1764327‑M2/50‑1 (MQ=255)
                                                ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            >  1:1852937‑M2/1‑32 (MQ=255)
                                                ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            <  2:1116411‑M2/76‑45 (MQ=255)
                                                ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            <  1:1627251‑M2/76‑45 (MQ=255)
                                                  aacccttaagttagtgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                          <  1:654877‑M2/52‑1 (MQ=255)
                                                    gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        <  2:1161887‑M2/76‑49 (MQ=255)
                                                    gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  2:1051409‑M2/1‑28 (MQ=255)
                                                     cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       <  1:2022903‑M2/76‑50 (MQ=255)
                                                     cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:1112749‑M2/1‑27 (MQ=255)
                                                      gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  2:2024883‑M2/1‑26 (MQ=255)
                                                      gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  2:1651851‑M2/1‑26 (MQ=255)
                                                       ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     >  2:1071364‑M2/20‑76 (MQ=255)
                                                       cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     <  1:1261738‑M2/76‑52 (MQ=255)
                                                          tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  >  1:178213‑M2/1‑22 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  >  2:1625764‑M2/17‑76 (MQ=255)
                                                           gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 >  1:48344‑M2/16‑76 (MQ=255)
                                                            ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                <  1:309158‑M2/76‑57 (MQ=255)
                                                            ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                >  1:745284‑M2/15‑76 (MQ=255)
                                                              agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa              >  2:1867092‑M2/13‑76 (MQ=255)
                                                               tGCAGCCGGGGCCGCAGcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat             >  2:814504‑M2/1‑17 (MQ=255)
                                                               gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             <  2:1752058‑M2/65‑1 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:605519‑M2/1‑15 (MQ=255)
                                                                  cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          >  1:1056287‑M2/9‑76 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:1334516‑M2/8‑76 (MQ=255)
                                                                    tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        <  2:298782‑M2/70‑1 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        <  1:17566‑M2/76‑65 (MQ=255)
                                                                     cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  2:1740804‑M2/1‑11 (MQ=255)
                                                                      tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      >  2:1129974‑M2/5‑76 (MQ=255)
                                                                      ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac      <  1:273105‑M2/76‑67 (MQ=255)
                                                                       gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact     >  1:232148‑M2/1‑9 (MQ=255)
                                                                        ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    >  2:1108883‑M2/1‑8 (MQ=255)
                                                                        ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCagggg    <  2:2023167‑M2/74‑2 (MQ=255)
                                                                         gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc   <  1:1504271‑M2/75‑1 (MQ=255)
                                                                          ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  <  1:1869198‑M2/76‑71 (MQ=255)
                                                                          ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg  >  1:21385/1‑76 (MQ=255)
                                                                          |                                                                           
ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597

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.