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

CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562
                                   |                                   
tggttcaacccttaagttagcgcttatgggCCGCaa                                     <  1:4521772‑M2/6‑1 (MQ=255)
tggttcaacccttaagttagcgcttatgggCCGCaa                                     <  1:1071643‑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)
   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)
      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)
        gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                             >  1:2769149‑M2/1‑28 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             >  1:3619508‑M2/23‑36 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             <  1:4101367‑M2/14‑1 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             <  1:4958451‑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:981573‑M2/1‑26 (MQ=255)
          cttaagttagcgcttatgggCCGCAACCCCGGGtgc                           <  1:2929066‑M2/16‑1 (MQ=255)
          cttaagttagcgcttataggCCGCAACCCCGGGtgc                           <  1:4717369‑M2/16‑1 (MQ=255)
            taagttagcgcttatgggCCGCAACCCCGGGtgctg                         <  1:549587‑M2/18‑1 (MQ=255)
            taagttagcgcttatgggCCGCAACCCCGGGtgctg                         >  1:4935571‑M2/19‑36 (MQ=255)
            taagttaccgcttatgggCCGCAACCCCGGGtgctg                         <  1:5247906‑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: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)
              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)
                 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)
                   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)
                     cAACCGGGGCCGCAAcccataagcgctaacttaagg                <  1:1666676‑M2/36‑22 (MQ=255)
                     cAACCGGGGCCGCAAcccataagcgctaacttaagg                <  1:3339580‑M2/36‑22 (MQ=255)
                      cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa               >  1:5347883‑M2/9‑36 (MQ=255)
                      cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa               <  1:4362459‑M2/28‑1 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:700731‑M2/36‑23 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:73382‑M2/36‑23 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:2048939‑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:5482594‑M2/8‑36 (MQ=255)
                       ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt              <  1:271752‑M2/29‑1 (MQ=255)
                       ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt              >  1:5216527‑M2/8‑36 (MQ=255)
                       ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt              >  1:5102061‑M2/8‑36 (MQ=255)
                       ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt              <  1:6909647‑M2/29‑1 (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)
                           gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa          >  1:7035624‑M2/4‑36 (MQ=255)
                            ggCCGCAAcccataagcgctaacttaagggttgtgg         >  1:1441030‑M2/1‑8 (MQ=255)
                            ggCCGCAACCCCGGGTGCTGCTAATGCGGCTGGGAc         >  1:3056011/1‑36 (MQ=25)
                             gCCGCAAcccataagcgctaacttaagggttgtggt        <  1:5389097‑M2/36‑30 (MQ=255)
                             gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc        >  1:1256085‑M2/2‑36 (MQ=255)
                              ccGCAAcccataagcgctaacttaagggttgtggta       <  1:1860565‑M2/36‑31 (MQ=255)
                              ccGCAAcccataagcgctaacttaagggttgtggta       >  1:7377165‑M2/1‑6 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      <  1:275637‑M2/36‑32 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      <  1:3100521‑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     >  1:5059818/1‑36 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc     >  1:4327599/1‑36 (MQ=255)
                                 cAAcccataagcgctaacttaagggttgtggtatta    >  1:1708353‑M2/1‑3 (MQ=255)
                                 cAAcccataagcgctaacttaagggttgtggtatta    <  1:878329‑M2/36‑34 (MQ=255)
                                 cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa    <  1:6923434/36‑1 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   <  1:1774763‑M2/36‑35 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   <  1:85986/36‑1 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  1:1811709/1‑36 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  1:5540540/1‑36 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  >  1:639954/1‑36 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  <  1:7458128/36‑1 (MQ=255)
                                   |                                   
CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562

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