Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 91.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)30 (0.730) 19/64 NT 85.7% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 5 (0.110)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)35 (0.850) 27/64 NT 97.2% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 1 (0.020)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAC  >  REL606/4524487‑4524555
                                   |                                 
tcagatggttcaacccttaagttagcgcttatgggc                                   >  1:4028586‑M2/36‑36 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggcc                                  >  1:1173358‑M2/35‑36 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggcc                                  <  1:2848554‑M2/2‑1 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggcc                                  <  1:4400817‑M2/2‑1 (MQ=255)
 aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                  >  1:3650443/1‑36 (MQ=255)
  agatggttcaacccttaagttagcgcttatgggCCg                                 <  1:4272475‑M2/3‑1 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                >  1:2240867‑M2/33‑36 (MQ=255)
   gatggttcaacccttaagttagcgcttatgggCCGc                                >  1:1683132‑M2/33‑36 (MQ=255)
     tggttcaacccttaagttagcgcttatgggCCGCaa                              <  1:146696‑M2/6‑1 (MQ=255)
      gTTCCCGGCGCGTTATTTTACAACCGGGGCCGCAAc                             <  1:854934/36‑1 (MQ=37)
       gttcaacccttaagttagcgcttatgggCCGCAAcc                            <  1:2901081‑M2/8‑1 (MQ=255)
       gttcaacccttaagttagcgcttatgggCCGCAAcc                            >  1:3446285‑M2/29‑36 (MQ=255)
        ttcaacccttaagttagcgcttatgggCCGCAAccc                           <  1:3685803‑M2/9‑1 (MQ=255)
        tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                           >  1:1270740‑M2/1‑33 (MQ=255)
        ctcagGCGCGTTATTTTGCAACCGGGGCCGCAAccc                           >  1:492484‑M2/3‑33 (MQ=255)
         ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                          >  1:1907308‑M2/1‑32 (MQ=255)
          caacccttaagttagcgcttatgggCCGCAACCCCg                         <  1:3969227‑M2/11‑1 (MQ=255)
          cAGGCGCCTTATTTTGCAACCGGGGCCGCAAcccat                         >  1:2935418‑M2/1‑31 (MQ=255)
            acccttaagttagcgcttatgggCCGCAACCCCggg                       >  1:1661376‑M2/24‑36 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGCCGCAAcccataag                      >  1:1302775‑M2/1‑28 (MQ=255)
             gcgcGTTATTTTGCAACCGGGGCAGCAAcccataag                      <  1:494805‑M2/36‑9 (MQ=255)
             gcgcGTTATTTTGAAACCGGGGCCGCAAcccataag                      <  1:2895391‑M2/36‑9 (MQ=255)
             cccttaatttagcgcttatgggCCGCAACCCCGGGt                      <  1:5735946‑M2/14‑1 (MQ=255)
             cccttaagttagcgcttatgggCCGCAACCCCGGGt                      <  1:1521643‑M2/14‑1 (MQ=255)
                ttaagttagcgcttatgggCCGCAACCCCGGGtgct                   <  1:1811718‑M2/17‑1 (MQ=255)
                ttaagttagcgcttatgggCCGCAACCCCGGGtgct                   <  1:3687801‑M2/17‑1 (MQ=255)
                ttaagttagcgcttatgggCCGCAACCCCGGGggct                   >  1:5772552‑M2/20‑36 (MQ=255)
                 taagttagcgctaatgggCCGCAACCCCGGGtgctg                  <  1:4334252‑M2/18‑1 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                 >  1:4000619‑M2/18‑36 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                 >  1:3237707‑M2/18‑36 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                 <  1:2697847‑M2/19‑1 (MQ=255)
                  aagttagcgcttatgggCCGCAACCCCGGGtgctgc                 <  1:2438814‑M2/19‑1 (MQ=255)
                   tATTTTGCAACCGGGGCCGCAAcccataagcgccaa                <  1:1912072‑M2/36‑15 (MQ=255)
                   agttagcgcttatgggCCGCAACCCCGGGTgctgct                <  1:1607389‑M2/20‑1 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa               <  1:5530837‑M2/21‑1 (MQ=255)
                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTa               <  1:4565293‑M2/21‑1 (MQ=255)
                    aTTTTGCAACCGGGGCCGCAAcccataagcgctaac               <  1:2124379‑M2/36‑16 (MQ=255)
                    aTTTTGCAACCGGGGCCGCAAcccataagcgctaac               <  1:4349874‑M2/36‑16 (MQ=255)
                      tttGCAGCCGGGGCCGCAAcccataagcgctaactt             >  1:4791943‑M2/1‑19 (MQ=255)
                      tttGCAACCGGGGCCGCAAcccataagcgctaactt             <  1:245067‑M2/36‑18 (MQ=255)
                      tagcgcttatgggCCGCAGCCCCGGGTGCTGCTAAt             <  1:2496134‑M2/23‑1 (MQ=255)
                      tagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt             >  1:1416885‑M2/14‑36 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta            >  1:3205912‑M2/1‑18 (MQ=255)
                       ttGCAACCGGGGCCGCAAcccataagcgctaactta            >  1:3946199‑M2/1‑18 (MQ=255)
                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg            <  1:1790892‑M2/24‑1 (MQ=255)
                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg            >  1:2611123‑M2/13‑36 (MQ=255)
                        tGCTACCGGGGCCGCAAcccataagcgctaacttaa           >  1:1061476‑M2/1‑17 (MQ=255)
                        gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc           <  1:5665567‑M2/25‑1 (MQ=255)
                        cctcttatgggCCGCAACCCCGGGTGCTGCTAATGc           <  1:5775243‑M2/25‑1 (MQ=255)
                         gCAACCGGGGCCGCAAcccataagcgctaacttaag          <  1:4167911‑M2/36‑21 (MQ=255)
                         cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg          <  1:1608590‑M2/26‑1 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg        <  1:939733‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg        <  1:3122010‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg        <  1:1584354‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccataagcgctaacttaaggg        <  1:1414326‑M2/36‑23 (MQ=255)
                           aaCCGGGGCCGCAAcccagaagcgctaacttaaggg        <  1:5076047‑M2/36‑23 (MQ=255)
                            aCCGGGGCCGCAAcccataagcgctaacttaagggt       <  1:3199915‑M2/36‑24 (MQ=255)
                              atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc     <  1:1316326‑M2/31‑1 (MQ=255)
                              atgggCCGCAACCCCGGGCGCTGCTAATGCGGATGc     >  1:4843386‑M2/6‑36 (MQ=255)
                                gggCCGCAAcccataagcgctaacttaagggttgtg   <  1:3356146‑M2/36‑28 (MQ=255)
                                gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa   >  1:1558318‑M2/4‑36 (MQ=255)
                                 ggCCGCAAcccataagcgctaacttaagggttgtgg  >  1:2154019‑M2/1‑8 (MQ=255)
                                   |                                 
CATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAC  >  REL606/4524487‑4524555

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