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)264 (1.160) 108/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)240 (1.060) 107/144 0.1 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 520 total aligned reads displayed.

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596
                                                                           |                                                                          
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                                                            >  2:169655/1‑76 (MQ=255)
 gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                           >  2:4548251‑M2/75‑76 (MQ=255)
 aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                           >  2:2456459/1‑76 (MQ=255)
 aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                           >  2:1128046/1‑76 (MQ=255)
   ggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGc                                                                         >  1:3984004‑M2/73‑76 (MQ=255)
    tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                        <  1:4576874/76‑1 (MQ=255)
     aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                       >  2:6675607/1‑76 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                      >  1:1804078‑M2/70‑76 (MQ=255)
       tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                     <  2:6400244‑M2/8‑1 (MQ=255)
        cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                    >  1:3874512‑M2/68‑76 (MQ=255)
         gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                   <  1:6768312‑M2/10‑1 (MQ=255)
         gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                   >  2:3686737‑M2/67‑76 (MQ=255)
          tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                  >  2:3456442‑M2/66‑76 (MQ=255)
          aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                  >  1:4101896‑M2/1‑71 (MQ=255)
              acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                              <  1:3545613‑M2/76‑10 (MQ=255)
              aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                              >  1:3384977‑M2/62‑76 (MQ=255)
                  cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                          <  2:2760640‑M2/76‑14 (MQ=255)
                  aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                          <  1:5297025‑M2/19‑1 (MQ=255)
                    cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                        >  2:1843246‑M2/1‑61 (MQ=255)
                    cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                        >  1:485449‑M2/1‑61 (MQ=255)
                    cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                        <  1:5832331‑M2/76‑16 (MQ=255)
                     tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                       >  1:2665442‑M2/1‑60 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                      <  1:3227775‑M2/76‑18 (MQ=255)
                      accgaggtgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                      >  1:2470336‑M2/54‑76 (MQ=255)
                      accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                      >  1:1347296‑M2/54‑76 (MQ=255)
                      accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                      <  2:7120229‑M2/23‑1 (MQ=255)
                       aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                     >  2:6104660‑M2/1‑58 (MQ=255)
                       aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                     <  2:4321817‑M2/76‑19 (MQ=255)
                        aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                    <  2:6989256‑M2/76‑20 (MQ=255)
                         tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                   >  1:2394243‑M2/1‑56 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  2:2549264‑M2/76‑22 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                <  1:5392240‑M2/29‑1 (MQ=255)
                            cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAActcataagcgctaacttaagggt                                                <  2:5727058‑M2/76‑24 (MQ=255)
                             cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                               <  1:1022269‑M2/30‑1 (MQ=255)
                              gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                              >  2:3214810‑M2/1‑51 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             >  1:5224347‑M2/1‑50 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             >  1:5487490‑M2/1‑50 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             >  1:991785‑M2/1‑50 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             <  2:3788046‑M2/76‑27 (MQ=255)
                                 gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                           >  2:2619524‑M2/43‑76 (MQ=255)
                                 gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                           >  1:1926992‑M2/43‑76 (MQ=255)
                                 gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                           <  1:6105360‑M2/34‑1 (MQ=255)
                                   attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                         >  2:5097767‑M2/41‑76 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                         >  1:4243526‑M2/1‑46 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                         <  2:5971284‑M2/76‑31 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                         <  1:4677281‑M2/76‑31 (MQ=255)
                                    ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                        >  1:7147156‑M2/40‑76 (MQ=255)
                                    attccattccGTTCCAGGCGCGTTATTTGGAAACCGGGGCAGCAAcccataagcgctaacttaagggttgtggtat                                        >  1:1250878‑M2/1‑45 (MQ=255)
                                      aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                      >  1:503428‑M2/1‑43 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                     >  2:61788‑M2/1‑42 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                     >  1:2300195‑M2/1‑42 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                    <  2:394566‑M2/76‑36 (MQ=255)
                                         aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                   >  2:758552‑M2/1‑40 (MQ=255)
                                         aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                   <  2:1467382‑M2/76‑37 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                 <  2:5785346‑M2/76‑39 (MQ=255)
                                            ccGTTCGAGGCACGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                >  2:3359947‑M2/1‑37 (MQ=255)
                                             tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                               <  1:2587821‑M2/46‑1 (MQ=255)
                                             cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                               <  2:7269341‑M2/76‑41 (MQ=255)
                                              gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                              <  1:2192298‑M2/76‑42 (MQ=255)
                                              gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                              >  1:309051‑M2/1‑35 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                             >  2:4416176‑M2/1‑34 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                             <  1:372201‑M2/76‑43 (MQ=255)
                                               gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                             >  1:7073742‑M2/29‑76 (MQ=255)
                                                ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                            <  1:3460519‑M2/49‑1 (MQ=255)
                                                tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                            >  1:6843807‑M2/1‑33 (MQ=255)
                                                tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                            >  2:1281503‑M2/1‑33 (MQ=255)
                                                 tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                           <  1:7291928‑M2/50‑1 (MQ=255)
                                                   aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                         <  2:3003349‑M2/76‑47 (MQ=255)
                                                   aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                         >  1:5136400‑M2/1‑30 (MQ=255)
                                                     gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                       >  2:3131775‑M2/1‑28 (MQ=255)
                                                     gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                       >  2:4722320‑M2/1‑28 (MQ=255)
                                                     cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                       <  2:1200749‑M2/54‑1 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                      >  1:275605‑M2/1‑27 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                     >  2:1610271‑M2/21‑76 (MQ=255)
                                                          aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                  <  2:5395957‑M2/59‑1 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 >  1:5735223‑M2/17‑76 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 <  2:6493413‑M2/60‑1 (MQ=255)
                                                            aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                <  2:3583667‑M2/76‑56 (MQ=255)
                                                              tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg              <  1:6502377‑M2/76‑58 (MQ=255)
                                                              tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc              >  1:6382212‑M2/14‑76 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat            >  1:786890‑M2/1‑17 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg           >  1:6899372‑M2/1‑16 (MQ=255)
                                                                 cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg           <  1:4398665‑M2/66‑1 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga          >  1:6180412‑M2/1‑15 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga          >  2:5543990‑M2/1‑15 (MQ=255)
                                                                   cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt         <  2:781365‑M2/68‑1 (MQ=255)
                                                                   cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt         >  1:1198784‑M2/9‑76 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc        >  2:49047‑M2/8‑76 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc        >  1:5896880‑M2/8‑76 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt       >  1:7432667‑M2/1‑12 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt       >  1:1674366‑M2/1‑12 (MQ=255)
                                                                      cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta      <  2:5301767‑M2/76‑66 (MQ=255)
                                                                      atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg      >  2:2844804‑M2/6‑76 (MQ=255)
                                                                      atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg      >  2:2153280‑M2/6‑76 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     <  2:4796748‑M2/72‑1 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     >  2:2399266‑M2/5‑76 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     >  2:1340864‑M2/5‑76 (MQ=255)
                                                                       ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac     >  1:146277‑M2/1‑10 (MQ=255)
                                                                         ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa   <  2:2245826‑M2/74‑1 (MQ=255)
                                                                          gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct  >  2:4139377‑M2/1‑7 (MQ=255)
                                                                           |                                                                          
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596

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