Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 68.8% 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)16 (0.290) 14/64 NT 61.5% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 10 (0.230)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)28 (0.500) 20/64 NT 73.7% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 10 (0.250)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAC  >  REL606/4524487‑4524555
                                   |                                 
tcagatggttcaacccttaagttagcgcttatgggc                                   <  1:1773295‑M2/1‑1 (MQ=255)
cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                   <  1:6741431/36‑1 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggcc                                  <  1:7427683‑M2/2‑1 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                  <  1:6736530/36‑1 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                  >  1:4249733/1‑36 (MQ=255)
  ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                 >  1:3076360/1‑36 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                 >  1:521101‑M2/34‑36 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                 >  1:2712266‑M2/34‑36 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                 >  1:2472392‑M2/34‑36 (MQ=255)
   tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                <  1:5265971/36‑1 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                >  1:3356954‑M2/33‑36 (MQ=255)
    ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                               >  1:6451797/1‑36 (MQ=255)
    atggttcaacccttaagttagcgctgatgggCCGCa                               >  1:7647249‑M2/32‑36 (MQ=255)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                             <  1:1535038‑M2/36‑2 (MQ=255)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                             <  1:3848795‑M2/36‑2 (MQ=255)
      gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                             <  1:3284932‑M2/36‑2 (MQ=255)
        ttcaacccttaagttagcgcttatgggCCGCAAccc                           <  1:3294312‑M2/9‑1 (MQ=255)
         ttaacccttaagttagcgcttatgggCCGCAAcccc                          >  1:6835638‑M2/27‑36 (MQ=255)
         tcaacccttaagttagcgcttatgggCCGCAAcccc                          <  1:1223911‑M2/10‑1 (MQ=255)
         tcaacccttaagttagcgcttatgggCCGCAAcccc                          <  1:490759‑M2/10‑1 (MQ=255)
          caacccttaagttagcgcttatgggCCGCAACCCCg                         <  1:7235455‑M2/11‑1 (MQ=255)
          cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                         <  1:2359398‑M2/36‑6 (MQ=255)
           aacccttaagttagcgcttatgggCCGCAACCCCgg                        >  1:2776955‑M2/25‑36 (MQ=255)
           aacccttaagttagcgcttatgggCCGCAACCCCgg                        >  1:3652032‑M2/25‑36 (MQ=255)
            acccttaagttagcgcttatgggCCGCAACCCCggg                       <  1:5619271‑M2/13‑1 (MQ=255)
            acccttaagttagcgcttatgggCCGCAACCCCggg                       <  1:3000885‑M2/13‑1 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                      <  1:4529168‑M2/36‑9 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                      <  1:4500986‑M2/14‑1 (MQ=255)
                cgTTATTTTGCAACCGGGGCCGCAAcccataagcgc                   >  1:1004166‑M2/1‑25 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                  <  1:4457500‑M2/18‑1 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                  >  1:1650628‑M2/19‑36 (MQ=255)
                 taagttagcgcttatgggCCGCAACCCCGGGtgctg                  >  1:3288968‑M2/19‑36 (MQ=255)
                  ttATTTTGCAACCGGGGCCGCAACCCCGGGtgctgc                 >  1:275870/1‑36 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAACCCCGGGTgctgct                <  1:383289/36‑1 (MQ=255)
                   agttagcgcttatgggCCGCAACCCCGGGTgctgct                <  1:750232‑M2/20‑1 (MQ=255)
                    aTTTTGCAACCGGGGCCGCAAcccataagcgctaac               <  1:1277342‑M2/36‑16 (MQ=255)
                    aTTTTGCAACCGGGGCCGCAAcccataagcgctaac               >  1:5267936‑M2/1‑21 (MQ=255)
                     ttttGCAACCGGGGCCGCAAcccataagcgctaact              >  1:3026884‑M2/1‑20 (MQ=255)
                     ttttGCAACCGGGGCCGCAACCCCGGGNNCTGCTaa              >  1:7237176/1‑36 (MQ=255)
                      tttGCAACCGGGGCCGCAAcccataagcgctaactt             >  1:7105572‑M2/1‑19 (MQ=255)
                      tttGCAACCGGGGCCGCAAcccataagcgctaactt             >  1:233194‑M2/1‑19 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGTTGCTGCTAAt             <  1:153577‑M2/23‑1 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt             <  1:6135238‑M2/23‑1 (MQ=255)
                       ttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATg            <  1:3646407/36‑1 (MQ=255)
                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg            <  1:1802188‑M2/24‑1 (MQ=255)
                          cAACCGGGGCCGCAAcccataagcgctaacttaagg         >  1:650778‑M2/1‑15 (MQ=255)
                           cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa        <  1:7318405‑M2/28‑1 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg        <  1:5732424‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGa        >  1:7045446/1‑36 (MQ=255)
                            ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt       <  1:1986219‑M2/29‑1 (MQ=255)
                             tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg      <  1:5794824‑M2/30‑1 (MQ=255)
                             ccGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATg      <  1:5291096/36‑1 (MQ=255)
                             ccGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATg      >  1:4395226/1‑36 (MQ=255)
                              cGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGc     <  1:1615924/36‑1 (MQ=255)
                               tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg    <  1:7084599‑M2/32‑1 (MQ=255)
                               ggggCCGCAAcccataagcgctaacttaagggttgt    <  1:5018394‑M2/36‑27 (MQ=255)
                               ggggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg    >  1:3120911/1‑36 (MQ=255)
                               ggggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg    <  1:3404155/36‑1 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa   >  1:4732735‑M2/4‑36 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa   >  1:4945899‑M2/4‑36 (MQ=255)
                                 ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc  >  1:3569619‑M2/3‑36 (MQ=255)
                                 ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc  >  1:1464678‑M2/3‑36 (MQ=255)
                                   |                                 
CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAC  >  REL606/4524487‑4524555

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