Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 84.3% 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 78.4% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 11 (0.180)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)46 (0.810) 31/64 NT 90.2% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 5 (0.080)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 106 total aligned reads displayed.

CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562
                                   |                                   
tggttcaacccttaagttagcgcttatgggCCGCaa                                     >  1:3606800‑M2/31‑36 (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:3566822‑M2/36‑2 (MQ=255)
 gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                    <  1:3306910‑M2/36‑2 (MQ=255)
  gttcaacccttaagttagcgcttatgggCCGCAAcc                                   >  1:1177990‑M2/29‑36 (MQ=255)
  gttcaacccttaagttagcgcttatgggCCGCAAcc                                   >  1:3704833‑M2/29‑36 (MQ=255)
   tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                  >  1:4537539‑M2/1‑33 (MQ=255)
   tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                  >  1:4394492‑M2/1‑33 (MQ=255)
   tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                  <  1:3676004‑M2/36‑4 (MQ=255)
    ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                 <  1:3420464‑M2/36‑5 (MQ=255)
     cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                >  1:8252873‑M2/1‑31 (MQ=255)
     cAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCg                                <  1:1315894/36‑1 (MQ=255)
      aacccttaagttagcgcttatgggCCGCAACCCCgg                               >  1:3199496‑M2/25‑36 (MQ=255)
      aacccttaagttagcgcttatgggCCGCAACCCCgg                               <  1:2225958‑M2/12‑1 (MQ=255)
      aacccttaagttagcgcttatgggCCGCAACCCCgg                               >  1:4840964‑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:5405410‑M2/1‑29 (MQ=255)
       ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                              >  1:4524361‑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:1631082‑M2/20‑36 (MQ=255)
           ttaagttagcgcttatgggCCGCAACCCCGGGtgct                          <  1:344172‑M2/17‑1 (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:1679557‑M2/17‑36 (MQ=255)
              agttagcgcttatgggCCGCAACCCCGGGTgctgct                       <  1:1084723‑M2/20‑1 (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)
                ttttGCAACCGGGGCCGCAAcccataagcgctaact                     <  1:4002364‑M2/36‑17 (MQ=255)
                ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                     >  1:1557447‑M2/15‑36 (MQ=255)
                ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                     >  1:4883119‑M2/15‑36 (MQ=255)
                ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                     >  1:3527162‑M2/15‑36 (MQ=255)
                ttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                     >  1:2681902‑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:5611100‑M2/23‑1 (MQ=255)
                 tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                    >  1:2217040‑M2/14‑36 (MQ=255)
                 tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                    <  1:2524820‑M2/23‑1 (MQ=255)
                  ttGCAACCGGGGCCGCAAcccataagcgctaactta                   >  1:481770‑M2/1‑18 (MQ=255)
                  ttGCAACCGGGGCCGCAAcccataagcgctaactta                   <  1:4365296‑M2/36‑19 (MQ=255)
                  agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                   >  1:1416796‑M2/13‑36 (MQ=255)
                  agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                   >  1:2636411‑M2/13‑36 (MQ=255)
                  agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                   <  1:1309081‑M2/24‑1 (MQ=255)
                   tGCAACCGGGGCCGCAAcccataagcgctaacttaa                  >  1:7783020‑M2/1‑17 (MQ=255)
                   gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                  <  1:2025934‑M2/25‑1 (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:6794420‑M2/1‑14 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:498486‑M2/36‑23 (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:782183‑M2/2‑36 (MQ=255)
                             gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc        >  1:4353371‑M2/2‑36 (MQ=255)
                             gCCGCAACCCCGGGTGCTGCTAATGCGCATGCGAcc        >  1:7528593‑M2/2‑36 (MQ=255)
                              ccGCAAcccataagcgctaacttaagggttgtggta       <  1:4647177‑M2/36‑31 (MQ=255)
                              ccGCAAcccataagcgctaacttaagggttgtggta       <  1:4212182‑M2/36‑31 (MQ=255)
                              ccGCAAcccataagcgataacttaagggttgtggta       <  1:3784558‑M2/36‑31 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      >  1:6013391‑M2/1‑5 (MQ=255)
                               cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt      >  1:1895251/1‑36 (MQ=255)
                               cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt      <  1:5017486/36‑1 (MQ=255)
                                gCAAcccataagcgctaacttaagggttgtggtatt     >  1:4523484‑M2/1‑4 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGCTGCGACCTTc     <  1:5344929/36‑1 (MQ=255)
                                 cAAcccataagcgctaacttaagggttgtggtatta    <  1:2407216‑M2/36‑34 (MQ=255)
                                 cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa    >  1:2499131/1‑36 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   <  1:7565401‑M2/36‑35 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   <  1:5612767‑M2/36‑35 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGTATGCGACCTTCaa   >  1:425164/1‑36 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  1:1578575/1‑36 (MQ=255)
                                   acccataagcgctaacttaagggttgtggtattacg  <  1:5356340‑M2/36‑36 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  <  1:2251997/36‑1 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  >  1:328205/1‑36 (MQ=255)
                                   |                                   
CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562

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