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)54 (1.070) 35/64 0.1 94.7% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 3 (0.220)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)45 (0.890) 31/64 0.2 95.7% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 2 (0.160)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 108 total aligned reads displayed.

CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT  >  REL606/4524487‑4524557
                                   |                                   
cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                     <  2:60182/36‑1 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                   >  2:300782‑M2/34‑36 (MQ=255)
   tctgTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                  <  2:3722508/33‑1 (MQ=255)
   tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                  >  2:2163251/1‑36 (MQ=255)
   tCCGTCCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                  <  1:3200168/36‑1 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                  >  1:1265936‑M2/33‑36 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                  >  2:2085010‑M2/33‑36 (MQ=255)
    atggttcaacccttaagttagcgcttatgggCCGCa                                 <  2:527668‑M2/5‑1 (MQ=255)
     tggttcaacccttaagttagcgcttatgggCCGCaa                                <  1:1071643‑M2/6‑1 (MQ=255)
     tggttcaacccttaagttagcgcttatgggCCGCaa                                <  2:782562‑M2/6‑1 (MQ=255)
     cGTTCCGGGCGCGTTATTTTGCAACCGGGGCCGCaa                                <  2:2611034/36‑1 (MQ=255)
     cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                >  2:1400296/1‑36 (MQ=255)
      gTTCCAGGCGCGTTATTTTGTAACCGGGGCGGGAAc                               >  1:3095831/1‑36 (MQ=25)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                               >  2:3099080‑M2/1‑35 (MQ=255)
       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                              <  2:3709184‑M2/36‑3 (MQ=255)
       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                              >  2:122553‑M2/1‑34 (MQ=255)
       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                              <  2:659965‑M2/36‑3 (MQ=255)
       ttCCAAGCGCGTTTTTTTGCAACCGGGGCAGCAAcc                              <  1:914940/36‑1 (MQ=25)
       gttcaacccttaagttagcgcttatgggCCGCATcc                              >  2:2746390‑M2/29‑36 (MQ=255)
       gttcaacccttaagttagcgcttatgggCCGCAAcc                              <  1:1343633‑M2/8‑1 (MQ=255)
        ttcaacccttaagttagcgcttatgggCCGCACccc                             >  2:439343‑M2/28‑36 (MQ=255)
        ttcaacccttaagttagcgcttatgggCCGCAAccc                             >  1:2341753‑M2/28‑36 (MQ=255)
         tcaacccttaagttagcgcttatgggCCGCAAcccc                            <  1:2463458‑M2/10‑1 (MQ=255)
         ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                            >  1:642942‑M2/1‑32 (MQ=255)
          cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                           <  2:1448006‑M2/36‑6 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          >  1:2145058‑M2/1‑30 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          <  2:1071419‑M2/36‑7 (MQ=255)
            acccttaagttagcgcttatgggCCGCAACCCCggg                         <  2:510864‑M2/13‑1 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGTCGCAAcccataag                        >  2:1843808‑M2/1‑28 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                        <  2:3638539‑M2/36‑9 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        <  2:1219241‑M2/14‑1 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        <  2:362157‑M2/14‑1 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        >  1:3619508‑M2/23‑36 (MQ=255)
              cgcgTTATTTTGCAACCGGGGCCGCAAcccataagc                       >  1:2834966‑M2/1‑27 (MQ=255)
              ccttaagttagcgcttaagggCCGCAACCCCGGGtg                       <  2:3619240‑M2/15‑1 (MQ=255)
               gcgTTATTTTGCAACCGGGGCCGCAAcccataagcg                      <  2:1440961‑M2/36‑11 (MQ=255)
               gcgTTATTTTGCAACCGGGGCCGCAAcccataagcg                      >  1:81132‑M2/1‑26 (MQ=255)
               gcgTTATTTTGCAACCGGGGCCGCAAcccataagcg                      >  1:981573‑M2/1‑26 (MQ=255)
               cttaagttagcgcttatgggCCGCAACCCCGGGtgc                      <  2:3166362‑M2/16‑1 (MQ=255)
               cttaagttagcgcttatgggCCGCAACCCCGGGtgc                      <  1:2929066‑M2/16‑1 (MQ=255)
                cgTTATTTTGCAACCGCGGCCGCAAcccataagcgc                     <  1:1935812‑M2/36‑12 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                    <  1:549587‑M2/18‑1 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                    >  2:1196361‑M2/19‑36 (MQ=255)
                 taagttaccgcttatgggCCGCAACCCCGGGtgctg                    <  2:3327668‑M2/18‑1 (MQ=255)
                 taagttaccgcttatgggCCGCAACCCCGGGtgctg                    <  2:1508696‑M2/18‑1 (MQ=255)
                 gTTATTTTGCAACCGGGGCCGCAAcccataagcgct                    >  1:1405983‑M2/1‑24 (MQ=255)
                  ttATTTTGCAACCGGGGGCGCAAcccataagcgcta                   >  1:2464260‑M2/1‑23 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                   <  1:258081‑M2/19‑1 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  >  1:2248073‑M2/1‑22 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  <  1:2864217‑M2/36‑15 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  <  2:2749401‑M2/36‑15 (MQ=255)
                   tATTTAGCAACCGGGGCCGCAAcccataagcgctaa                  >  2:2931582‑M2/1‑22 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                 <  1:853643‑M2/21‑1 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                 >  1:1620754‑M2/16‑36 (MQ=255)
                     ttttTCAACCGGGGCCGCAAcccataagcgctaact                >  1:1230548‑M2/1‑20 (MQ=255)
                     ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                <  1:287128‑M2/22‑1 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCtaac               >  1:860956‑M2/14‑35 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt               <  1:3706716‑M2/23‑1 (MQ=255)
                      tagcacttatgggCCGCAACCCCGGGTGCTGCTAAt               >  2:2397088‑M2/14‑36 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              >  1:943054‑M2/1‑18 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              <  2:3483499‑M2/36‑19 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              <  1:321364‑M2/36‑19 (MQ=255)
                        tGCAACCGGGGCCGCAAcccataagcgctaacttaa             <  1:3641616‑M2/36‑20 (MQ=255)
                        gcgctaatgggCCGCAACCCCGGGTGCTGCTAATGc             <  2:2280031‑M2/25‑1 (MQ=255)
                         gacaccgGGGCCGCAAcccataagcgctaacttaag            >  2:3163985‑M2/4‑16 (MQ=255)
                         gCAACCGGGGCCGCAAcccataagcgctaactttag            <  2:184165‑M2/36‑21 (MQ=255)
                         gCAACCGGGGCCGCAAcccataagcgctaacttaag            >  1:2097099‑M2/1‑16 (MQ=255)
                         cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg            >  2:876436‑M2/11‑36 (MQ=255)
                         cgcttatgggCCGCAACCCAGGGTGCTGCTAATGCg            >  2:2994208‑M2/11‑36 (MQ=255)
                          gcttatgggCCCCAACCCCGGGTGCTGCTAATGCgg           <  1:170787‑M2/27‑1 (MQ=255)
                          cAACCGGGGCCGCAAcccataagcgctaacttaagg           <  1:1666676‑M2/36‑22 (MQ=255)
                          cAACCGGGGCCGCAAcccataagcgctaacttaagg           <  1:3339580‑M2/36‑22 (MQ=255)
                           cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa          <  2:623249‑M2/28‑1 (MQ=255)
                           cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa          >  2:1608673‑M2/9‑36 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          >  1:1489061‑M2/1‑14 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:2048939‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:73382‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:700731‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:316591‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacctaaggg          <  1:1058174‑M2/36‑23 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         <  1:271752‑M2/29‑1 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         >  2:1743384‑M2/8‑36 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         >  2:1477317‑M2/8‑36 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         >  2:1362851‑M2/8‑36 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         <  2:3170437‑M2/29‑1 (MQ=255)
                            aCCGGGGCCGCAAcccataagcgctaacttaagggt         >  2:473297‑M2/1‑13 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg        <  1:878675‑M2/30‑1 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg        <  2:168132‑M2/30‑1 (MQ=255)
                               ggggCCGAAAcccataagcgctaacttaagggttgt      >  2:3208455‑M2/1‑10 (MQ=255)
                                gggCCGCAAcccataagcgctaacttaagggttgtg     >  2:319296‑M2/1‑9 (MQ=255)
                                gggCCGCAAcccataagcgctaacttaagggttgtg     >  1:136726‑M2/1‑9 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa     >  2:3296414‑M2/4‑36 (MQ=255)
                                 ggCCGCAAcccataagcgctaacttaagggttgtgg    >  1:1441030‑M2/1‑8 (MQ=255)
                                 ggCCGCAACCCCGGGTGCTGCTAATGCGGCTGGGAc    >  1:3056011/1‑36 (MQ=38)
                                 ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc    <  2:78871‑M2/34‑1 (MQ=255)
                                 gcccgcAAcccataagcgctaacttaagggttgtgg    >  1:2749072‑M2/3‑8 (MQ=255)
                                  gCCGCAAcccataagcgctaacttaagggttgtggt   <  2:1649887‑M2/36‑30 (MQ=255)
                                  gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc   >  1:1256085‑M2/2‑36 (MQ=255)
                                   ccGCAAcccataagcgctaacttaagggttgtggta  >  2:3637955‑M2/1‑6 (MQ=255)
                                   ccGCAAcccataagcgctaacttaagggttgtggta  <  1:1860565‑M2/36‑31 (MQ=255)
                                   |                                   
CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT  >  REL606/4524487‑4524557

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG