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 annotation gene product
* ? REL606 15386 =NA (NA)54 (1.080) 35/64 0.1 noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 57 (1.040)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)6 (0.120) 31/64 0.2 noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 57 (1.040)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 108 total aligned reads displayed.

CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT  >  REL606/4524487‑4524557
                                   |                                   
  agatggttcaacccttaagttagcgcttatgggCCg                                   >  1:4039992‑M2/34‑36 (MQ=255)
   tctgTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                  <  1:7461718/33‑1 (MQ=255)
   tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                  >  1:5902461/1‑36 (MQ=255)
   tCCGTCCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                  <  1:3200168/36‑1 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                  >  1:1265936‑M2/33‑36 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                  >  1:5824220‑M2/33‑36 (MQ=255)
    ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGcc                                 >  1:6770324/1‑35 (MQ=255)
    atggttcaacccttaagttagcgcttatgggCCGCa                                 <  1:4266878‑M2/5‑1 (MQ=255)
     tggttcaacccttaagttagcgcttatgggCCGCaa                                <  1:1071643‑M2/6‑1 (MQ=255)
     tggttcaacccttaagttagcgcttatgggCCGCaa                                <  1:4521772‑M2/6‑1 (MQ=255)
     cGTTCCGGGCGCGTTATTTTGCAACCGGGGCCGCaa                                <  1:6350244/36‑1 (MQ=255)
     cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                >  1:5139506/1‑36 (MQ=255)
      gTTCCAGGCGCGTTATTTTGTAACCGGGGCGGGAAc                               >  1:3095831/1‑36 (MQ=21)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                               >  1:6838290‑M2/1‑35 (MQ=255)
       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                              >  1:3861763‑M2/1‑34 (MQ=255)
       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                              <  1:7448394‑M2/36‑3 (MQ=255)
       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                              <  1:4399175‑M2/36‑3 (MQ=255)
       ttCCAAGCGCGTTTTTTTGCAACCGGGGCAGCAAcc                              <  1:914940/36‑1 (MQ=21)
       gttcaacccttaagttagcgcttatgggCCGCAAcc                              <  1:1343633‑M2/8‑1 (MQ=255)
        ttcaacccttaagttagcgcttatgggCCGCACccc                             >  1:4178553‑M2/28‑36 (MQ=255)
        ttcaacccttaagttagcgcttatgggCCGCAAccc                             >  1:2341753‑M2/28‑36 (MQ=255)
        tCCAGGCGCGTTATTTTGCAACCGGGGCCGCACccc                             >  1:5546262‑M2/1‑33 (MQ=255)
         tcaacccttaagttagcgcttatgggCCGCAAcccc                            <  1:2463458‑M2/10‑1 (MQ=255)
         ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                            >  1:642942‑M2/1‑32 (MQ=255)
          cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                           <  1:5187216‑M2/36‑6 (MQ=255)
           aacccttaagttagcgcttatgggCCGCAACCCCgg                          >  1:4860245‑M2/25‑36 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          >  1:2145058‑M2/1‑30 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          <  1:4810629‑M2/36‑7 (MQ=255)
            acccttaagttagcgcttatgggCCGCAACCCCggg                         <  1:4250074‑M2/13‑1 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGTCGCAAcccataag                        >  1:5583018‑M2/1‑28 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                        <  1:7377749‑M2/36‑9 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        >  1:3619508‑M2/23‑36 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        <  1:4958451‑M2/14‑1 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        <  1:4101367‑M2/14‑1 (MQ=255)
              cgcgTTATTTTGCAACCGGGGCCGCAAcccataagc                       >  1:2834966‑M2/1‑27 (MQ=255)
              ccttaagttagcgcttaagggCCGCAACCCCGGGtg                       <  1:7358450‑M2/15‑1 (MQ=255)
               gcgTTATTTTGCAACCGGGGCCGCAAcccataagcg                      >  1:81132‑M2/1‑26 (MQ=255)
               gcgTTATTTTGCAACCGGGGCCGCAAcccataagcg                      <  1:5180171‑M2/36‑11 (MQ=255)
               cttaagttagcgcttatgggCCGCAACCCCGGGtgc                      <  1:6905572‑M2/16‑1 (MQ=255)
               cttaagttagcgcttatgggCCGCAACCCCGGGtgc                      <  1:2929066‑M2/16‑1 (MQ=255)
               cttaagttagcgcttataggCCGCAACCCCGGGtgc                      <  1:4717369‑M2/16‑1 (MQ=255)
                cgTTATTTTGCAACCGCGGCCGCAAcccataagcgc                     <  1:1935812‑M2/36‑12 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                    <  1:549587‑M2/18‑1 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                    >  1:4935571‑M2/19‑36 (MQ=255)
                 taagttaccgcttatgggCCGCAACCCCGGGtgctg                    <  1:7066878‑M2/18‑1 (MQ=255)
                 gTTATTTTGCAACCGGGGCCGCAAcccataagcgct                    >  1:1405983‑M2/1‑24 (MQ=255)
                  ttATTTTGCAACCGGGGGCGCAAcccataagcgcta                   >  1:2464260‑M2/1‑23 (MQ=255)
                  ttATTTTGCAACCGGGGCCGCAAcccataagcgcta                   >  1:7415322‑M2/1‑23 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                   <  1:258081‑M2/19‑1 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  <  1:6488611‑M2/36‑15 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  >  1:2248073‑M2/1‑22 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  <  1:2864217‑M2/36‑15 (MQ=255)
                   tATTTAGCAACCGGGGCCGCAAcccataagcgctaa                  >  1:6670792‑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               >  1:6136298‑M2/14‑36 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              <  1:7222709‑M2/36‑19 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              >  1:943054‑M2/1‑18 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              <  1:321364‑M2/36‑19 (MQ=255)
                        tGCAACCGGGGCCGCAAcccataagcgctaacttaa             <  1:3641616‑M2/36‑20 (MQ=255)
                        gcgctaatgggCCGCAACCCCGGGTGCTGCTAATGc             <  1:6019241‑M2/25‑1 (MQ=255)
                         gacaccgGGGCCGCAAcccataagcgctaacttaag            >  1:6903195‑M2/4‑16 (MQ=255)
                         gCAACCGGGGCCGCAAcccataagcgctaactttag            <  1:3923375‑M2/36‑21 (MQ=255)
                         gCAACCGGGGCCGCAAcccataagcgctaacttaag            >  1:2097099‑M2/1‑16 (MQ=255)
                         cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg            >  1:4615646‑M2/11‑36 (MQ=255)
                         cgcttatgggCCGCAACCCAGGGTGCTGCTAATGCg            >  1:6733418‑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          >  1:5347883‑M2/9‑36 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:2048939‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:700731‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          >  1:1489061‑M2/1‑14 (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         <  1:6909647‑M2/29‑1 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         >  1:5216527‑M2/8‑36 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         >  1:5482594‑M2/8‑36 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt         >  1:5102061‑M2/8‑36 (MQ=255)
                            aCCGGGGCCGCAAcccataagcgctaacttaagggt         >  1:4212507‑M2/1‑13 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg        <  1:878675‑M2/30‑1 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg        <  1:3907342‑M2/30‑1 (MQ=255)
                               ggggCCGCAAcccataagcgctaacttaagggttgt      >  1:434581‑M2/1‑10 (MQ=255)
                               ggggCCGAAAcccataagcgctaacttaagggttgt      >  1:6947665‑M2/1‑10 (MQ=255)
                                gggCCGCAAcccataagcgctaacttaagggttgtg     >  1:4058506‑M2/1‑9 (MQ=255)
                                gggCCGCAAcccataagcgctaacttaagggttgtg     >  1:136726‑M2/1‑9 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa     >  1:7035624‑M2/4‑36 (MQ=255)
                                 ggCCGCAAcccataagcgctaacttaagggttgtgg    >  1:1441030‑M2/1‑8 (MQ=255)
                                 ggCCGCAACCCCGGGTGCTGCTAATGCGGCTGGGAc    >  1:3056011/1‑36 (MQ=25)
                                 ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc    <  1:3818081‑M2/34‑1 (MQ=255)
                                 gcccgcAAcccataagcgctaacttaagggttgtgg    >  1:2749072‑M2/3‑8 (MQ=255)
                                  gCCGCAAcccataagcgctaacttaagggttgtggt   <  1:5389097‑M2/36‑30 (MQ=255)
                                  gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc   >  1:1256085‑M2/2‑36 (MQ=255)
                                   ccGCAAcccataagcgctaacttaagggttgtggta  >  1:7377165‑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