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

CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562
                                   |                                   
tggttcaacccttaagttagcgcttatgggCCGCaa                                     <  1:146696‑M2/6‑1 (MQ=255)
 gTTCCCGGCGCGTTATTTTACAACCGGGGCCGCAAc                                    <  1:854934/36‑1 (MQ=37)
  gttcaacccttaagttagcgcttatgggCCGCAAcc                                   >  1:3446285‑M2/29‑36 (MQ=255)
  gttcaacccttaagttagcgcttatgggCCGCAAcc                                   <  1:2901081‑M2/8‑1 (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:3687801‑M2/17‑1 (MQ=255)
           ttaagttagcgcttatgggCCGCAACCCCGGGtgct                          <  1:1811718‑M2/17‑1 (MQ=255)
           ttaagttagcgcttatgggCCGCAACCCCGGGggct                          >  1:5772552‑M2/20‑36 (MQ=255)
            taagttagcgctaatgggCCGCAACCCCGGGtgctg                         <  1:4334252‑M2/18‑1 (MQ=255)
             aagttagcgcttatgggCCGCAACCCCGGGtgctgc                        <  1:2438814‑M2/19‑1 (MQ=255)
             aagttagcgcttatgggCCGCAACCCCGGGtgctgc                        >  1:3237707‑M2/18‑36 (MQ=255)
             aagttagcgcttatgggCCGCAACCCCGGGtgctgc                        >  1:4000619‑M2/18‑36 (MQ=255)
             aagttagcgcttatgggCCGCAACCCCGGGtgctgc                        <  1:2697847‑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:3946199‑M2/1‑18 (MQ=255)
                  ttGCAACCGGGGCCGCAAcccataagcgctaactta                   >  1:3205912‑M2/1‑18 (MQ=255)
                  agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                   >  1:2611123‑M2/13‑36 (MQ=255)
                  agcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                   <  1:1790892‑M2/24‑1 (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:1414326‑M2/36‑23 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:1584354‑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)
                              cggcaAcccataagcgctaacttaagggttgtggta       >  1:4263118‑M2/3‑6 (MQ=255)
                               cGCAAcccataagcgctaacttaagggttgtggtat      <  1:2673679‑M2/36‑32 (MQ=255)
                               cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt      <  1:3997067/36‑1 (MQ=255)
                               cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt      <  1:3984827/36‑1 (MQ=255)
                                 cAAcacataagcgctaacttaagggttgtggtatta    >  1:4862711‑M2/1‑3 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   <  1:1758086‑M2/36‑35 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   >  1:4699319‑M2/1‑2 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggaattac   <  1:64882‑M2/36‑35 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   <  1:2418071/36‑1 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  1:5215706/1‑36 (MQ=255)
                                  aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa   >  1:3505549/1‑36 (MQ=255)
                                   acccataagcgctaacttaagggttgtggtattacg  <  1:2337785‑M2/36‑36 (MQ=255)
                                   acccataagcgctaacttaagggttgtggtattacg  <  1:616437‑M2/36‑36 (MQ=255)
                                   |                                   
CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562

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