Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 93.5% 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 15386 =NA (NA)40 (0.710) 29/64 NT 93.0% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 3 (0.150)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)46 (0.810) 31/64 NT 93.9% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 3 (0.160)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 101 total aligned reads displayed.

CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT  >  REL606/4524487‑4524557
                                   |                                   
tcagatggttcaacccttaagttagcgcttatgggc                                     >  1:3924413‑M2/36‑36 (MQ=255)
 cagatggttaaacccttaagttagcgcttatgggcc                                    >  1:2840441‑M2/35‑36 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                    <  1:1195799/36‑1 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                    <  1:4745113/36‑1 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                    <  1:4707265/36‑1 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                    >  1:4699591/1‑36 (MQ=255)
  ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                   >  1:7961822/1‑36 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                   <  1:800679‑M2/3‑1 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                   >  1:1757542‑M2/34‑36 (MQ=255)
     tggttcaacccttaagttagcgcttatgggCCGCaa                                >  1:3606800‑M2/31‑36 (MQ=255)
     cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                <  1:3020991/36‑1 (MQ=255)
     cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                <  1:7556021/36‑1 (MQ=255)
      ggttcaacccttaagttagcgcttatgggCCGCAAc                               >  1:1536771‑M2/30‑36 (MQ=255)
      ggttcaacccttaagttagcgcttatgggCCGCAAc                               >  1:1859660‑M2/30‑36 (MQ=255)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                               <  1:3306910‑M2/36‑2 (MQ=255)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                               <  1:3566822‑M2/36‑2 (MQ=255)
       gttcaacccttaagttagcgcttatgggCCGCAAcc                              >  1:3704833‑M2/29‑36 (MQ=255)
       gttcaacccttaagttagcgcttatgggCCGCAAcc                              >  1:1177990‑M2/29‑36 (MQ=255)
        tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                             <  1:3676004‑M2/36‑4 (MQ=255)
        tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                             >  1:4394492‑M2/1‑33 (MQ=255)
        tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                             >  1:4537539‑M2/1‑33 (MQ=255)
         ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                            <  1:3420464‑M2/36‑5 (MQ=255)
         ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                            >  1:4594175‑M2/1‑32 (MQ=255)
          cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                           >  1:3685538‑M2/1‑31 (MQ=255)
          cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                           >  1:8252873‑M2/1‑31 (MQ=255)
          cAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCg                           <  1:1315894/36‑1 (MQ=255)
           aacccttaagttagcgcttatgggCCGCAACCCCgg                          <  1:2225958‑M2/12‑1 (MQ=255)
           aacccttaagttagcgcttatgggCCGCAACCCCgg                          >  1:4840964‑M2/25‑36 (MQ=255)
           aacccttaagttagcgcttatgggCCGCAACCCCgg                          >  1:3199496‑M2/25‑36 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          >  1:3571826‑M2/1‑30 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          >  1:460962‑M2/1‑30 (MQ=255)
           aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                          <  1:6205484‑M2/36‑7 (MQ=255)
            ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                         >  1:4524361‑M2/1‑29 (MQ=255)
            ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                         >  1:5405410‑M2/1‑29 (MQ=255)
            acccttaagttagcgcttatgggCCGCAACCCCggg                         <  1:3421608‑M2/13‑1 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        <  1:7569420‑M2/14‑1 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                        <  1:1578810‑M2/14‑1 (MQ=255)
              cgcgTTATTTTGCAACCGGGGCCGCAAcccataagc                       >  1:6764757‑M2/1‑27 (MQ=255)
                ttaagttagcgcttatgggCCGCAACCCCGGGtgct                     <  1:344172‑M2/17‑1 (MQ=255)
                ttaagttagcgcttatgggCCGCAACCCCGGGtgct                     >  1:1631082‑M2/20‑36 (MQ=255)
                 gTTATTTTGCAACCGGGGCCGCAAcccataagcgct                    <  1:3474530‑M2/36‑13 (MQ=255)
                 gTTATTTTGCAACCGGGGCCGCAAcccataagcgat                    <  1:1616650‑M2/36‑13 (MQ=255)
                  ttATTTTGCAACCGGGGCCGCAAcccataagcgcta                   >  1:4806779‑M2/1‑23 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                   >  1:4209795‑M2/18‑36 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                   <  1:1880996‑M2/19‑1 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                  >  1:4282991‑M2/1‑22 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAACCCCGGGTgctgct                  <  1:6723916/36‑1 (MQ=255)
                   agttagcgcttatgggCCGCAACCCCGGGTgctgct                  <  1:7336645‑M2/20‑1 (MQ=255)
                   agttagcgcttatgggCCGCAACCCCGGGTgctgct                  <  1:1084723‑M2/20‑1 (MQ=255)
                   agttagcgcttatgggCCGCAACCCCGGGTgctgct                  >  1:1679557‑M2/17‑36 (MQ=255)
                   agttagcgcttatgggCCGCAACCCCGGGTgctgct                  >  1:2782757‑M2/17‑36 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                 <  1:7967347‑M2/21‑1 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                 >  1:7191779‑M2/16‑36 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                 >  1:1685063‑M2/16‑36 (MQ=255)
                     ttttGCAACCGGGGCCGCAAcccataagcgctaact                <  1:4002364‑M2/36‑17 (MQ=255)
                     ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                >  1:2681902‑M2/15‑36 (MQ=255)
                     ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                >  1:3527162‑M2/15‑36 (MQ=255)
                     ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                >  1:4883119‑M2/15‑36 (MQ=255)
                     ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                >  1:1557447‑M2/15‑36 (MQ=255)
                     ttagcgcttatgggACGCAACCCCGGGTGCTGCTaa                >  1:4343052‑M2/15‑36 (MQ=255)
                      tttTCAACCGGGGCCGCAAcccataagcgctaactt               >  1:3711554‑M2/1‑19 (MQ=255)
                      tttGCAACCGGGGCCGCAAcccataagcgctacctt               >  1:2357130‑M2/1‑19 (MQ=255)
                      tttGCAACCGGGGCCGCAAcccataagcgctaactt               <  1:1563108‑M2/36‑18 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt               <  1:1090672‑M2/23‑1 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt               <  1:2524820‑M2/23‑1 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt               <  1:5611100‑M2/23‑1 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt               >  1:2217040‑M2/14‑36 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              <  1:4365296‑M2/36‑19 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta              >  1:481770‑M2/1‑18 (MQ=255)
                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg              >  1:1416796‑M2/13‑36 (MQ=255)
                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg              <  1:1309081‑M2/24‑1 (MQ=255)
                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg              >  1:2636411‑M2/13‑36 (MQ=255)
                        tGCAACCGGGGCCGCAAcccataagcgctaacttaa             >  1:7783020‑M2/1‑17 (MQ=255)
                        gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc             <  1:2025934‑M2/25‑1 (MQ=255)
                         gCAACCGGGGCCGCAAcccataagcgctaacttaag            <  1:605835‑M2/36‑21 (MQ=255)
                          gcttatgggCCGCAACCCCGGGTGCTGTTAATGCgg           >  1:5632960‑M2/10‑36 (MQ=255)
                          cAACCGGGGCCGCAAcccataagcgctaacttaagg           <  1:4119728‑M2/36‑22 (MQ=255)
                           cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa          >  1:2745663‑M2/9‑36 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          <  1:498486‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg          >  1:6794420‑M2/1‑14 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaagat          <  1:7597946‑M2/36‑23 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg        >  1:1379534‑M2/7‑36 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAAAGCGGATg        >  1:5329693‑M2/7‑36 (MQ=255)
                              cGGGGCCGCAAcccataagcgctaacttaagggttg       <  1:72552‑M2/36‑26 (MQ=255)
                              atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc       >  1:3848683‑M2/6‑36 (MQ=255)
                               tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg      <  1:2069515‑M2/32‑1 (MQ=255)
                                gggCCGCAAcccataagcgctaacttaagggttgtg     >  1:3513277‑M2/1‑9 (MQ=255)
                                gggCCGCAACTCGGGGTGCTGCTAATTCGGATGCGa     <  1:1132585/36‑1 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa     <  1:3518265‑M2/33‑1 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa     >  1:3451859‑M2/4‑36 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa     >  1:3623406‑M2/4‑36 (MQ=255)
                                  gCCGCAAcccataagcgctaacttaagggttgtggt   >  1:2011765‑M2/1‑7 (MQ=255)
                                  gCCGCAAcccataagcgctaacttaagggttgtggt   >  1:1409268‑M2/1‑7 (MQ=255)
                                  gCCGCAAcccataagcgctaacttaagggttgtggt   >  1:5358454‑M2/1‑7 (MQ=255)
                                  gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc   >  1:4353371‑M2/2‑36 (MQ=255)
                                  gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc   >  1:782183‑M2/2‑36 (MQ=255)
                                  gCCGCAACCCCGGGTGCTGCTAATGCGCATGCGAcc   >  1:7528593‑M2/2‑36 (MQ=255)
                                   ccGCAAcccataagcgctaacttaagggttgtggta  <  1:4212182‑M2/36‑31 (MQ=255)
                                   ccGCAAcccataagcgctaacttaagggttgtggta  <  1:4647177‑M2/36‑31 (MQ=255)
                                   ccGCAAcccataagcgataacttaagggttgtggta  <  1:3784558‑M2/36‑31 (MQ=255)
                                   |                                   
CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT  >  REL606/4524487‑4524557

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG