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)

AGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG...................................  <  REL606/15423‑15386
...................................GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAAC  <  REL606/4524527‑4524493
                                                                         
   CGTAATACCACATCCCTTAAGTTAGCGCTTATGGGT                                    <  1:3574979/36‑1 (MQ=0)
   CGTAATACCACAACCCTTAAGTTAGCGCTTATGGGT                                    <  1:1040859/36‑1 (MQ=1)
    GTAATACCACAACCCTTAAGTTAGCGCTTATGGGTT                                   <  1:7551936/36‑1 (MQ=1)
     TAATACCACAACCCTTAAGTTAGCGCTTATGGGTTG                                  >  1:2905070/1‑36 (MQ=1)
      AATACCACAACCCTTAAGTTAGCGCTTATGGGTTGC                                 <  1:5168277/36‑1 (MQ=1)
            ACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCC                           <  1:5018394/36‑1 (MQ=1)
                CCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTT                       <  1:5732424/36‑1 (MQ=1)
                 CCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTG                      >  1:650778/1‑36 (MQ=1)
                     AAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAA                  >  1:233194/1‑36 (MQ=1)
                     AAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAA                  >  1:7105572/1‑36 (MQ=1)
                      AGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAA                 >  1:3026884/1‑36 (MQ=1)
                       GTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAAT                <  1:1277342/36‑1 (MQ=1)
                       GTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAAT                >  1:5267936/1‑36 (MQ=1)
                           GCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACG            >  1:1004166/1‑36 (MQ=1)
                              CTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGC         <  1:4529168/36‑1 (MQ=1)
                                 ATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTG      <  1:2359398/36‑1 (MQ=1)
                                                                         
AGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG...................................  <  REL606/15423‑15386
...................................GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAAC  <  REL606/4524527‑4524493

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG