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)61 (0.910) 48/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)73 (1.090) 54/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 134 total aligned reads displayed.

GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602
                                                                          |                                                                           
atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                            >  1:2187998‑M2/70‑76 (MQ=255)
atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                            >  1:2156802‑M2/70‑76 (MQ=255)
 tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                           <  2:658565‑M2/8‑1 (MQ=255)
   gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         >  2:1731153‑M2/67‑76 (MQ=255)
   cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         <  2:220097‑M2/76‑5 (MQ=255)
     gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGAAAcccata                                                                       >  1:1588375‑M2/1‑70 (MQ=255)
      aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      <  1:2037626‑M2/76‑8 (MQ=255)
       aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     <  2:338004‑M2/76‑9 (MQ=255)
         caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                   <  1:1395402‑M2/76‑11 (MQ=255)
            cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                <  1:1235695‑M2/76‑14 (MQ=255)
            cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                <  1:75314‑M2/76‑14 (MQ=255)
             ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               <  1:1866220‑M2/76‑15 (MQ=255)
              gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                              <  2:1019147‑M2/21‑1 (MQ=255)
              cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                              <  1:292292‑M2/76‑16 (MQ=255)
                gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            >  1:380939‑M2/1‑59 (MQ=255)
                accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                            <  1:2010962‑M2/23‑1 (MQ=255)
                  aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                          <  1:1156972‑M2/76‑20 (MQ=255)
                   tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                         <  1:304751‑M2/76‑21 (MQ=255)
                   gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         <  1:2117902‑M2/26‑1 (MQ=255)
                    aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                        >  2:1273151‑M2/50‑76 (MQ=255)
                    aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccacaagcgctaacttaagg                                                        <  1:1647700‑M2/76‑22 (MQ=255)
                     aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                       <  1:1341216‑M2/76‑23 (MQ=255)
                       ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     >  2:385894‑M2/1‑52 (MQ=255)
                       cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                     <  2:1734451‑M2/30‑1 (MQ=255)
                        gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                    >  1:1057524‑M2/1‑51 (MQ=255)
                          cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  2:1677787‑M2/44‑76 (MQ=255)
                          aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                  >  1:122930‑M2/1‑49 (MQ=255)
                           gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                 <  1:1478470‑M2/34‑1 (MQ=255)
                            cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                >  1:611963‑M2/42‑76 (MQ=255)
                            cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                <  2:254977‑M2/35‑1 (MQ=255)
                             attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                               <  1:362043‑M2/36‑1 (MQ=255)
                             actacCATTCCGTTCCGGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               >  2:443908‑M2/3‑46 (MQ=255)
                             aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               <  2:612464‑M2/76‑31 (MQ=255)
                              ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                              >  1:1406923‑M2/40‑76 (MQ=255)
                               tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                             >  1:2196009‑M2/39‑76 (MQ=255)
                                tccaTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:2210659‑M2/2‑43 (MQ=255)
                                cttcagatcgttcagcccttaaattagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                            <  1:1052585‑M2/39‑1 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  1:101884‑M2/1‑43 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  1:2208234‑M2/1‑43 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:645630‑M2/1‑43 (MQ=255)
                                 ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           <  1:960512‑M2/40‑1 (MQ=255)
                                  tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                          <  1:1978365‑M2/41‑1 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  1:132617‑M2/1‑41 (MQ=255)
                                   cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                         >  2:1109216‑M2/35‑76 (MQ=255)
                                    agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                        <  2:1675933‑M2/43‑1 (MQ=255)
                                      ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      <  1:1283818‑M2/76‑40 (MQ=255)
                                        gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggggtgtggtgttatggctgat                                    <  1:458057‑M2/76‑42 (MQ=255)
                                         ttCCAGGCGCGTTATTTTGCAACCTGCGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  2:1537869‑M2/1‑34 (MQ=255)
                                          tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  >  2:366762‑M2/1‑33 (MQ=255)
                                           tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                 <  1:169809‑M2/50‑1 (MQ=255)
                                           ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 <  1:1706336‑M2/76‑45 (MQ=255)
                                             aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               <  2:1474623‑M2/52‑1 (MQ=255)
                                             aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                               <  2:2194442‑M2/52‑1 (MQ=255)
                                              ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatt                              >  2:1652346‑M2/1‑29 (MQ=255)
                                               gcgcGTTATTTTGCAGCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             >  2:1753977‑M2/1‑28 (MQ=255)
                                               gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                             >  2:1552384‑M2/1‑28 (MQ=255)
                                               cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                             >  2:124271‑M2/23‑76 (MQ=255)
                                               cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGGTCAGTATCAAt                             >  1:1789765‑M2/23‑76 (MQ=255)
                                                ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            >  1:551896‑M2/22‑76 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           <  2:1182591‑M2/56‑1 (MQ=255)
                                                   taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         >  2:400280‑M2/19‑76 (MQ=255)
                                                   taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         >  2:902003‑M2/19‑76 (MQ=255)
                                                    ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        >  2:287715‑M2/1‑23 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:443148‑M2/1‑22 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       <  1:979156‑M2/76‑55 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:2150155‑M2/60‑1 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       >  2:2070508‑M2/17‑76 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       >  2:478530‑M2/17‑76 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:1502241‑M2/60‑1 (MQ=255)
                                                      gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      <  1:1546788‑M2/61‑1 (MQ=255)
                                                      gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      <  1:935171‑M2/61‑1 (MQ=255)
                                                       ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                     <  1:860398‑M2/62‑1 (MQ=255)
                                                          tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  >  1:1771465‑M2/1‑17 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 <  1:729643‑M2/76‑61 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 <  2:236056‑M2/76‑61 (MQ=255)
                                                            gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                >  2:1606226‑M2/10‑76 (MQ=255)
                                                            gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                >  1:688035‑M2/10‑76 (MQ=255)
                                                             aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa               <  1:1911724‑M2/76‑63 (MQ=255)
                                                              aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat              <  2:1349853‑M2/76‑64 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  1:983000‑M2/7‑76 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  1:1326072‑M2/70‑1 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  1:1168166‑M2/70‑1 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  2:864989‑M2/70‑1 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  1:2045160‑M2/7‑76 (MQ=255)
                                                               ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             >  2:1417609‑M2/1‑12 (MQ=255)
                                                                 tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCCGGTTCAgg           >  2:397808‑M2/5‑76 (MQ=255)
                                                                 ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           >  1:935799‑M2/1‑10 (MQ=255)
                                                                  ggcccGCACCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  1:793588/73‑1 (MQ=255)
                                                                      cGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctca      <  1:1717967‑M2/76‑72 (MQ=255)
                                                                      cGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctca      >  1:549965‑M2/1‑5 (MQ=255)
                                                                      cGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctca      >  2:1530647‑M2/1‑5 (MQ=255)
                                                                       gCAAcccataagcgctaacttaagggttgtggtattacgcatgatatgatttaacgtgccgatgaattactctcac     <  2:1762968‑M2/76‑73 (MQ=255)
                                                                        cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    <  2:2041750‑M2/76‑74 (MQ=255)
                                                                        cAAcccataagcgctaacttaagggttgtggaattacgcctgatatgatttaacgtgccgatgaattactctcacg    <  1:1154849‑M2/76‑74 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:252092/1‑76 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  1:265810/1‑76 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:1811913/1‑76 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:169572/1‑76 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:1267858/1‑76 (MQ=255)
                                                                          aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  <  2:1138507/76‑1 (MQ=255)
                                                                          |                                                                           
GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602

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