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 120 total aligned reads displayed.

CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562
                                   |                                   
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:659965‑M2/36‑3 (MQ=255)
  ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                   >  2:122553‑M2/1‑34 (MQ=255)
  gttcaacccttaagttagcgcttatgggCCGCATcc                                   >  2:2746390‑M2/29‑36 (MQ=255)
  gttcaacccttaagttagcgcttatgggCCGCAAcc                                   <  1:1343633‑M2/8‑1 (MQ=255)
   ttcaacccttaagttagcgcttatgggCCGCACccc                                  >  2:439343‑M2/28‑36 (MQ=255)
   tCCAGGCGCGTTATTTTGCAACCGGGGCCGCACccc                                  >  2:1807052‑M2/1‑33 (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)
      aacccttaagttagcgcttatgggCCGCAACCCCgg                               >  2:1121035‑M2/25‑36 (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)
        gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                             <  2:3638539‑M2/36‑9 (MQ=255)
        gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                             >  1:2769149‑M2/1‑28 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             >  1:3619508‑M2/23‑36 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             <  2:362157‑M2/14‑1 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             <  2:1219241‑M2/14‑1 (MQ=255)
         cgcgTTATTTTGCAACCGGGGCCGCAAcccataagc                            >  1:2834966‑M2/1‑27 (MQ=255)
         ccttaagttagcgcttaagggCCGCAACCCCGGGtg                            <  2:3619240‑M2/15‑1 (MQ=255)
          gcgTTATTTTGCAACCGGGGCCGCAAcccataagcg                           >  1:981573‑M2/1‑26 (MQ=255)
          cttaagttagcgcttatgggCCGCAACCCCGGGtgc                           <  1:2929066‑M2/16‑1 (MQ=255)
          cttaagttagcgcttatgggCCGCAACCCCGGGtgc                           <  2:3166362‑M2/16‑1 (MQ=255)
            taagttagcgcttatgggCCGCAACCCCGGGtgctg                         <  1:549587‑M2/18‑1 (MQ=255)
            taagttaccgcttatgggCCGCAACCCCGGGtgctg                         <  2:1508696‑M2/18‑1 (MQ=255)
            taagttaccgcttatgggCCGCAACCCCGGGtgctg                         <  2:3327668‑M2/18‑1 (MQ=255)
            gTTATTTTGCAACCGGGGCCGCAAcccataagcgct                         >  1:1405983‑M2/1‑24 (MQ=255)
             ttATTTTGCAACCGGGGGCGCAAcccataagcgcta                        >  1:2464260‑M2/1‑23 (MQ=255)
             ttATTTTGCAACCGGGGCCGCAAcccataagcgcta                        >  2:3676112‑M2/1‑23 (MQ=255)
             aagttagcgcttatgggCCGCAACCCCGGGtgctgc                        <  1:258081‑M2/19‑1 (MQ=255)
              tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                       <  1:2864217‑M2/36‑15 (MQ=255)
              tATTTTGCAACCGGGGCCGCAAcccataagcgctaa                       >  1:2248073‑M2/1‑22 (MQ=255)
              tATTTAGCAACCGGGGCCGCAAcccataagcgctaa                       >  2:2931582‑M2/1‑22 (MQ=255)
              agttagcgcttatgggCCGCAACCCCGGGTgctgct                       <  1:1823146‑M2/20‑1 (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)
                   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)
                     cAACCGGGGCCGCAAcccataagcgctaacttaagg                <  1:1666676‑M2/36‑22 (MQ=255)
                     cAACCGGGGCCGCAAcccataagcgctaacttaagg                <  1:3339580‑M2/36‑22 (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:316591‑M2/36‑23 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:700731‑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)
                       aCCGGGGCCGCAAcccataagcgctaacttaagggt              >  2:473297‑M2/1‑13 (MQ=255)
                        tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg             <  1:878675‑M2/30‑1 (MQ=255)
                          ggggCCGAAAcccataagcgctaacttaagggttgt           >  2:3208455‑M2/1‑10 (MQ=255)
                           gggCCGCAAcccataagcgctaacttaagggttgtg          >  2:319296‑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)
                             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)
                               cGCAAcccataagcgctaacttaagggttgtggtat      >  1:983902‑M2/1‑5 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      <  2:628678‑M2/36‑32 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      <  1:3100521‑M2/36‑32 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      <  1:275637‑M2/36‑32 (MQ=255)
                               cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt      <  1:1998395/36‑1 (MQ=255)
                                gCAAcccataagcgctaacttcagggttgtggtatt     <  1:3671232‑M2/36‑33 (MQ=255)
                                gCAAcccataagcgctaacttaagggttgtggtatt     <  1:2417509‑M2/36‑33 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc     >  2:1716170/1‑36 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc     >  2:1320608/1‑36 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc     >  2:588389/1‑36 (MQ=255)
                                 cAAcccataagcgctaacttaagggttgtggtatta    <  1:878329‑M2/36‑34 (MQ=255)
                                 cAAcccataagcgctaacttaagggttgtggtatta    >  1:1708353‑M2/1‑3 (MQ=255)
                                 cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa    <  2:3184224/36‑1 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   <  1:1774763‑M2/36‑35 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGTATGCGACCTTCaa   <  2:2191319/36‑1 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   <  1:85986/36‑1 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  1:1811709/1‑36 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  2:1801330/1‑36 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  <  2:3718918/36‑1 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  >  1:639954/1‑36 (MQ=255)
                                   |                                   
CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562

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